{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "480765be",
   "metadata": {},
   "outputs": [],
   "source": [
    "from astropy.io import fits\n",
    "import HAWCtoHerschel as HtoH"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "d5477674",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Source information\n",
    "SourceName = 'Fil1'\n",
    "FullName = 'Filament 1'\n",
    "# Loading the fits files\n",
    "SourceT = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Tdust_10_10.fits')\n",
    "SourceN = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_N_10_10.fits')\n",
    "SourcePol = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Gal.fits')\n",
    "SourcePlanck = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Planck.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "8fbf97a2",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Plot parameters\n",
    "ScaleLength = 0.0136419 # 0.0139746 # 1 pc in degrees for the target\n",
    "\n",
    "# Region location, width, and height\n",
    "l_center = 26.9076450 # Position in galactic longitude (degrees)\n",
    "b_center = -0.3039280 # Position in galactic latitude (degrees)\n",
    "l_width = 0.2645292 # Width in galactic longitude (degrees)\n",
    "b_height = 0.1365312 # Height in galactic latitude (degrees)\n",
    "Region = [l_center,b_center,l_width,b_height]\n",
    "\n",
    "# Intensity contours for the plot\n",
    "Icontours = [10, 20, 30, 40, 50, 60] # in mJy/arcsec^2\n",
    "# Ranges for plots\n",
    "Nscale = [0,1.2e+22] # Range for column density in cm^-2\n",
    "Tscale = [15.0,30.0] # Range for temperature map in K\n",
    "Iscale = [0.0,20.0] # Range for intensity map in mJy/arcsec^2\n",
    "# Vector scaling fraction in polarization plots\n",
    "Pvscale = 0.6\n",
    "\n",
    "# Figure size\n",
    "MapFigSize = [10.67,6] # in inches\n",
    "\n",
    "# Effective HAWC+ beam size to use\n",
    "HAWCBeam = 0.00519 # 18.7'' in Band E"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "a2a163b6",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Beam Size of polarization data\n",
      "18.2\n",
      "Beam size of target pojection\n",
      "36.395999999999994\n",
      "Pixel size of the polarization data\n",
      "4.550000000000001\n",
      "Pixel size of the target projection\n",
      "14.0000004\n",
      "Target standard deviation for smoothing in arcseconds\n",
      "13.384761981071023\n",
      "Target standard deviation for smoothing in pixels\n",
      "2.941705929905719\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: nan_treatment='interpolate', however, NaN values detected post convolution. A contiguous region of NaN values, larger than the kernel size, are present in the input array. Increase the kernel size to avoid this. [astropy.convolution.convolve]\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:96: RuntimeWarning: invalid value encountered in divide\n",
      "  dPI = ((StokesQ*dQ)**2.0 + (StokesU*dU)**2.0)**0.5/PI_biased\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:106: RuntimeWarning: invalid value encountered in divide\n",
      "  P_biased = 100.0*PI_biased/StokesI\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:109: RuntimeWarning: invalid value encountered in divide\n",
      "  dP = np.absolute(P_biased*((dPI/PI_biased)**2.0 + (dI/StokesI)**2.0)**0.5)\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:120: RuntimeWarning: invalid value encountered in divide\n",
      "  dO = (0.5*180.0/math.pi)*((StokesQ*dU)**2.0 +\n"
     ]
    }
   ],
   "source": [
    "# Reprojecting the HAWC+ Stokes data to the Herschel pixel scale\n",
    "NewPol = HtoH.ReprojectCube(SourcePol,SourceN, Smooth=True)\n",
    "#NewPol.writeto('test2.fits', overwrite=True)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "31eb9087",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:228: RuntimeWarning: invalid value encountered in divide\n",
      "  imask_01 = np.where(Pdata[0].data/Pdata[1].data < IdI) # Total intensity SNR threshold\n"
     ]
    }
   ],
   "source": [
    "# Creating downsampled vector catalogs (independent vectors)\n",
    "# Creating the mask for the region\n",
    "DownsampleMask = HtoH.DownsampleVectors(NewPol,Step=2)\n",
    "# Creating the new catalogs for the region\n",
    "CatI, CatdI, CatQ, CatdQ, CatU, CatdU, CatP, CatdP, CatO, CatB,CatdO, CatPI, CatdPI, CatN, CatT = HtoH.MakePCats(SourceN, SourceT, SourcePol, NewPol, CatMask=DownsampleMask)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "2c4fb0cb",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "===\n",
      "Curve fit test\n",
      "Log scale fit: (-0.840020379906213,25.623764523064075)\n",
      "Direct power law fit: (-0.7580390821914124,22.82043835759203)\n",
      "===\n",
      "\n",
      "Power law index: -0.840020379906213 ± 0.0451406865958307\n",
      "Coefficient: 25.623764523064075 ± 2.4984965492426383\n",
      "Chi-Squared: 172.37677611253156\n",
      "Number of elements: 123\n",
      "Reduced Chi-Squared: 1.4246014554754673\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAegAAAF2CAYAAAClPAs/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABuLElEQVR4nO3deVhU1RvA8e8Mq8iigAgouOe+b5mVlqaZv1wys1VEtDRNjVyyLHPPFjWXtFwwS3OpbM8scylzz31fUHADFQFxQWDO749pJoZ1ZhiYAd7P89zH4c6957734vDOufcsGqWUQgghhBAORWvvAIQQQgiRnSRoIYQQwgFJghZCCCEckCRoIYQQwgFJghZCCCEckCRoIYQQwgFJghZCCCEckCRoIYQQwgFJghZCCCEckCToAkpNTeWdd94hNTXV3qEIITKRz6Yo7hwiQc+fP59GjRrh7e2Nt7c3bdq04ZdffjG+f+fOHYYMGYKfnx+enp706tWLuLg4O0b8n9TUVCZMmCB/BIRwMPLZFMWdQyToypUr8+6777Jnzx52797Nww8/TPfu3Tl8+DAAr776Kj/88ANr1qxh8+bNXLx4kSeeeMLi48ybN69QtxemisP1K+oYC/N4tiq7oOVYs798NotWcbh+9oixsI5pdbnKQZUvX14tWrRIJSYmKhcXF7VmzRrje0ePHlWA2rZtm0Vl1q1b1+bbJyUlKUAlJSVZVHZpYOn1toeijrEwj2ersgtajjX7y2ezaMlns2iPaW25zjb8kmATGRkZrFmzhps3b9KmTRv27NlDWloaHTt2NG5Tp04dQkND2bZtG/fee2+O5aSmpprc2tLpdKSmppKUlIRGozE7luTk5Dy3Mbyf33alkTnXz96KOsbCPJ6tyi5oOdbsb+k+8tksGPlsFt0xlVKkpqaSmJiIVvvfTWs3Nzfc3Nzy3dkhHDhwQJUtW1Y5OTkpHx8f9dNPPymllFq+fLlydXXNtn3Lli3V6NGjcy1v/PjxCpBFFllkkUUWh1vGjx+fb150mBp07dq12bdvH0lJSXz11VeEhYWxefNmq8sbO3YskZGRxp+TkpIIDQ0lNjYWb09PaNYMoqNh9mwIC7PFKQghhBAmkpOTCQkJISYmBh8fH+P6fGvPgMMkaFdXV2rWrAlA8+bN2bVrFx999BF9+vTh7t27JCYmUq5cOeP2cXFxBAYG5lpebrcPDC3FGTIERo6EJUtg6FAw87a3EEIIYSkfHx997rGAQ7TizonhmXHz5s1xcXFhw4YNxveOHz9OTEwMbdq0sf4A/fqBuzvs2wc7dhQ4XiGEEMKWHKIGPXbsWLp06UJoaCg3btxgxYoVbNq0iV9//RUfHx8iIiKIjIzE19cXb29vXnnlFdq0aZNrAzGz+PnB00/D0qUwbx4UpCwhhBDCxhyiBh0fH0/fvn2pXbs2HTp0YNeuXfz666888sgjAMycOZP//e9/9OrViwcffJDAwEC++eabgh/45Zf1/65eDVeuFLw8IYQQwkY0Sill7yCKQnJyMj4+PiQlJZk+B2jVCnbtgnffhTFj7BegEEKIEifX3GMGh6hB25WhFr1gAWRk2DcWIYQQ4l+SoPv0gfLl4exZyDT+txBCCGFPkqDLlIH+/fWvP/7YvrEIIYQQ/5IEDTBokP7fdevg9Gn7xiKEEEIgCVqvZk3o3BmUgk8+sXc0QgghhCRooyFD9P8uXgy3b9s3FiGEEKWeJGiDxx6D0FBISND3ixZCCCHsSBK0gZPTf8+ipbGYEEIIO5MEnVlEBLi4wM6dsHu3vaMRQghRikmCziwgAHr31r+WWrQQQgg7kgSdlWFksS+/1D+PFkIIIexAEnRW990HjRvDnTv6ma6EEEIIO5AEnZVG818teu5cSE+3bzxCCCFKJUnQOXn+ef180dHRYItpLYUQQggLSYLOiYfHfwOXfPCBfoQxIYQQoghJgs7NkCHg7q6fK3rLFntHI4QQopSRBJ2bgAAIC9O//uAD+8YihBCi1JEEnZfISH2jsR9/hCNH7B2NEEIUyNKlS9FoNGg0GuO69u3bo9Fo6Nevn/0CEzmSBJ2Xe+6B7t31rz/80L6xCCFEHgyJNqfl22+/BaBChQq0bt2a1q1b2zdYG+rXrx8ajYb27dubtf0rr7xC48aNcXZ2RqPREBgYWLgBFoCzvQNweKNGwbffwhdfwOTJEBRk74iEECJXrq6uNG3a1GSdr68vAF27dqVr1672CMthfP7557i6uuLr68uVK1fsHU6epAadn/vu0y9378KcOfaORggh8hQUFMT27dtNlgcffBDI+RZ3TsLCwqhVqxZeXl64urpSpUoVhg0bRnJysnEbQ821atWqfPbZZ4SEhODt7c2IESO4ffs2I0aMwMfHh9DQUBYsWGBS/sWLF+nfvz/BwcG4urpSvXp1Jk2aRHqmcScMdwT69u3L+PHjCQoKonz58jz//PPcuHEDwHhsgM2bNxvPbdOmTbme28GDB4mPj+exxx6z6LragyRoc4wcqf93/nxISbFvLEIIUci+++47rl+/To0aNQgJCSEmJoY5c+YQERGRbdtLly7x8ssv4+bmxo0bN/joo49o3rw5UVFR+Pj4EBsby5AhQzh69CgA165d49577yUqKoqUlBTq1q1LbGwsb7/9Ni+++GK28leuXMnMmTMpU6YMiYmJLF++nHfffReApk2b4u/vD4CXl5fx9r23t3eu5xYSEmKLS1QkJEGbo1s3qFULEhNh8WJ7RyOEELk6d+5ctmfQltq8eTNXr15l3759nD59mjfffBOAb7/9ljt37phse/fuXdavX8+JEyeMye/UqVPs37+fI0eO4O7ujk6nY/PmzQDMnTuX2NhYKlasyOnTp9m/fz9fffUVoK/hnzp1yqR8d3d3jh49yqlTp2jevDkAGzZsAGDt2rXGW/bNmjUz3jFo1qyZxefsiCRBm8PJCV57Tf965kwZ/lMI4bBcXV2NNUlrG4T9/vvvNGjQgDJlyqDRaJgyZQoA6enp2Z7bli9fnrZt26LVagkNDQWgQYMGVK1aFU9PTwICAgCIi4sDYOfOncafAwIC0Gg09OjRAwClFDt27DAp/+GHH6ZSpUpotVrq1KljUlZJJ43EzNW3L7z1Fpw7B199BU8/be+IhBAiG8MzaGstX76ckf8+1gsKCiIkJISrV69y5swZADIyMky2z3w72dnZOds6Qw1eZRmR0cvLi3r16mU7voeHh8nP5cqVy1Z+1rJKKqlBm6tMGRg6VP/6/fdl+E8hRIlkSO5eXl5ER0ezY8cOOnXqZLPyW7ZsCeiT7cqVK423pX/77TdefvllevbsaVF5hoR+8+ZNm8XoKCRBW+Lll/WJ+p9/II9WgkIIUVw1atQIgBs3blC9enWqV6/O6tWrbVb+kCFDqFSpEtevX6d27do0adKEGjVq4OfnR5hh9EYLGG577969m4YNG3Lvvfdy+/btXLdv3749NWvW5Jt/J0K6evUqNWvWpGbNmtlur9ubJGhL+PtDeLj+9fvv2zcWIYQoBBEREURGRuLv78+NGzdo3749EydOtFn5FSpUYPv27YSHh+Pn58fhw4e5ffs2DzzwADNnzrS4vP79+9OrVy98fHw4dOgQO3bsyHYbPrOzZ89y+vRpY1etjIwMTp8+zenTp/NM7PagUaXkZn5ycjI+Pj4kJSXl2QQ/X6dP60cY0+ng4EFo0MB2QQohhChRCpJ7pAZtqRo14Ikn9K9l+E8hhBCFRBK0NQwDlyxfDhcu2DcWIYQQJZIkaGu0bg0PPABpaTB7tr2jEUIIUQJJgrbWqFH6fxcsgEzj0wohhBC2IAnaWl27Qp06+uS8aJG9oxFCCFHCSIK2llb73/Cfs2bpb3cLIYQQNiIJuiCefx4qVoTYWPjyS3tHI4QQogSRBF0Q7u7w6qv615MnyyQaQgghbEYSdEENGQJ+fnDyJKxcae9ohBBClBAOkaCnTZtGy5Yt8fLyIiAggB49enD8+HGTbdq3b59tjtNBgwbZKeJMPD3/exY9eTLkMcScEEIUtXnz5lG1alXc3d1p3bq1cbrH3GRkZPDWW29RrVo1ypQpQ40aNZg0aVKOM0i9++67aDQaRowY4RCxm7PPO++8ky2XGMbzdjjKAXTu3FlFRUWpQ4cOqX379qnHHntMhYaGqpSUFOM27dq1UwMHDlSXLl0yLklJSWYfIykpSQEW7WO25GSlfH2VAqVWrLB9+UIIYYWVK1cqV1dXtWTJEnX48GE1cOBAVa5cORUXF5frPlOmTFF+fn7qxx9/VNHR0WrNmjXK09NTffTRRybb7dy5U1WtWlU1atRIDR8+3OyY2rVrp6KiogoldnP2GT9+vKpfv75JLrly5YrZ8VuqILnHIRJ0VvHx8QpQmzdvNq5r166dRf8JsirUBK2UUpMn6xN03bpKpacXzjGEEA5hx44dql27dsrd3V3Vrl1b7dq1S33yySfq8ccft3doJlq1aqWGDBli/DkjI0MFBweradOm5bpP165dVf/+/U3WPfHEE+q5554z/nzjxg1Vq1Yt9dtvv1n8t9ncBG1N7ObsM378eNW4cWOz4y2oguQeh7jFnVVSUhIAvr6+JuuXL1+Ov78/DRo0YOzYsdy6dSvXMlJTU0lOTjZZCtUrr0C5cnD0KHz1VeEeSwhhN9u3b6ddu3Z07dqVAwcOULduXSZOnMj06dOZMGGCybZTp07F09MzzyUmJqZQ4rx79y579uyhY8eOxnVarZaOHTuybdu2XPe777772LBhAydOnABg//79/PXXX3Tp0sW4zZAhQ+jatatJ2faO3ZJ9Tp48SXBwMNWrV+e5554rtN9BZlnzUWpqav47FcIXhgLJyMhQXbt2VW3btjVZ/8knn6h169apAwcOqC+++EJVqlRJ9ezZM9dyxo8fr4BsS6HVoJVSasIEfS26fn2lMjIK7zhCCLtp06aNeuGFF4w/r1q1Smm12hz/Hl27dk2dPHkyzyUtLa1Q4rxw4YIC1N9//22yftSoUapVq1a57peRkaHGjBmjNBqNcnZ2VhqNRk2dOtX4/pdffqkaNGigbt++rZTK/+7mlClTVNmyZY2LVqtVbm5uJuvOnTtX4NjN3efnn39Wq1evVvv371fr1q1Tbdq0UaGhoSo5OTnXcygIQw066zJ+/Ph893UulK8KBTBkyBAOHTrEX3/9ZbL+xRdfNL5u2LAhQUFBdOjQgdOnT1OjRo1s5YwdO5bIyEjjz8nJyYSEhBRe4ADDhsGMGXD4MHz9NfTuXbjHE0IUqfPnz7Nt2zY++OAD4zpnZ2eUUtlqz6C/C5j1TqClXn/9daZPn57nNkePHrVZQ6fVq1ezfPlyVqxYQf369dm3bx8jRowgODiYhx9+mOHDh/Pbb7/h7u5uVnmDBg3iqaeeMv783HPP0atXL54wzAoIBAcH2yR2c2S+E9CoUSNat25NlSpVWL16NREREYV23NjYWJPpJt3c3PLdx6ES9NChQ/nxxx/ZsmULlStXznPb1q1bA3Dq1KkcE7Sbm5tZF8CmypWDESNgwgSYOBF69dKPOCaEKBGOHj0KQLNmzYzrjh8/TqtWrWjYsGG27adOncrUqVPzLPPIkSOEhobm+v5rr71Gv3798iyjevXq2db5+/vj5OREXFycyfq4uDgCAwNzLWvUqFG8/vrrPP3004C+QnTu3DmmTZuGj48P8fHxJuefkZHBli1bmDt3LqmpqTg5OZmUl/VLSpkyZQgICKBmzZq5xmBN7Naeb7ly5bjnnns4depUrtvYgre3t8XzQTtEglZK8corr7B27Vo2bdpEtWrV8t1n3759AAQFBRVydBYaPhxmzoRDh2DtWn2SFkKUCElJSTg5OaHRaABISEjggw8+oHHjxjlun7X2mJP8ao8VKlSgQoUKFsfq6upK8+bN2bBhAz169ABAp9OxYcMGhg4dmut+t27dQpulYuHk5IROp6NDhw4cPHjQ5L3w8HDq1KnDmDFjsiVna1kTu7Xnm5KSwunTp3nhhRdsErtNFcpNdwsNHjxY+fj4qE2bNpk0fb9165ZSSqlTp06piRMnqt27d6vo6Gj13XffqerVq6sHH3zQ7GMUeivuzN56S/8sulEjeRYtRAly8uRJBaiJEyeqo0ePqk6dOqlmzZqpihUrqrNnz9o7vGxWrlyp3Nzc1NKlS9WRI0fUiy++qMqVK6cuX75s3GbOnDnq4YcfNv4cFhamKlWqZOxm9c033yh/f381evToHI+R3zPoGzdumPxdz2lJz6HnizWxm7PPa6+9pjZt2qSio6PV1q1bVceOHZW/v7+Kj48365paqth3syKHB+iAsSl+TEyMevDBB5Wvr69yc3NTNWvWVKNGjXKcftBZXbumlJeXPkl/803hH08IUWQmTpyo/Pz8lLu7u+rXr5+6evWqatasmapTp469Q8vRnDlzVGhoqHJ1dVWtWrVS27dvN3l//PjxqkqVKsafk5OT1fDhw1VoaKhyd3dX1atXV2+++aZKTU3Nsfz8EnRuDXYzL9HR0TaJ3Zx9+vTpo4KCgpSrq6uqVKmS6tOnjzp16lSu8RdUQXKPRqkchocpgZKTk/Hx8SEpKcni5wBWGTcOpkyBJk3gn3/g31tiQgghSo+C5B5pwVRYXn1VPwzovn3w/ff2jkYIIUQxU6AEnZaWRmxsLMePHychIcFWMZUMfn76wUtA36q7dNyoEEIIYSMWJ+gbN24wf/582rVrh7e3N1WrVqVu3bpUqFCBKlWqMHDgQHbt2lUYsRY/kZFQtizs3Qs//mjvaIQQQhQjFiXoGTNmULVqVaKioujYsSPffvst+/bt48SJE2zbto3x48eTnp5Op06dePTRRzl58mRhxV08+Pvrp6MEqUULIYSwiEWNxJ555hnGjRtH/fr189wuNTWVqKgoXF1d6d+/f4GDtIUibyRmcOUKVK0Kt27pa9FduxbdsYUQQthVkTUS+/LLL/NNzqAfxWvQoEEOk5ztqkIFePll/WupRQshLLB3715cXFxo3759oZQ/bdo0WrZsiZeXFwEBAfTo0YPjx4+bbGPt/MkXLlzg+eefx8/PjzJlytCwYUN2795d5PFnVZzmg7ZJK+6rV6/y008/8f3333Pp0iVbFFmyjBqlfxa9a5d+dDEhhDDDsGHDGDlyJPv37y+U8jdv3syQIUPYvn07v/32G2lpaXTq1ImbN2+abFe/fn0uXbpkXLLOlZDV9evXadu2LS4uLvzyyy8cOXKEDz/8kPLly+cbU/v27Vm6dKlN48/K0vOxlwIP9fn1118TERHBPffcQ1paGsePH2fevHmEh4fbIr6SISBA3+1q8mR44w3o1g2cHWKUVSGEg1qxYgXly5dnyJAhvPvuu5w9e5aqVava9Bjr1q0z+Xnp0qUEBASwZ88eHnzwQeN6Z2fnPMezzmr69OmEhIQQFRVlXGfOEM6WMjf+rCw9H3uxuAadkpJi8vOECRPYuXMnO3fuZO/evaxZs4Y333zTZgGWGKNG6bteHT8OZn47FEKUTjdv3uSNN95g+vTpVK5cGR8fH+P8A1nZcs7ppKQkgGwzcFk6f/L3339PixYt6N27NwEBATRt2pSFCxeaFUNB5BZ/VvaYD9oaFifo5s2b89133xl/dnZ2Jj4+3vhzXFwcrq6utomuJPH2BsMXl3fegdu37RqOEMJxTZ06lUcffZS6desCUK9evVwT9KBBg9i3b1+eiznTOep0OkaMGEHbtm1p0KCBcX3r1q1ZunQp69atY/78+URHR/PAAw9w48aNXMs6c+YM8+fPp1atWvz6668MHjyYYcOG8dlnn+V4rpm/TPz5558MGjTI4i8YucWflTXnYzeWjg0aHR2tHnvsMdWjRw914cIFtWXLFhUcHKwqVqyo/Pz8VIUKFdRPP/1k8Zijha1Ix+LOze3bSoWG6sfonj7dfnEIIRzW6dOnlZ+fn7p06ZJx3cCBA1X37t0L9biDBg1SVapUUbGxsXlud/36deXt7a0WLVqU6zYuLi6qTZs2JuteeeUVde+992bb9tq1a+rkyZPGpVWrVmr69Okm69LS0mwWvzXnUxAFyT0WPwitWrUqP/30E19++SXt2rVj2LBhnDp1ilOnTpGRkUGdOnXMnsi71HF3h0mTICwMpk2DgQPBjEYTQojS49VXX+XatWtUrlzZuE6n0+U6Z7Qt5pweOnQoP/74I1u2bDE5bk7MmT85KCiIevXqmayrW7cuX3/9dbZtrZkvuiDxZ1VU80Fbw+pW3M888wy7du1i//79tG/fHp1OR5MmTSQ55+e556BBA0hMhOnT7R2NEMKBrF+/nq1bt7J3716TW9SLFy/m3LlzJCYmZtunILe4lVIMHTqUtWvX8scff5jVkMswf3JQUFCu27Rt2zZbd6cTJ05QpUqVfMu3hDXxZ2XO+diNNVX2n376SX3wwQfqt99+U0optWnTJlW/fn01cuRI4xzOjsYhbnEb/PCD/ja3u7tS58/bOxohhAO4e/euqlu3rnr33XezvXfu3DkFqI0bN9r0mIMHD1Y+Pj5q06ZNJvMzZ/47nt/8yVnnZFZKqZ07dypnZ2c1ZcoUdfLkSbV8+XLl4eGhvvjii2wxWDtftLnxZ42vRM8HHRkZqQICAlSfPn1U9erV1cSJE5VSSqWmpqpx48ap2rVrq59//tniQAqbQyVonU6p++/XJ+mBA+0djRDCAcyYMUMFBASolJSUbO/pdDrl4eGhZs6cadNjksv8zFFRUcZt8ps/Oac5mZVS6ocfflANGjRQbm5uqk6dOurTTz/NMYaCzBdtTvxZ4yvR80H7+fmxfv16mjdvTkJCAvfeey8nTpwwvn/kyBFeeukl/vzzT2sr9YXCbkN95ubvv6FtW3BygsOHoXZte0ckhBDCxop0PuiyZcsSHR0NQGxsbLZnzvXq1XO45OyQ7rtPP2BJRsZ/3a+EEEKIf1mcoKdNm0bfvn0JDg6mXbt2TJo0qTDiKh2mTgWtFr7+GnbutHc0QgghHIjFt7gBrl27xpkzZ6hVqxblypUrhLBsz+FucRuEh+tHFnvoIdiwATQae0ckhBDCRorsFrdhNBc/Pz9atmyZb3K+cOGCRcGUShMmgJsbbNwI69fbOxohhBAOwqIE3bJlS1566SV27dqV6zZJSUksXLiQBg0a5NgpXWQRGgpDhuhfv/466HT2jUcIIYRDsOgW97Vr15gyZQpLlizB3d2d5s2bExwcjLu7O9evX+fIkSMcPnyYZs2a8dZbb/HYY48VZuwWcdhb3ADXrkH16pCcDMuXw7PP2jsiIYQQNlCQ3GPVM+jbt2/z008/8ddff3Hu3Dlu376Nv78/TZs2pXPnznkOVG4vDp2gAaZMgXHj9DXqY8egTBl7RySEEKKAijxBF0cOn6Bv3YI6dSA2Vj9e97hx9o5ICCFyFBsbywsvvEB8fDzOzs689dZb9O7d295hOSRJ0GZw+AQNsHIlPPMMeHjAiRNQqZK9IxJCiGwuXbpEXFwcTZo04fLlyzRv3pwTJ05QtmxZe4fmcIp0oBJRiPr00Y8udusWjB1r72iEECJHQUFBNGnSBIDAwED8/f1JSEiwb1AlkCRoR6LRwKxZ+teffw47dtg1HCGEyM+ePXvIyMggJCTEpuVu2bKFxx9/nODgYDQaDd9++61F+7/77rtoNBpGjBhhsn7atGm0bNkSLy8vAgIC6NGjR7aZtxyFJGhH06IF9Ounfz1iBJSOJxBCiGIoISGBvn378umnn9q87Js3b9K4cWPmzZtn8b67du3ik08+oVGjRtne27x5M0OGDGH79u389ttvpKWl0alTJ27evGmLsG3KqgR97tw5fvjhB3bK8JSFY+pU8PSE7dthxQp7RyOEKGV69+5NhQoVTBLvjh07cHV1Zf2/AyqlpqbSo0cPXn/9de677z6bx9ClSxcmT55Mz549LdovJSWF5557joULF1K+fPls769bt45+/fpRv359GjduzNKlS4mJiWHPnj22Ct1mLE7QX375Jffccw/du3enTZs2tGjRgitXrhRGbKVXUBC88Yb+9Zgx4IDf7IQQJdfs2bPp1asXEydOBPRJ7/nnn2fw4MF06tQJpRT9+vXj4Ycf5oUXXsi1nKlTp+Lp6ZnnYhih0laGDBlC165d6dixo1nbJyUlAeDr62vTOGzB4lbcderUoU2bNowdO5bY2FjGjBlD48aNWbx4cWHFaBPFohV3ZnfuQL16EB0Nb7+tHxJUCCGKyLFjx6hbty5Xr15l1KhR7Ny5k927d+Pu7s5ff/3Fgw8+aHIL+fPPP6dhw4YmZSQkJOTbeKxq1ao4OzvnuY1Go2Ht2rX06NEjz+1WrlzJlClT2LVrF+7u7rRv354mTZowy9C2JwudTke3bt1ITEzkr7/+yrNsaxUo91g6gbSLi4vJ5NlHjx5VHh4eFk9EXdQKMmm23Xz1lVKglLu7UufO2TsaIUQpkpGRoTw8PNTLL7+sXF1d1b59++wWC6DWrl2b5zYxMTEqICBA7d+/37iuXbt2avjw4bnuM2jQIFWlShUVGxtro0izK0jusfgWd3p6Oh4eHsaf69Spg06n4/Lly5YWJfLzxBPQrp2+Nj1mjL2jEUKUIlqtloYNG/Lxxx8zefJkGjdubHEZRXmLe8+ePcTHx9OsWTOcnZ1xdnZm8+bNzJ49G2dnZzIyMky2Hzp0KD/++CMbN26kcuXKNonB1vK+r5CLzz77jLZt29KoUSM8PT1xdnbm1q1bto5NGLpdNWumH8Rk6FB9P2khhChk6t+nn82aNeO1117LdbuMjAycnJxyfG/QoEE89dRTeR4nODjY+iAz6dChAwcPHjRZFx4eTp06dRgzZowxRqUUr7zyCmvXrmXTpk1Uq1bNJscvDBYn6AceeIDJkydz48YNtFot1apV486dOyxevJiOHTvSokULvLy8CiPW0qlJExgwABYuhOHDYedO0ErvOCFE4Zo1axY7duygSZMmaLP8zenWrRuVK1dm165dvPTSSwwYMCDHMnx9fa1ufJWSksKpU6eMP0dHR7Nv3z58fX0JDQ1l7ty5rF27lg0bNgDg5eWVbR6IsmXL4ufnZ7J+yJAhrFixgu+++w4vLy/j3V8fHx/KONocCNbeVz9x4oT68ssv1ahRo1SHDh1U+fLllUajUU5OTqpOnTrWFltoiuUzaIO4OKW8vPTPo5cssXc0QogS7sCBA8rNzc34/DktLc3k/apVq6pZs2YVagwbN25UQLYlLCxMKaXU+PHjVZUqVfIsI6dn0DmVCaioqKhCOY+C5B6bjsUdHR3N7t272bt3L1OnTrVVsTZR7FpxZ/XBBzBqFFSoAMePQw79+4QQoqDu3LlDy5YtadasGbNmzcLX15eDBw8aa6E3btygfv36Nu8eVVI5zFjc1apVo3fv3g6XnEuEYcP03a6uXPmvj7QQQtjY66+/zs2bN5k7dy7ly5enSpUqzJo1i4sXLwJw+PDhQhmYRGTnEA8zzRkb9c6dOwwZMgQ/Pz88PT3p1asXcXFxdorYDlxd4eOP9a8/+UT/LFoIIWxo/fr1zJs3jy+++MLYlmjcuHF8++23DBkyBICDBw9m6+8sCodDJGhzxkZ99dVX+eGHH1izZg2bN2/m4sWLPPHEE3aM2g7atYO+ffXjcw8aBOnp9o5ICFGCdOrUibS0NJMa8oABA7h69Spr164FJEEXJYecD/rKlSsEBASwefNmHnzwQZKSkqhQoQIrVqzgySefBP4b5Wbbtm3ce++9+ZZZ7J9BG8THQ+3akJgIH32kv/UthBDCITnMM2hbyTo26p49e0hLSzMZW7VOnTqEhoaybdu2HMtITU0lOTnZZCkRAgJg2jT963Hj4NIl+8YjhBAiX1nzUWpqar77WDVQicGGDRvYsGED8fHx6HQ6k/eWLFliVZk6nY4RI0bQtm1bY6vBy5cv4+rqSrly5Uy2rVixYq4jmE2bNo0JJXX86hdfhKgo/XPoyEj48kt7RySEECIPWefLHj9+PO+8806e+1idoCdMmMDEiRNp0aIFQUFBaDQaa4syMWTIEA4dOlTggcvHjh1LZGSk8efk5GSbTyhuN1otLFignzt65Uro3x8eecTeUQkhhMhFbGysyS1uNze3fPexOkEvWLCApUuX5jnVmKUMY6Nu2bLFZGzUwMBA7t69S2JiokktOi4ujsDAwBzLcnNzM+sCFFtNm+qH/pw9G4YMgQMHwN3d3lEJIYTIgbe3d9E9g757967N+sIppRg6dChr167ljz/+yDY2avPmzXFxcTEO6QZw/PhxYmJiaNOmjU1iKJYmToTAQDh5Et57z97RCCGEsCGrE/SAAQNYsWKFTYIYMmQIX3zxBStWrDCOjXr58mVu374N6MdIjYiIIDIyko0bN7Jnzx7Cw8Np06aNWS24SywfH5g5U/966lTING6tEEKI4s3qblbDhw9n2bJlNGrUiEaNGuHi4mLy/owZM8wPIpfn11FRUfTr1w/QD1Ty2muv8eWXX5Kamkrnzp35+OOPc73FnVWJ6WaVlVLQqRP8/jt07gy//KKfBUsIIYTdFST3WJ2gH3roodwL1Wj4448/rCm20JTYBA1w4gQ0bAh378Ly5fDss/aOyGHodDqSk5Px9vbONiOPEEIUNrsk6OKmRCdogMmT4a23wM8PjhzR95cWpKenc/78eSpXroyzc4F6FQohhMVK3EAlwgpjxkCjRnDtmn7eaCGEEMVagaoUiYmJLF68mKNHjwJQr149IiIi8PHxsUlwwgIuLrBkCbRqpe8b/cwz0K2bvaMSQghhJatr0Lt376ZGjRrMnDmThIQEEhISmDlzJjVq1OCff/6xZYzCXM2bw8iR+teDB+vH6xZCCFEsWf0M+oEHHqBmzZosXLjQ+GwvPT2dAQMGcObMGbZs2WLTQAuqxD+DNrh9Gxo31veNHjAAFi60d0R2Jc+ghRD2ZJdn0Lt372bMmDEmf/ScnZ0ZPXo0u3fvtrZYUVBlysCiRfrXixZBpsFdhBBCFB9WJ2hvb29iYmKyrY+NjTVO9C3s5MEH4eWX9a8HDoRM82oLIYQoHqxO0H369CEiIoJVq1YRGxtLbGwsK1euZMCAATzzzDO2jFFY4913ISQEoqP13a+EEEIUK1Y/lPvggw/QaDT07duX9PR0AFxcXBg8eDDvvvuuzQIUVvLygk8/hS5dYNYseOopKM3DogohRDFT4IFKbt26xenTpwGoUaMGHh4eNgnM1kpNI7GswsJg2TKoWxf++afUzXgljcSEEPZk14FKPDw8aNiwIQ0bNnTY5FyqzZgBFSvC0aNyq1sIIYoRi6oUkZGRTJo0ibJlyxIZGZnntpZMliEKkZ+fvqtVt27w4Yfw+OP6RmRCCCEcmkUJeu/evaSlpRlf5ya32amEnTz+OEREwOLF+lve+/dDabrNL4QQxZDVz6BjYmKoXLlythmClFLExsYSGhpqkwBtpdQ+gzZITtYPYHL2rD5ZG/pKl3DyDFoIYU92eQZdrVo1rl69mm19QkIC1apVs7ZYUVi8vWHpUv1c0YsXww8/2DsiI51OR2JiIjqdzt6hCCGEw7A6QedW8U5JScG9lLUULjbatQND24EBA+DKFfvG8y9J0EIIkZ3F9/wMjcM0Gg1vv/22ScvtjIwMduzYQZMmTWwWoLCxyZNh3To4fBgGDYKvvtLXqoUQQjgUixO0oXGYUoqDBw/i6upqfM/V1ZXGjRsz0jCjknA87u7w+ef6aSm/+QaWL4fnn7d3VEIIIbKwOEFv3LgRgPDwcD766KPS2eCquGvaFN55B8aNg6FD9d2uHKxRnxBClHZWP4OOioqS5FycjRkDrVtDUhL07QsZGfaOSAghRCZWJ+hp06axZMmSbOuXLFnC9OnTCxSUKALOzvpb3WXLwubNMG2avSMSQgiRidUJ+pNPPqFOnTrZ1tevX58FCxYUKChRRGrVgo8/1r9+5x34+2+7hiOEEOI/Vifoy5cvExQUlG19hQoVuHTpUoGCEkXohRfguef0t7iffRYSE+0dkRBCCAqQoENCQti6dWu29Vu3biU4OLhAQYkipNHoa9HVq8O5c/Dii1CwCc6EEELYgNVjHw4cOJARI0aQlpbGww8/DMCGDRsYPXo0r732ms0CFEXA2xu+/BLatoU1a6BzZ/1woEIIIezG6gQ9atQorl27xssvv8zdu3dRSlGmTBnGjBnD66+/bssYRVFo1Uo/iMnrr8OwYfpknUMbAyGEEEXD6skyDFJSUjh69ChlypShVq1auLm52So2myr1k2WYQ6fT155//10/scb27fqBTQpZYU5oIZNlCCHsqSC5p8B/sWJiYrh27Rp3797lzJkzxvXdunUraNGiqGm1sGwZNGqkn5Jy1CiYM8feUQkhRKlkdYI+c+YMPXv25ODBg2g0GuPkGYa5oDNk4IviKSgIPvsMunaFuXP1o4z17m3vqIQQotSxuhX38OHDqVatGvHx8Xh4eHDo0CG2bNlCixYt2LRpkw1DFEXuscf0z6JB31jsxAn7xiOEEKWQ1Ql627ZtTJw4EX9/f7RaLU5OTtx///1MmzaNYcOG2TJGYQ+TJumnp7xxA558Em7dsndEQghRqlidoDMyMvDy8gLA39+fixcvAlClShWOHz9um+iE/Tg767teVawIBw/qJ9UQQghRZKxO0A0aNGD//v0AtG7dmvfee4+tW7cyceJEqlevbrMAhR0FBemTtFYLUVH6RQghRJGwOkGPGzcOnU4HwMSJE4mOjuaBBx7g559/Zvbs2TYLUNjZQw/BxIn61y+/DAcO2DceIYQoJQrcDzqzhIQEypcvb2zJ7UikH3QB6HTwv//BL7/oJ9jYvVs/+piNSD9oIURJVZDcY1UNOi0tjQ4dOnDy5EmT9b6+vg6ZnEUBabX6qSlDQuDkSf380f/ePRFCCFE4rErQLi4uHJBbnaWLnx98/TW4ucF338GUKfaOSAghSjSrn0E///zzLF682JaxCEfXsiXMn69/PX48/PijfeMRQogSzOoEnZ6ezvz582nRogUvvfQSkZGRJosltmzZwuOPP05wcDAajYZvv/3W5P1+/fqh0WhMlkcffdTa0EVBhIfrG4sppZ9HWgYxEUKIQmF1q5lDhw7RrFkzAE5k+SNt6XPomzdv0rhxY/r3788TTzyR4zaPPvooUZm6+TjqpBylwsyZ+tbcf/0FPXroJ9WwQ8M7nU5HcnIy3t7eaLVWf9cUQgiHZHWC3rhxo82C6NKlC126dMlzGzc3NwIDA212TFEArq76eaObN4ejRyEsTP98uoiTpE6nIzExEU9PT0nQQogSx+K/amfOnMGGPbPMtmnTJgICAqhduzaDBw/m2rVreW6fmppKcnKyySJsKDAQvvlGn6y//VYajQkhRB6y5qPU1NR897E4QdeqVYsrV64Yf+7Tpw9xcXGWFmORRx99lGXLlrFhwwamT5/O5s2b6dKlS54zZk2bNg0fHx/jEhISUqgxlkqtW8PHH+tfv/22vhYthBAim5CQEJOcNG3atHz3sXigEq1Wy+XLlwkICADAy8uL/fv322x4T41Gw9q1a+nRo0eu25w5c4YaNWrw+++/06FDhxy3SU1NNfmGkpycTEhIiAxUUhiGD4fZs8HDA/78E/5tm2AuawcTMWc/GahECGFPhoFKYmNjTXKPm5tbvm2piuWDu+rVq+Pv78+pU6dy3cbNzQ1vb2+TRRSSDz+Ezp31M1516waXLtk7IpsxPOfWycAsQogCyJqPzGnobHGCNnRzyrquKJ0/f55r164RFBRUpMcVuXB2hlWroE4duHABuneH27ftHRUAFy9e5KOPPjLOtmYpSdBCCHux+J6fUop+/foZs/+dO3cYNGgQZcuWNdnum2++MbvMlJQUk9pwdHQ0+/btw9fXF19fXyZMmECvXr0IDAzk9OnTjB49mpo1a9K5c2dLwxeFxccHfvhB/1x61y7o3x9WrAA7D/26ePFijh07xuLFi5kwYYJdYxFCCEtYnKDDwsJMfn7++ecLHMTu3bt56KGHjD8bBjoJCwtj/vz5HDhwgM8++4zExESCg4Pp1KkTkyZNkr7QjqZmTX1DsUcegZUroW5dfeMxO4qIiCA5OZmIiAi7xiGEEJay6WxWjkxmsypCixbBwIH61198oR9xLA+O3EhMGpkJIQqiyGezEiJPAwbAyJH61+HhsGmTXcMRQojiSBK0KBzTp8OTT0JaGvTsqR9xTAghhNkkQYvCodXCsmXQpg0kJsJjj0EhD2gjhBAliSRoUXjKlIHvv9c3Hjt7Fv73P7h5095RCSFEsSAJWhQuf3/45Rfw84Pdu+HppyE93WSTgvZVFkKIkqhAzVI3bNjAhg0biI+PzzaQw5IlSwoUmChBatbU16Q7dIAff4SXXtK39P63j7T0VRZCiOysrkFPmDCBTp06sWHDBq5evcr169dNFiFM3Hefvm+0VgtLlsC4cca3IiIiqFOnjvRVFkKITKyuQS9YsIClS5fywgsv2DIeUZJ17w6ffKLvIz11KlSsCMOGERwczPDhwwkODrZ3hEII4TCsrkHfvXuX++67z5axiNJgwACYPFn/esQIfa1aCCFENlYn6AEDBrBixQpbxiJKizfegKFDQSno2xfNb7/ZOyIhhHA4Vt/ivnPnDp9++im///47jRo1wsXFxeT9GTNmFDg4UUJpNDBrFsTHw+rVaJ98ErfPPoPKle0dmRBCOAyrE/SBAwdo0qQJAIcOHTJ5r6innxTFkJOTfiCTpCQ0v/5Kxf790VWpAi1b2jsyIYRwCFYn6I0bN9oyDlEaubnBN9+gOnVCu3Urmq5dYcsW/bzSQghRyslAJcK+PDzI+O47Uhs0QHPlin6qynPn7B2VEELYXYEGKklMTGTx4sUc/XcihHr16hEREYGPj49NghOlhI8PcUuXEvL882iOHdMPaLJpkzyTFkKUalbXoHfv3k2NGjWYOXMmCQkJJCQkMHPmTGrUqME///xjyxhFKaDz8yNj3TqoVg1On4aHH4YLF+wdlgxDKoSwG6sT9Kuvvkq3bt04e/Ys33zzDd988w3R0dH873//Y8SIETYMUZQalSrBxo1QtSqcPAkPPQR2ToyZhyEtbXQ6HYmJidmG8RVCFI0C1aDHjBmDs/N/d8mdnZ0ZPXo0u3fvtklwohSqUkWfpKtU+S9JX7pkt3BK8zCkkqCFsC+rE7S3tzcxMTHZ1sfGxuLl5VWgoEQpV7WqPkmHhsKJE/okffmyXUJx1GFIJXkKUfJZnaD79OlDREQEq1atIjY2ltjYWFauXMmAAQN45plnbBmjKI2qVdMn6ZAQOH5cn6Tj4kw2Kc3PhyVBC1HyWd2K+4MPPkCj0dC3b1/S/53f18XFhcGDB/Puu+/aLEBRilWvrk/S7dvDsWP6hmN//KGfZAOZplIIUbJZXYN2dXXlo48+4vr16+zbt499+/YZW3K7ubnZMkZRmtWooU/SlSrBkSP6JB0fD5g+H86pRim1TCFEcVbggUo8PDxo2LAhDRs2xMPDwxYxCWGqZk19kg4O/i9Jx8WZPB/Omox1Oh0JCQlcv34dpZSdTyA7+fIghMiPRbe4IyMjmTRpEmXLliUyMjLPbWWyDGFTtWrpBy9p1w4OH4YHH4RffgFtzt8xdTodSUlJDpmc4b8E7enpiTaXcxBClG4WJei9e/eSlpZmfJ0bmSxDFIpatWDzZujYEU6cwKl9e5yXLpURx4QQJZJFCTrzBBmfffYZlStXzvbtXylFbGysbaITIqtateCvv6BjRzQnThD01FOo9euhUSN7RyaEEDZl9b21atWqcfXq1WzrExISqFatWoGCEiJPISHw55+oxo1xunYNpw4dSPn9dzIyMuS5rhCixLA6Qef2bC8lJQV3d3erAxLCLAEBZPz+O3eaN0eTmIh3r164/vUXSUlJkqCFECWCxf2gDY3DNBoNb7/9tknL7YyMDHbs2EGTJk1sFqAQuSpXjrhly6g8bBhOGzZQaeBA4mbPhgED7B2ZkU6nIzk5GW9vb2kMJoSwiMUJ2tA4TCnFwYMHcXV1Nb7n6upK48aNGTlypO0iFCIPysODtK+/5vZTT+G5fj2BQ4eiK1MG9dxz9g4NkNbaQgjrWZygDQ3FwsPD+eijj/D29rZ5UEJYxM2NuNmz0Y0Zg/d33+HUrx8Z8fHQs6e9IxNCCKtZPdRnVFQUAEeOHCEmJoa7d++avN+tW7eCRSaEJZydufzuu+jKl6fc0qU4jR6N3+HDxI8ZA05O9o5OCCEsZnWCjo6OpkePHhw8eBCNRmNsNGboA52RkWGbCIUwl1bLtTffxKtePZxGj6bcZ5/hdPkyVz780N6RCSGExax+KDZs2DCqVatGfHw8Hh4eHD58mC1bttCiRQs2bdpkwxCFsIBGg3r1VTK++ALl4oLXr78SFBYG167ZO7JipzTPFiaEI7A6QW/bto2JEyfi7++PVqtFq9Vy//33M23aNIYNG2bLGIWwmOrTh4tRUWR4eVFmzx6c2rWDs2fz3EfGxzaVebYwIUTRszpBZ2Rk4OXlBYC/v7/xW3aVKlU4fvy4baITogDu3HsvsV9+SVpQEJrjx6FNG/jnn1y3lwRtKvNsYUKIomd1gm7QoAH79+8HoHXr1rz33nts3bqViRMnUr16dZsFKERB3K1Viwtr1qAaNoTLl/WTbPzwg73DypcjfFnIPFuYEKLoWZ2gx40bZ/zjMXHiRKKjo3nggQf4+eefmT17ts0CFKKgMgIDydi4ETp0gJs3oXt3eP99cNCZriD/BC3Ph4Uo+axO0J07d+aJJ54AoGbNmhw7doyrV68SHx/Pww8/bFFZW7Zs4fHHHyc4OBiNRsO3335r8r5SirfffpugoCDKlClDx44dOXnypLWhi9LIx0c/PeVLL+kT8+jR0L8/pKbmuZtWq8Xb25vk5GSb1mYLmmDl+bAQJZ9Nhzby9fW1aqrJmzdv0rhxY+bNm5fj+++99x6zZ89mwYIF7Nixg7Jly9K5c2fu3LlT0JBFaeLiAvPnw+zZ+nmkly7VT1155UquuxRWgi5ogrXX82FHuPUuRGlhdYJ++OGHmTBhQrb1169ft7gG3aVLFyZPnkzPHEZ+Ukoxa9Ysxo0bR/fu3WnUqBHLli3j4sWL2WraQuRLo4FXXoGff9bXqv/6C1q1gkOHijSMgiZYez0flgQtRNGxOkFv2rSJuXPn0qNHD27evGlcf/fuXTZv3myT4EA/IMrly5fp2LGjcZ2Pjw+tW7dm27Ztue6XmppKcnKyySKEUefOsG0b1Kih737Vpg2an34qssNLAywhSpes+Sg1n8drUMBb3L///juXL1/m3nvv5Ww+fUytdfnyZQAqVqxosr5ixYrG93Iybdo0fHx8jEtISEihxCeKsbp1YccOaN8eUlLQ9uyJz+zZUIJqh1LjFcIxhISEmOSkadOm5btPgRJ0UFAQmzdvpmHDhrRs2dKhRhAbO3YsSUlJxiU2NtbeIQlH5OcHv/4KgwahUYryM2ei7dULkpJMNnOERGdNDI4QtxACYmNjTXLS2LFj893H6gRtaAzm5ubGihUrGD58OI8++igff/yxtUXmKDAwEIC4uDiT9XFxccb3cuLm5oa3t7fJIkSOXF1h/nwyPv0U5eqK9scfoWVLOHzYuIlSiqSkpGKXoIUQjiFrPnJzc8t3H6sTtMrSh3TcuHEsX76cD208MUG1atUIDAxkw4YNxnXJycns2LGDNm3a2PRYwj60Wi3lypWz+XzJcXFxzJ07l9u3b5tVtgoP59KaNaiQEDh5Elq3hjVrChSDTqcjJSWlVCdV+WIhhHWs/osYHR2Nv7+/ybpevXqxfft2lixZYlFZKSkp7Nu3j3379hnL3rdvHzExMWg0GkaMGMHkyZP5/vvvOXjwIH379iU4OJgePXpYG75wIIWRoC9evEhkZCQHDhxgzZo1Zpd9t1EjMnbsgIcf1g9q8tRTaMeMgbQ0q+KQBC0JWghrWT3dZJUqVXJc36BBAxo0aGBRWbt37+ahhx4y/hwZGQlAWFgYS5cuZfTo0dy8eZMXX3yRxMRE7r//ftatW4e7u7u14YsSbvHixbi7u3Pnzh3LuzJVqKB/Lv3GG/D++2hnzCB482biZs2CXP7flxbyhUOIomNRgo6MjGTSpEmULVvWmERzM2PGDLPLbd++fbZb5plpNBomTpzIxIkTzS5TlG4REREkJibSp08f67oyOTvDe+9Bq1ao/v0ps2cPId26weefQ9eutg84E0eucUqCFqLoWJSg9+7dS9q/t/r27t2b63bWjCYmhKXyujUeHBxsm2lPn3ySjAYNyOjVC7cjR+B//4OxY2HiRH0SLwQ6nY6kLK3IC5tOpyM5ORlvb2+btwUQQljHor8wGzduBCAtLQ2tVsuCBQuoVatWoQQmRH4MCTo9PZ24uDhWrlzJ008/nevjF6vVrMmFNWvwmzoVn+XLYdo0/QhkX34JlSrZ9lh2Yqi1e3p65pmgL168yOLFixk7dqzMWidEIbPqq7KLiwsHDhywdSxCWG3VqlUcP36cVatWFUr5ys2NqxMmkLFiBXh5wZ9/QpMm8N13hXK8omTJLfWoqChOnTpFVFRUEUQmROlm9b2s559/XmbSEQ7h4sWLpKSkUKlSJfr06VOox1K9e8M//0DTpnD1KvTooZ8hK9Nwt/ZQkOfWhlvqebUDMQgPD6dmzZqEh4dbE6YQwgJWJ+j09HTmz59PixYteOmll4iMjDRZhCgqixcv5sKFC4C+Jl3ocyTXrKkfx3vkSP3kG59+Cs2awe7d2TY13BI2NyZDor106RKzZ8+2eD9bNt7K6Rl/cHAwERERFjW8k7mrhbCO1Qn60KFDNGvWDC8vL06cOMHevXuNi6E/sxBFISIigtq1awNw/Pjxormz4+YG778Pv/+ufw594gS0aaN/Pp2RYdzM0lvChtrs6tWrbXIuBUmOtuqfLnNXC2Edq5uhGhqMCWFvhhbbFy9e5KuvviraOZIffhgOHIBBg/Sjjr3xBqxbB8uWQZUqhIeHc/nyZYtvCffu3dvsc8mr9pw5OeY0PWxRiIiIIDk5ucjnrhaiuJP+FKLEqFixIsOGDSv6KRx9fWHVKli6FDw9YcsWaNwYvvwy11vC+dVOLTmXvJ4hF3TeaVuQqTWFsE6BO3IeOXKEmJgY7t69a7K+W7duBS1aCLsy9A1OTExk9uzZ9OnTJ/cuXBoNhIXB/ffD88/D9u3w7LO49OyJ84AB2TY3JOjCJslRiOLL6gR95swZevbsycGDB9FoNMZv74ZBSjIyPYcTojgwPK999dVXCQ0NNd46Xrx4sbELV6tWrfIupEYNfResKVNg8mSc1q6l8R9/oObM0SfuTIP4FGRwkMw1cBnVS4iSyepb3MOHD6datWrEx8fj4eHB4cOH2bJlCy1atHCoeaGFMFdujZnCw8OpXbu2+V24nJ1h/HjYuRNdo0a4JCXh2rcv9OoFly8bNytIy+vCmgFMCOE4rP50b9u2jYkTJ+Lv749Wq0Wr1XL//fczbdo02wyxKEQRy+15raERWsWKFS0rsGlT7v75J7EDBqCcnWHtWqhfH5YvBzP6HBeluLg4i7p1CSEKn9UJOiMjAy8vLwD8/f2NH+wqVapw/Phx20QnRBHK63mt1QnM1ZULERHc3bpVP7hJQoL+VnePHnDpkm0CtwHDSGyO2BXKkScPEaIwWZ2gGzRowP79+wFo3bo17733Hlu3bmXixIkyRq8ocQqawFSjRrBjB0yaBC4u8P33ODVsiNfy5eAAiadPnz7Url0739beWq023/G6bU0StCitrP6UjRs3zviBmThxItHR0TzwwAP8/PPPzJ4922YBCuEIzE1geXJxgXHjYM8eaNECTVISfuPG4fTgg/q+1HZy8eJFVq1aZdbUnPZI0EKUVlZ/yjp37swTTzwBQM2aNTl27BhXr14lPj6ehx9+2GYBCuEIbNrHumFD2L6djJkz0Xl6otmxQz9U6OjRxjG9c5p32ZZDZmaulWZupV4YpEGbENax6SfG19dX5oIWIpNcx+J2ckINHcqF335D98QT+uFB338fp8aNSV67lrlz53L27FmTBJ3XkJmWPiPPnKANQ6UW1kQjkqCFsI5F/aAtmQRjxowZFgcjREmTeSzuSZMmZXs/IzAQ3apVaNetgyFD0Jw7x/ejRhETEMBPV67QuXNn47Z5DZmZ+Rm5pUN6BgYG0r9/fzw8PCw/QSFEobEoQe/du9es7aQWLYSe2WNx/+9/8NBD6MaPZ8DMmRAfT/9vv8Vp6lQYOxbKlMmzlXmfPn1YuXIlvXv3RqfTWVRbNdxOL1OmTJ7bJCcnF5skXpBBYIRwFBYlaJkgQwjL5Dc9o+FWs6+vL9qyZdG9+y66hx9m5DvvUHbXLn2r788/hxkz9Ek8FxUrVmTo0KE4OTmZJGhb3V42xJmQkMDs2bMZMGAA9erVs3nyyymxZh3hzZJ4pUGbKM4K9D83MTGRDz/8kAEDBjBgwABmzpxJUlKSrWITosRTSpGUlGTyrPlu7drELF3K8QkTUJUqwdmz8MQTaB97DJdTp7KVcfHiRWbPnk1cXFy292z9/DcqKopjx46xcOHCQun2lFOXKpmuUpRWVn9qd+/eTY0aNZg5cyYJCQkkJCQwY8YMatSowT///GPLGIWwikajwcfHp3jWoDQarnboQOr+/fDmm+Dqivb33wnu0gVtZCRcu2bc1NAKe82aNYUelmHY0969exd632RDsg4PD7f7jFxC2IPVf7leffVVunXrxtmzZ/nmm2/45ptviI6O5n//+x8jRoywYYhCWE6r1eLr66u/dVwME7RGo9HfnvXygsmT4cgRdP/7H5r0dLRz5ugn5Xj/fbhzx9gKu3fv3jaPI2uN1tCgzN/fP1vNv7COHRgYKDNyiVKpQDXoMWPG4Oz832NsZ2dnRo8eze7du20SnBDWKq5dewy3q69cuWL6/LRGDXRr13J52TJUw4aQlKTvN12nDpU2b2bY0KGWjxVuhvT0dM6fP096ejrwX4OynOaeFkLYltV/vby9vYmJicm2PjY21jhGtxDCMvndrr7zwANk7NoFUVFQqRKcO4dT375U6tWLMjt3Wn3c3B4H5DRgij3JsJ+iNLE6Qffp04eIiAhWrVpFbGwssbGxrFy5kgEDBvDMM8/YMkYhSg2zblc7OUG/fnDiBEyZgvL0xP3gQUJeeIHAF1+Eo0ctPq5Wq801Qd++fTvfhFjYidNQfnp6uiRoUWpYnaA/+OADnnjiCfr27UvVqlWpWrUq/fr148knn2T69Om2jFEIu7Dl0JrmsmhqSw8PeOMNMo4dI+nZZ1FOTpT94w+cmjSBsDA4c8bqODInREdI0Dm1dheipLM6Qbu6uvLRRx9x/fp19u3bx759+0hISGDmzJm4ubnZMkYh7MLQvScqKqpQys+re5RFKlbk6sSJnP3hB1IeeQSNTgfLlkHt2jBoEJw/n+OxDV8+chqOtLjcSi4ucQphDYsTtE6nY/r06bRt25aWLVsyceJEatasScOGDYvNKENCmCMiIoI6derkPwqYlWw9SUVajRrEzZ9P+t9/Q+fOkJ4On3wCNWvCq69CfLzJsQ19i5cuXcrZs2dZunSpWce5evUq8+bNK/gXCzMZbr/nNEKhJGhRklmcoKdMmcIbb7yBp6cnlSpV4qOPPmLIkCGFEZsQdmUYWrNy5cq5JoiCyGuSCksnvzDRsiWsWwebN8MDD0BqKsyaBdWrwxtvwLVrxi8fERERDBgwgAYNGjBgwACziv/+++85cOAAkZGRhXr731DLv3z5cqG0yJfkLhydxf/jly1bxscff8yvv/7Kt99+yw8//MDy5cvlP7mwm7xqWLYqvzASRF7Pm9esWWOc/MJqDz6oT9K//qpP2jdvwrRpUKUKlWfPJvLZZwkODs5zjG+DS5cu8dFHH3Hp0iW6devGnTt3cHd3Nys+axNhTiOIFdaUm0I4Iov/4sTExPDYY48Zf+7YsSMajaZIG9IIkZlhUJLy5cs73EQtWq3WqvGge/fuTe3atc0ePUuj0XDr1i3mzJlj+lnUaKBTJ9ixA779Fpo2hZs30X74IZUeeADt8OGQQ3fJrJYtW8axY8dYunQp/v7+vP/++zRu3Nis+KxNhOHh4VSvXp3w8HBjYp49e7ZZw37ao4GfELZmcYJOT0/H3d3dZJ2LiwtpaWk2C0oISxX2yGHW1qKtTdAVK1Zk2LBhZo+epdVq+frrr3OvdWs00L077NkDP/6Iat0abWoq2o8/xqlOHfxefx1On861/L59+1KnTh369esHQEBAgEXxWSMwMJAXXniBwMBAY20aMGvYTxm/W5QEFs1mBfruDv369TNpqX3nzh0GDRpE2bJljeu++eYb20QohJkMSbS4lW0r4eHhpKSk5J28NBro2pWMTp24umYNAQsXot20Ca9Vq1Br1sAzz+ifU99zD6C/tb1ixQrefPNNhg8fTlpaGkuWLKFv375Urly5iM7sv7mwX3nlFbRabb5fDPKaO1uI4sLiqkZYWBgBAQH4+PgYl+eff57g4GCTdUKIomXOs2QjjYY7992H7rffSN+8mVsPPaTvnrV8OdSvj7Z7d9y3b2fevHls2bKFefPmAbB06VJOnz5daBNz5HZr2qJzs2J7A3kuLRyJxTXowuoTKoTQ02g0eHt7c+vWraI54H33Eb9kCZXj43F+7z345hu0P/9M4M8/4+rvj7NOB/8mrH79+nH58uUcRzozmdvayscMmW9Njx07tkCnZQ1z55HOad5qIWxN/mcJ8S97TLCRU3cqu0300awZfPUVHD+ObtAgdO7ujLx6lbCEBF5fuxavTz+lbHo6/fr1IyAgINvuthjtK3P3L0cmNW1RFCRBC/EveyTGVatWFbg7lc3jrlUL3Zw5nN+6lfJvvsnr5cpR+fJl/KZNo2LLllSbOxeXLC2/cxsVzdJEZrg1HRgYKEN7ilKvWCTod955B41GY7LUqVPH3mEJka/8kmefPn0s6k5lzTGscfHiRWZ+/jk7H3mEAY8+yulJk7hbowZON29SefVqanTpQuCAAWjWrUOXns68efM4duxYtlHRdDodSUlJeU5PmVMSN9xCLlu2rMl52fJcpSuWcHTFIkED1K9fn0uXLhmXv/76y94hCZGv/BKKJd2prElOBR0kZNy4cZyJiWHRzZvE/vILMfPnk9C6NRqlKLtpE06PP46mTh1euHWLOlWqWDUqWtYYDUk9Li6OJUuWcPny5VyvQUGSrHTFEo6u2CRoZ2dnAgMDjYu/v7+9QxKiSBVWgs6pXMOz4MmTJ1OtWjX69u0LWi0327Xj6Icfcnv/fu68/DLKxwfN6dM0mDuXmZ9/ToOPPoJ9+0zKz3obP7+YDM+yV69ebdxPq9Xi7e1NcnKyyX4FSbI5Pe+WWrVwJMUmQZ88eZLg4GCqV6/Oc889R0w+ox+lpqaSnJxssgghssspQRueBTdu3Jhnn32WoKAgk30099yD+7x5aC5cIGP+fFLr1EF75w7eq1fj3LIl6a1bo1uyBG7ezHYb39xafebR1HJL0AVpVJZTV6xFixZx6NAhFi1aZHF5QuQlaz5KTU3Nd59ikaBbt27N0qVLWbduHfPnzyc6OpoHHniAGzdu5LrPtGnTTPplh4SEFGHEQjg2a5/lajQa0y5IZcuiBgzg/A8/ELN8OTe6dkU5O+O8cyfaiAicQkKoP3cuox56iOAsST4/5tz+t7a/c2769etnnNs+L1LTFpYKCQkxyUnTpk3Ld59ikaC7dOlC7969adSoEZ07d+bnn38mMTGR1atX57rP2LFjSUpKMi6xsbFFGLEQjs3mQ5dqNNxp0YL4jz4i9cQJro8ejapRA82NG/h8+SWVe/bEqUULmDsXrl+34ZmYKmj3p+DgYCIiIrIl/KzlZr21Lt2uRH5iY2NNcpI5/fyLRYLOqly5ctxzzz2cOnUq123c3Nzw9vY2WYRwdHbrA50Pc2YMMzYGU4qkwYPJOHKEjN9+40a3buhcXdEcOACvvIJTaCj+I0ag+e03yMiwaZy5tQgvaPLMWkbWW+uSoEV+suajzMNl58ax/gqYKSUlhdOnT2d7LiZEcedoCVqr1VKmTBmcnZ0pV64cTk5O+Pj45BifYYpM42iDWi2qfXviZ8zg3N9/kzFzJjRsiObOHTy/+w6nxx6DkBC0o0fjcvSoSVlxcXHMnTs3W79qS+I2PLNOT083K3lakmRtfWtdiJw4xl+BfIwcOZLNmzdz9uxZ/v77b3r27ImTkxPPPPOMvUMTJZg9k6U1xy6MeA0J2lBmXscwNOoKDw83rjMMYHIpNRU1dCjs30/61q0kP/88ytcXLl1CO3MmlR57DKdmzdDMmIFTXJwx2WfuV23Jc9/cGpXlRafTcfLkSRYvXizPloVDKBYJ+vz58zzzzDPUrl2bp556Cj8/P7Zv306FChXsHZoowSRBWzZdZk6NuhYvXmySaHVKcczbm7f9/IjduRPWruX8o4/yjpMTlw4exGnMGFzuvx/N119T7fZtnuna1aSswu63vGbNGk6fPi1zDgiHUCwS9MqVK7l48SKpqamcP3+elStXUqNGDXuHJUSBmPNc194yJ+jcarC5DfMJ+me1tWvXNg5gotPpWLhwoT7RLlsGPXqwsFUrdrRrxyePPYZq04ZFSnH+4kWC/vyT1o8/ju7xx9F98QURffoU+jjdTz31FHXr1s12DGm1LeyhWCRoIUoiR3venJ/carCGWnJOU1AGBwczbNgwKlasaFyXtV90REQEdRo1ImL+fDK2bKHLqlXUbNWKsCpV0Ny9i+u6dWhfeIGQ5s2ZfOwYlbZuhZs3C3w+ma+/4flzxYoVGTlyZLa5rs2tvUtjMWFLxeMvgxDCKras+eU2KIihlpzTFJSg7zuduWFZ1lvhWRtc+TVvTsTnn5Px+++k7tpF4pAhqJo10dy5Q9l163B69lkICEA99RS3Pv8clWk8hJzON7drkDVBnzhxgjlz5uR6JyDzuec25Oj58+clQQubkQQtRAlmy+e2ubVczlpLzpoQrb5ToNGgGjQgceRIMo4cIX3nThIHD0ZVqwa3bqFZswbv/v2p0rIl2p49YckSFs+ebTxfQxyzM63Li2FI0qwTfuR07lnPyXCdlyxZQkpKiiRoYROSoIUowawdCjNzAjInsRpqyRqNhqioKKu/FOT6XF6jgaZNSRw9mozjx2HXLnSRkaSFhKBNTUX7448QEcGLM2bQ4sABBmRksPiDDzh27BiAWdegV69eVK5cmV69epmsN+e2teE6h4WFWZyg5ba4yI2zvQMQQhQea/vrGhJ0enq6WYnDkFhv375NeHg4KSkpVjXmMhw3KSkp9400GmjRAl2TJsQMHozryZME79qF0/ffU2nPHiZdvAhTpvASoPz8GFCzJs4dOhAQGJjnsStUqEC3bt2y9Q4xzK6VF8N1TktLY/HixYwdO5bq1aubdc6GBO3p6Qnox2z29vYuNm0TROGR/wFC2FFxayhmjsIaxCOna6XRainTsiWaceNg927ST5/m2jvvoHvoIYKcnJh47Rqhc+cS3LMnTpUrQ9++sHIlJCSYlH3x4kXmzZtHQkJCgZLj0qVLOXXqlNXdtHQ6HQkJCSQkJEiNWkiCFsKeSlKCNudcsnbJyu32btaGZbmVn21daCg3wsLQrV9PxsWLXPnwQzJ69kTn5YXmyhX4/HN45hmoUAEeeACmTYP9+1m8aBEnTpzgt99+w8vLK9sAJznFk5lOpyMlJYXOnTuTnJzMY489ZvH1MzBMtykJWhT/vwpCCIdgToI2dMn67rvvjK2nk5KSUErlWNbly5etb4Xu68vNJ54gbcUKYvbsIX3DBhg9Gho0AJ0O/voL3ngDmjThxU8/pXlsLH29vVFxcdm+NOR3boYE/csvv+Dt7c1PP/1k9nPlzI3qjCOvXbokz6WFJGghROHJmtgiIiKoW7cuQ4YMyTXZZd5n8eLFHD16lHnz5hWs4ZWLCzz4IEyfDgcPwrlzsGABdOsGHh5Uiouj97FjTP/kE07fcw/BXbqgHTUK1q3j4qlTZn9J6Nu3LzVr1iQsLMzsBJu5pb3hC8zq1aulFi2kkZgQovAYkq1BTs+nNRqNyXPfzPtERESQmJhI9+7d0el0xlp3cnIyHh4euR7XsE2uz5NDQ+Gll/TLnTtkbNzIsP79OXPlCq9lZLDx2DE4dgxmzWKJVku0jw+L9+1jwrRp0KIFOOf8pzMoKIiIiAiCgoK4evWqWdcoIiKC5ORkY6O6xMREnnzyScvOR5RI8psWQhRIQZ+j57V/TiORGWrH58+fz3UAkszlGmbiyjU+d3fUI4/w2scfE9ioEcOjohjdpQvn+/SB0FBe1Oloev06L23aBG3aoPz8UF26cPudd9Bt3QppaSbFWTp0a+YvLcHBwfTp04c1a9aYDJgiXbFKJ0nQQohcmTMSmb0auuXW3zprEjM3vvr16zNnzhx2nTzJQaVYWLs2nD1LwNGjDJ08mYo9epDh44MmORnNunWUmTAB7f33Q7lyXO3Qga+HDePqDz/AnTtcunTJ6mfneQ2YIkoXSdBCiFwVxQxS1goPD7fZ5BmG6Sk1Gg39+vX7r1yNBl316lzq3p3U5cuJ3bOH9B07yPjgA2527ozy84Nbt1j899+c2bWL5X370vKRR/iiQwdOf/89i994A/6t+WauAWf+OesXCMNY5YYJRiwlte2SQxK0ECJXmUciM7cmWlQ1alv0tzYkM8AYc1BQkEm5hhbaOp0OnJygWTPU8OFcWbCAjIsX4dAh+k+eTJ177mFAhQpo09J4KTaW5mfO8NLy5eDri6ZhQ5xeegm1aBEcOYIuPT3HBK3Varnnnnt46qmnWLVqlbEGfv78ed5//33Onz9vEndOSTi39yRxFz+SoIUowQqaLDMnwcJI0ObcQo+Li2P27NkW3S42d5KQrDVZc+e+NtJqoX59gsaMYeCvv+J37hznfv8d7eTJRPboQVCNGqAUmiNH8Fq9GqeXXoL69XGqWJGKYWFoJk2C9evh3y8Jhmu3Zs0ajh8/brxzsXjxYo4cOWL82ZpkKwm6+JEELUQJ5ugDoZhzC93wTNaS2+xRUVEWd89ydnamcuXKOGdpoZ1T4s51buxLl5jx/fecfPBB4j/4gIxjxyAujthFixjbogUXWrcGDw8uJSby3pYtxE2cCJ07g68v1KsHYWFo5szhhXr1qFOjBhEREeh0Op588klq1KhBeHi42ddAFH+O+akVQpQK5kzmkXX+6Lxqx4b3HnvsMWrXrm3snmWO3Fp9GxJ05kFTcvtiERUVlX1u7IAAFsXEsM/fn087d4akJD4ZOJDddevyaYMGUL06KAVHj8KyZThFRtJi0CBmLlhASNeuEB4On32GS1wc3Llj9jXI9uXBhlOPiqIh/aCFEEUma40+v+fIhmeymbcxJMeoqKhsid3w3s8//8ywYcMKFKOBIbEppYw1+ax9lw37GNb37t3bZGhQk+2dnYkYN47kmTPp/+qr+j7ZcXGwe7dxUbt3o7l8GQ4dQnPoEN8C8cCyn3/mnYYNWaLRcCo1lcWJiUz46CN0np7GftKZvzxMmDAh27XJul44LknQQogikzX5mbt9SkqKcZ0h2eV0uzfze1mHD7WWYTSzatWqGWvyuc0P7enpSZ8+fVi1ahWvvfYa/v7+wH9fRAIDA0lMTCQwMND0i0nFitC1q34BMtLSuLRnD8GXLqHZs4fnfv+djP37efHOHdi/nxeBhcDAo0dh6VI01arhUqsWtG7NgEqVSAkMJCLL9cn6paKoyCAr1pMELYQoVjInu8uXL+f63oULF/Ity1A7fvXVVwkNDc1xG8NoZk8++aRxMJG8ZH5mnrWmamio5e7unndgGg0ZgYGoFi240Lw5s5OT6frOO6R6eJARH0/gnj2M2rGDMidPwvnzaKKjKRsdDevXEwJ8CKivvoJGjaBxY2jcmOC6dXk1PNzms4zlJ/N0mpKgLSMJWghRaplz29cwmllGRoZZZRpq0HnVVE26buXz/uLFizlx4gQANYcPR913H6pHD+LPn9c3aktKIva335g3fz5DK1ak0qlTcPgwmpQU+Ptv/YL+j30ooIKD9Q3SDEv9+vp/fX3NOj9LmfMlSORMErQQoljKq4W6Vqvl9u3bLFq0iMjIyDxrx7a+7VuxYkWGDRtmUlPNOjvX0KFDzS4vIiKCpKQk+vXrl/Mwon5+LDp6lP0eHiysW5e3VqzgYkwMwTdu4Hz4MBw4APv3ow4fRhMbi+biRbh4EX7/PWvg/yVrw1K7tn69hcOXZibPvq0nCVoIkStH7qaV1/NsrVZr0pc4r9qxOYOdxMXFsXLlSp5++mmqVKmS63XRarX4+PiQnJycY6wzZ87k2LFjLF261Dj4S27xG24JBwcHM2LECHQ6HXPmzOHVV1+lcuXK2WYJy/xFQ/dv/2waN4ZnnwVA6XQknT+P94ULaI8dgyNH4PBh/b/nzukbqsXFwR9/mAbj5QX33JN9qVULfHzyvG45xSbMJwlaCJErSxt1FWYM5nxJyLytLRND5vGxW7VqZfV1yauBW2aZbwsbvjxkHnt8woQJ+c4SlpVWq6VcaKi+1XibNqZvpqToZ+8yJGzDcvYs3LgBe/bol6wqVjRN2PfcAzVrQrVq4OlpdmwiZxplq6aODi45ORkfHx+SkpLw9va2dzhCCKxr4WvuPunp6Zw3PKfNZXpIc7fZuXMnK1eu5MUXX6RevXq5Hjc9PZ2YmBi8vb3x9fXNtp3heJlrwzndfh8/fjw7d+6kVatWjB8/nuTkZBITE3N9lpv5PC5evMjMmTNt88w3NRXOnCHj6FGSdu2iXHw82lOn4MQJyNJAL5uAAKhWDVWtGqmVKuFapw7aGjX0/b4rV9YPm1oKFCT3SA1aCGE31tREi7pWb+iLPWLECEJDQ20y1GnW2nBWmWv/mbtwmVMTtdUzX51OR/Lt23jXro2qVYvkZs3wrlwZreGLTHIynDypT9YnTuhfHz8OZ85AQgLEx0N8PJodO8jWZt3FBapU0de0q1fXL1Wr6mv3oaEQGKgfRrWUkwQthCi1zEmmOfXFzk1+LZYzD2iSkpKS6+13S28LF8at/czdo3I8L29vaN5cv2SVmAjR0fpkbfjXsJw9q59D+9Qp/ZITFxd9LduQsENC/nttWLy8CnR+xYEkaCFEqWXr2nh+tVdLa8PmynwehfHM1+Jaebly0LSpfskqI0Pfijxz0j5zRt9QLSYGLlzQJ/DoaP2S1zEMybpyZahUCYKD//s3OFjfdawALdDtTRK0EKJEskcLdFvVXh2h9XzmWrNNW2I7OelrxCEh0K5d9vfT0/UJPCZGv8TG/vfasCQm/rccOJD7sdzc/kvWWZN35tf/NmhzNNJITAgh8mHLhmmWlFeYMeQnc0O1t956yyZl2kxy8n+J+9w5fa3b0L/b8PraNfPL8/TUt0jPvAQEZF9XsaL+1roFtXJpJCaEEIXI1rfCC7NxnK1q3zk1VHOY/vDe3vp+3vXr577NnTv6luaZk3fmBG5YbtzQdzNLSYHTp/M/trt7zgk8IAAqVPhv8ffX/1sAUoMWQggbsVXt1RGUpHPJ040b+kRuGKgl6xIf/99rMxoKZpXs6YlPSorUoIUQwp4crqYp8uflpV9q1cp/21u3ck7chp+vXIGrV//7Nz3dqqRuIAlaCCFsxBFGXrMV+bKRAw8Pfd/tatXy31ap/7qb5dQVzQySoIUQQmRTkr5s2IVGA+XLF2jENPlqJIQQQjggSdBCCCGEA5IELYQQQjigYpWg582bR9WqVXF3d6d169bs3LnT3iEJIYQQhaLYJOhVq1YRGRnJ+PHj+eeff2jcuDGdO3cmPj7e3qEJIYQQNldsEvSMGTMYOHAg4eHh1KtXjwULFuDh4cGSJUvsHZoQQghhc8Wim9Xdu3fZs2cPY8eONa7TarV07NiRbdu25bhPamoqqampxp+TkpIA/YhiQgghRFEw5BxDDjJwc3PDzc0tz32LRYK+evUqGRkZVKxY0WR9xYoVOXbsWI77TJs2Lcdp0UJCQgolRiGEECI3WecHHz9+PO+8806e+xSLBG2NsWPHEhkZafxZp9PRvHlz/vnnHzRmzkTSsmVLdu3alec2ycnJhISEEBsbK2N8Z2HO9bO3oo6xMI9nq7ILWo41+1u6j3w2C0Y+m0V3TKUUzZo1Y8+ePSajsuVXe4ZikqD9/f1xcnIiLi7OZH1cXByBgYE57pPT7QM3Nzd8fHzMPq6Tk5PZH2xvb2/5I5CFJdfPXoo6xsI8nq3KLmg51uxv6T7y2SwY+WwW7THd3NysGpWtWDQSc3V1pXnz5mzYsMG4TqfTsWHDBtq0aWN2OUOGDLHouJZuL0wVh+tX1DEW5vFsVXZBy7Fmf/lsFq3icP3sEWNhHdPacovNdJOrVq0iLCyMTz75hFatWjFr1ixWr17NsWPHsj2bLkoyjaUQjkk+m6K4Kxa3uAH69OnDlStXePvtt7l8+TJNmjRh3bp1dk3OoL91MX78eLOeJwghio58NkVxV2xq0EIIIURpUiyeQQshhBCljSRoIYQQwgFJghZCCCEckCRoIYQQwgFJgi5kP/74I7Vr16ZWrVosWrTI3uEIIf7Vs2dPypcvz5NPPmnvUITIkbTiLkTp6enUq1ePjRs34uPjQ/Pmzfn777/x8/Ozd2hClHqbNm3ixo0bfPbZZ3z11Vf2DkeIbKQGXYh27txJ/fr1qVSpEp6ennTp0oX169fbOywhBNC+fXu8vLzsHYYQuZIEnYctW7bw+OOPExwcjEaj4dtvv822zbx586hatSru7u60bt2anTt3Gt+7ePEilSpVMv5cqVIlLly4UBShC1GiFfSzKURxIAk6Dzdv3qRx48bMmzcvx/dXrVpFZGQk48eP559//qFx48Z07tyZ+Pj4Io5UiNJFPpuiNJAEnYcuXbowefJkevbsmeP7M2bMYODAgYSHh1OvXj0WLFiAh4cHS5YsASA4ONikxnzhwgWCg4OLJHYhSrKCfjaFKA4kQVvp7t277Nmzh44dOxrXabVaOnbsyLZt2wBo1aoVhw4d4sKFC6SkpPDLL7/QuXNne4UsRKlgzmdTiOKg2EyW4WiuXr1KRkZGtsk6KlasyLFjxwBwdnbmww8/5KGHHkKn0zF69GhpwS1EITPnswnQsWNH9u/fz82bN6lcuTJr1qyxaPpaIQqbJOhC1q1bN7p162bvMIQQWfz+++/2DkGIPMktbiv5+/vj5OREXFycyfq4uDgCAwPtFJUQQj6boqSQBG0lV1dXmjdvzoYNG4zrdDodGzZskNtkQtiRfDZFSSG3uPOQkpLCqVOnjD9HR0ezb98+fH19CQ0NJTIykrCwMFq0aEGrVq2YNWsWN2/eJDw83I5RC1HyyWdTlApK5Grjxo0KyLaEhYUZt5kzZ44KDQ1Vrq6uqlWrVmr79u32C1iIUkI+m6I0kLG4hRBCCAckz6CFEEIIByQJWgghhHBAkqCFEEIIByQJWgghhHBAkqCFEEIIByQJWgghhHBAkqCFEEIIByQJWgghhHBAkqCFEEIIByQJWgghhHBAkqCFEEIIByQJWgghBLGxsbRv35569erRqFEj1qxZY++QSj2ZLEMIIQSXLl0iLi6OJk2acPnyZZo3b86JEycoW7asvUMrtaQGLeyiffv2jBgxwt5hFImSdq72PJ9r164REBDA2bNn7XL8kiwoKIgmTZoAEBgYiL+/PwkJCcb3n376aT788EM7RVc6SYIWebpy5QqDBw8mNDQUNzc3AgMD6dy5M1u3bgWKR/LRaDR5Lu+8845Z5RTmufbr148ePXpYtI+9rv0333zDpEmTbB5Du3bt6N+/f57bTJkyhe7du1O1atUCHcua612a7Nmzh4yMDEJCQozrxo0bx5QpU0hKSrJjZKWLs70DEI6tV69e3L17l88++4zq1asTFxfHhg0buHbtmr1DM9ulS5eMr1etWsXbb7/N8ePHjes8PT3tEVax5evra/MylVLs3buX3r1757rNrVu3WLx4Mb/++qvNj18U7t69i6urq11jaNKkCenp6dnWr1+/nuDgYAASEhLo27cvCxcuNNmmQYMG1KhRgy+++IIhQ4YUSbylnhIiF9evX1eA2rRpU47vh4WFKcBkiY6OVkopdefOHfXKK6+oChUqKDc3N9W2bVu1c+dO477t2rVTw4cPN/78448/Km9vb/XFF18opZTKyMhQU6dOVVWrVlXu7u6qUaNGas2aNSbHX7NmjWrQoIFyd3dXvr6+qkOHDiolJSXPc4qKilI+Pj45vpdXzHmd6y+//KLatm2rfHx8lK+vr+ratas6depUruea27Xs3r27yT6vvPKKGjVqlCpfvryqWLGiGj9+vMn2OcWT33XLr1xzrqvhfHKL4bPPPlO+vr7qzp07JuV2795dPf/88zme//HjxxWg/v7771yv0Zo1a1SFChVM1rVr104NHTpUDR8+XJUrV04FBASoTz/9VKWkpKh+/fopT09PVaNGDfXzzz/ner2tiVcp837vQ4YMUcOHD1d+fn6qffv2Sin9/+3p06erGjVqKFdXVxUSEqImT55sPMe8rn1+v9+8yjbHnTt31AMPPKCWLVuW4/sTJkxQ999/v9nliYKRBC1ylZaWpjw9PdWIESOy/fFSSqnExETVpk0bNXDgQHXp0iV16dIllZ6erpRSatiwYSo4OFj9/PPP6vDhwyosLEyVL19eXbt2TSllmrSWL1+uvLy81A8//GAse/LkyapOnTpq3bp16vTp0yoqKkq5ubkZvyxcvHhROTs7qxkzZqjo6Gh14MABNW/ePHXjxo08zymvBJ1XzHmd61dffaW+/vprdfLkSbV37171+OOPq4YNG6qMjIxs55qbnBK0t7e3euedd9SJEyfUZ599pjQajVq/fn2e1z6/65ZfueZcV8P55BbDrVu3lI+Pj1q9erVxn7i4OOXs7Kz++OOPHM9/xYoVysnJSd28eTPXazRs2DD16KOPmqxr166d8vLyUpMmTVInTpxQkyZNUk5OTqpLly7q008/VSdOnFCDBw9Wfn5+JmVnvt7WxKuUeb93T09PNWrUKHXs2DF17NgxpZRSo0ePVuXLl1dLly5Vp06dUn/++adauHChWdc+v99vbmWbQ6fTqaeffjrbF7bMfvnlF+Xq6prj3wNhe5KgRZ6++uorVb58eeXu7q7uu+8+NXbsWLV//37j+zkln5SUFOXi4qKWL19uXHf37l0VHBys3nvvPZP95s6dq3x8fExq6Xfu3FEeHh7ZalMRERHqmWeeUUoptWfPHgWos2fPWnQ+uSVoS2LOz5UrVxSgDh48aPZ+OSXorDWVli1bqjFjxphsk7lcc65bfuWac10zHze3cxs8eLDq0qWL8ecPP/xQVa9eXel0uhzLHDlypKpXr16ux1RKX6Pt379/tlgyn096eroqW7aseuGFF4zrLl26pAC1bds247qs19vSeHOS0++9adOmJtskJycrNze3HJNmftc+v99vXmWb488//1QajUY1btzYuBw4cMBkm/3791v1uRPWkWfQIk+9evWia9eu/Pnnn2zfvp1ffvmF9957j0WLFtGvX78c9zl9+jRpaWm0bdvWuM7FxYVWrVpx9OhR47qvvvqK+Ph4tm7dSsuWLY3rT506xa1bt3jkkUdMyr179y5NmzYFoHHjxnTo0IGGDRvSuXNnOnXqxJNPPkn58uWtOk9zY87JyZMnefvtt9mxYwdXr15Fp9MBEBMTQ4MGDayKB6BRo0YmPwcFBREfH5/r9uZct/zKtdV1HThwIC1btuTChQtUqlSJpUuX0q9fPzQaTY7b//PPPzRv3jzPMm/fvo27u3u29ZnPx8nJCT8/Pxo2bGhcV7FiRYA8r52l8YJ5v/es53T06FFSU1Pp0KFDtvLyu/b5/X7zKtsc999/v/EcclOmTBlA3x5AFD5J0CJf7u7uPPLIIzzyyCO89dZbDBgwgPHjx+eaoM3VtGlT/vnnH5YsWUKLFi2MfwxTUlIA+Omnn6hUqZLJPm5uboD+D/Fvv/3G33//zfr165kzZw5vvvkmO3bsoFq1agWKy1KPP/44VapUYeHChQQHB6PT6WjQoAF3794tULkuLi4mP2s0mjz/gJpz3fIr11bXtWnTpjRu3Jhly5bRqVMnDh8+zE8//ZTr9v/88w+PP/54nmX6+/tz/fr1bOtzOp/M6wz/r/K6dpbGC+b93rP2ITYkuJzkd+3z+/0mJibmGa8tGLpdVahQodCPJaSblbBCvXr1uHnzJgCurq5kZGSYvF+jRg1cXV2NXbEA0tLS2LVrF/Xq1TPZbuPGjXz33Xe88sorJuW7ubkRExNDzZo1TZbM3T40Gg1t27ZlwoQJ7N27F1dXV9auXWvVOZkTc07neu3aNY4fP864cePo0KEDdevWzTGJFIas8Zh73fJjyXXN6ZoYDBgwgKVLlxIVFUXHjh1zjeHMmTMkJibmW4Nu2rQpR44cMfs8LGVuvGD9771WrVqUKVOGDRs25Ph+Xtc+v99vfmXbwqFDh6hcuTL+/v6FdgzxH6lBi1xdu3aN3r17079/fxo1aoSXlxe7d+/mvffeo3v37gBUrVqVHTt2cPbsWTw9PfH19aVs2bIMHjyYUaNG4evrS2hoKO+99x63bt0iIiLC5Bj33HMPGzdupH379jg7OzNr1iy8vLwYOXIkr776Kjqdjvvvv5+kpCS2bt2Kt7c3YWFh7Nixgw0bNtCpUycCAgLYsWMHV65coW7duladqzkx53Su5cuXx8/Pj08//ZSgoCBiYmJ4/fXXC3bhzZRTPPldt/xYel1zikGr1X/vf/bZZxk5ciQLFy5k2bJluR5zz549aDQa4yAZuencuTNjx47l+vXrVj/KyIu58QJW/97d3d0ZM2YMo0ePxtXVlbZt23LlyhUOHz5MgwYN8rz25nwucis76+fOWn/++SedOnWySVnCDPZ+CC4c1507d9Trr7+umjVrpnx8fJSHh4eqXbu2GjdunLp165ZSSt895t5771VlypQx6Xp0+/Zt9corryh/f3+zulkdOXJEBQQEqMjISKWUvkXprFmzVO3atZWLi4uqUKGC6ty5s9q8ebNx+86dOxu7RN1zzz1qzpw5+Z5TXq2484s5t3P97bffVN26dZWbm5tq1KiR2rRpkwLU2rVrczzXnOTUSCzrPt27d1dhYWF5xpPfdcuvXHOua+YycrsmBi+88EKOXZgye/3119U999yT5/UxaNWqlVqwYEGOsRhUqVJFzZw502Rd5t+HIa5evXplK9+ceA2s/b1nZGSoyZMnqypVqigXFxcVGhqqpk6data1z+/3m1vZtnD79m3l4+Nj0thOFC4Zi1sIUWg6dOhA/fr1mT17tk3K++mnnxg1ahSHDh0y1tSt8eijj1KzZk3mzp1rst7W8ZYk8+fPZ+3ataxfv97eoZQacotbCGFz169fZ9OmTWzatImPP/7YZuV27dqVkydPcuHCBYueq2eOa+vWrWzatIlBgwYVerwliYuLC3PmzLF3GKWK1KCFEDZXtWpVrl+/zltvvcXIkSPtHY5Rz5492bVrF2FhYUyePNnYwttR4xWlmyRoIYQQwgFJNyshhBDCAUmCFkIIIRyQJGghhBDCAUmCFkIIIRyQJGghhBDCAUmCFkIIIRyQJGghhBDCAUmCFkIIIRyQJGghhBDCAUmCFkIIIRyQJGghhBDCAUmCFkIIIRzQ/wGvPv4PhVOB8wAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "PvI_plot_Step2 = HtoH.PvI(CatP, CatdP, CatI, Iscale=[0.8,50.0], Pscale=[0.0,30.0], showfit='true', weighted='false', errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "59006181",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Power law index: -1.0895886309264438 ± 0.10703992344817052\n",
      "Coefficient: 1.9375896420061683e+24 ± 1.0361457477407953e+25\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAegAAAF2CAYAAAClPAs/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABjo0lEQVR4nO3deVhU5dvA8e8MqyigiAokqOCeivuSmmam2aKmmdnPDaEVTcWlaNHUjLJFRU3fSlFL0yy1RVuMXMpccs3dVAyUxRQBQUVgzvvHaSaGfYaBGeD+XNe5ZM6ceeaeQeaeZ9coiqIghBBCCJuitXYAQgghhMhPErQQQghhgyRBCyGEEDZIErQQQghhgyRBCyGEEDZIErQQQghhgyRBCyGEEDZIErQQQghhgyRBCyGEEDaoyiTozMxM3njjDTIzM60dihCiEpLPGGFpNpGgly5dSps2bXBzc8PNzY1u3brx/fffG+6/ffs2oaGh1K5dmxo1ajB06FCSkpJMeo7MzExmzZolfzxCiDIhnzHC0mwiQdevX5+3336bgwcPcuDAAfr06cOgQYM4ceIEAJMnT+bbb79lw4YN7Ny5k/j4eIYMGWLlqAXAkiVLrB1CqdnSayjPWMryuSxVtiXKMbcMW/p/UZVVht+D2a9BsVG1atVSPvnkEyUlJUVxcHBQNmzYYLjv1KlTCqDs2bOnxOWlpqYqgJKamloW4VZZLVq0sHYIpWZLr6E8YynL57JU2ZYox9wyTH2cfMaUDVv6+zSXua/B3oJfEiwiJyeHDRs2kJGRQbdu3Th48CBZWVn07dvXcE3z5s3x8/Njz549dO3atcByMjMzjZqaUlJSAEhNTS3T+KuanJwc0tLSrB1GqdjSayjPWMryuSxVtiXKMbcMUx+nv9ZW/i9VFrb092kORVHIzMwkJSUFrfa/RmsnJyecnJyKfbBN+PPPP5Xq1asrdnZ2iru7u7JlyxZFURRlzZo1iqOjY77rO3XqpEyfPr3Q8mbOnKkAcsghhxxyyGFzx8yZM4vNizZTg27WrBlHjhwhNTWVL7/8kjFjxrBz506zywsPDycsLMxwOzU1FT8/P+Li4nBzc7NEyEIIIUSR0tLS8PX1JTY2Fnd3d8P5YmvPgM0kaEdHRxo3bgxAhw4d+OOPP1i4cCHDhw/nzp07pKSkULNmTcP1SUlJeHl5FVpeYc0H+pHiQgghRHlxd3c3OffYxCjuguh0OjIzM+nQoQMODg5ER0cb7jtz5gyxsbF069bNihEKIYQQZccmatDh4eEMGDAAPz8/bty4wdq1a9mxYwc//vgj7u7uBAcHExYWhoeHB25ubkyYMIFu3boVOkBMCCGEqOhsIkFfuXKF0aNHk5CQgLu7O23atOHHH3/kgQceAGD+/PlotVqGDh1KZmYm/fv358MPP7Ry1EIIIUTZ0SiKolg7iPKQlpaGu7s7qamp0gcthBCiXJQm99hsH7QQQghRlUmCFkIIIWyQJGghhBDCBkmCFkIIIWyQJGghhBDCBkmCFkIIIWyQJGghhBDCBkmCFkIIIWyQJGghhBDCBkmCFkIIIWyQJGghhBDCBkmCFkIIIWyQJGghhBDCBkmCFkIIIWyQJGghhBDCBkmCFkIIIWyQJGghhBDCBkmCFkKIKmLlypVoNBo0Go3hXO/evdFoNIwdO9Z6gYkCSYIWQohKQJ9oCzo2b94MQJ06dejSpQtdunSxbrAWNHbsWDQaDb179y7R9RMmTCAwMBB7e3s0Gg1eXl5lG2Ap2Fs7ACGEEJbj6OhIu3btjM55eHgA8PDDD/Pwww9bIyyb8emnn+Lo6IiHhwf//POPtcMpktSghRCiEvH29mbv3r1Gx7333gsU3MRdkDFjxtCkSRNcXV1xdHSkQYMGvPjii6SlpRmu0ddcGzZsyKpVq/D19cXNzY1JkyZx69YtJk2ahLu7O35+fixbtsyo/Pj4eMaNG4ePjw+Ojo74+/szZ84csrOzDdfoWwRGjx7NzJkz8fb2platWowcOZIbN24AGJ4bYOfOnYbXtmPHjkJf27Fjx7hy5QoPPfSQSe+rNUgNWgghhJGvv/4ae3t7AgICuHHjBhcuXGDRokUkJCSwYcMGo2sTEhJ44YUX8Pb25saNGyxcuJCffvqJy5cv4+7uTlxcHKGhofTq1YsWLVpw7do1unbtSlxcHK6urrRo0YKTJ08yY8YMYmJiWLFihVH569atw9nZGU9PTxITE1mzZg0NGjRg7ty5tGvXjoyMDK5evYqrqystW7YEwM3NrdDX5uvra/k3rIxIDVoIISqRv//+O18ftKl27tzJ1atXOXLkCOfPn+fVV18FYPPmzdy+fdvo2jt37vDTTz9x9uxZQ/I7d+4cR48e5eTJkzg7O6PT6di5cycAixcvJi4ujnr16nH+/HmOHj3Kl19+Cag1/HPnzhmV7+zszKlTpzh37hwdOnQAIDo6GoBNmzYZmuzbt29vaDFo3769ya/ZFkkNWgghKpGC+qBN9fPPP/O///2P8+fPGyXk7Oxs/vnnH6NaaK1atejevTsAfn5+xMXF0apVKxo2bAhA3bp1iY2NJSkpCYD9+/cDkJSURN26dY2eV1EU9u3bR+PGjQ3n+vTpw1133QVA8+bNOXjwoKGsyk4StBBCVCL6PmhzrVmzhqlTpxrK8vX15erVq1y4cAGAnJwco+tzNyfb29vnO6evwSuKYvS43E3Subm4uBjdrlmzZr7y85ZVWUkTtxBCCAN9cnd1dSUmJoZ9+/bRr18/i5XfqVMnQE2269atMzRLb9u2jRdeeIHHHnvMpPL0CT0jI8NiMdoKqUELIYQwaNOmDQA3btzA398fJycnUlNTLVZ+aGgon3zyCZcvX6ZZs2a0aNGCGzduEBcXR1ZWFqNHjzapvObNmwNw4MABWrduTfXq1dm+fTvVqlUr8PrevXtz6dIlrly5AsDVq1cNTepr1qyxqTniUoMWQghhEBwcTFhYGJ6enty4cYPevXsze/Zsi5Vfp04d9u7dS1BQELVr1+bEiRPcunWLnj17Mn/+fJPLGzduHEOHDsXd3Z3jx4+zb9++fM3wuV28eJHz588bpmrl5ORw/vx5zp8/z61bt8x+XWVBo1SRxvy0tDTc3d1JTU0tcgi+EEIIYSmlyT1SgxZCCCFskCRoIYQQwgZJghZCCCFskCRoIYQQwgZJghZCCCFskCRoIYQQwgZJghZCCCFskCRoIYQQwgZVvQT9++/WjkAIIYQolk0k6IiICDp16oSrqyt169Zl8ODBnDlzxuia3r1759vj9LnnnjP9yd5/30JRCyGE7dq1axePPvooPj4+aDQaNm/eXKLHLVmyhIYNG+Ls7EyXLl0M20Pq3bhxg0mTJtGgQQOqVavGPffcwx9//FEGr6D4WPIqyWs2932xBptI0Dt37iQ0NNSwo0lWVhb9+vXLtzvJ008/TUJCguGYN2+e6U/2889w+LCFIhdCCNuUkZFBYGAgS5YsKfFj1q9fT1hYGDNnzuTQoUMEBgbSv39/w8YSACEhIWzbto1PP/2UY8eO0a9fP/r27cvly5eLLb93796sXLnSYrHkVZLXbM77YjWKDbpy5YoCKDt37jSc69WrlzJx4kSzy0xNTVUAJRUUZdgwC0QphKiKmjRponTt2lW5efOm4ZxOp1O6dOmivPzyy1aMrHCAsmnTpmKv69y5sxIaGmq4nZOTo/j4+CgRERGKoijKzZs3FTs7O+W7774zelz79u2VV199tdjye/XqpURFRZUo5uJiKU5JXnNJ35fSMOSe1FSTH2sTNei89FubeXh4GJ1fs2YNnp6etGrVivDwcG7evFloGZmZmaSlpRkdBl9+CXma0IUQoiTWr1/PoUOH2L17t+HcmjVr+Pvvv3nllVcM59566y1q1KhR5BEbG2uNl1CgO3fucPDgQfr27Ws4p9Vq6du3L3v27AEgOzubnJwcnJ2djR5brVo1fvvtt3KNpaLJm48yMzOLfYzNJWidTsekSZPo3r07rVq1Mpx/6qmn+Oyzz9i+fTvh4eF8+umnjBw5stByIiIicHd3Nxy+vr7qHQ8+CIoC5jSPCyGqvHbt2tG2bVtOnz4NwM2bNwkPD+fNN9/E1dXVcN1zzz3HkSNHijx8fHys9TLyuXr1Kjk5OdSrV8/ofL169UhMTATA1dWVbt26MWfOHOLj48nJyeGzzz5jz549JCQk5Csz75eUX3/9leeee67YLykliaWi8fX1NcpJERERxT7GvhziMkloaCjHjx/P923smWeeMfzcunVrvL29uf/++zl//jwBAQH5ygkPDycsLMxwOy0tTU3SU6fCDz/A6tXwxhugT9xCCFFCTZs2NQxknTdvHp6engQFBRld4+Hhka8V0FQvv/wy77zzTpHXnDp1iubNm5fqeUzx6aefMm7cOO666y7s7Oxo3749I0aM4ODBg/mufe6553jiiScMt//3v/8xdOhQhgwZYjhnS19SylJcXJzRdpNOTk7FPsamEvT48eP57rvv2LVrF/Xr1y/y2i5dugBw7ty5AhO0k5NTwW9Ap07Quzfs2KGO6F6woPSBCyGqlGbNmrFr1y4uXbrEu+++y5YtW9BqjRsk33rrLd56660iyzl58iR+fn6F3j9lyhTGjh1bZBn+/v4ljrsonp6e2NnZkZSUZHQ+KSkJLy8vw+2AgAB27txJRkYGaWlpeHt7M3z48ALjyPslpVq1atStW5fGjRtbJJaKxM3NzeT9oG0iQSuKwoQJE9i0aRM7duygUaNGxT7myJEjAHh7e5v+hK+8oibojz6CV1+FOnVML0MIUWU1bdqUjz/+mJdffpl+/frRu3fvfNfkrT0WpLjaY506dahTTp9Pjo6OdOjQgejoaAYPHgyoXY7R0dGMHz8+3/XVq1enevXqXL9+nR9//NG8WTUWiqWysokEHRoaytq1a/n6669xdXU19DG4u7tTrVo1zp8/z9q1a3nooYeoXbs2f/75J5MnT+bee++lTZs2pj9h377QoQMcPAiRkTBnjoVfkRCiMmvatClxcXF8+eWXHD9+vMBrLNHEXRrp6emcO3fOcDsmJoYjR47g4eFhqLUvXryYTZs2ER0dDUBYWBhjxoyhY8eOdO7cmQULFpCRkWHUfP/jjz+iKArNmjXj3LlzTJs2jebNm+dr4tfHkJ6ebri9bt06AKN+5Dp16mBnZ5fvsSWJJW/8JXnNJbnGZlh8TLkZgAIP/XD82NhY5d5771U8PDwUJycnpXHjxsq0adNMGraeb6j7V18pCiiKu7uimDH8XQhRdWVkZCgajUaZMmWKtUMp1Pbt2wv8XB0zZozhmpkzZyoNGjQwetyiRYsUPz8/xdHRUencubOyd+9eo/vXr1+v+Pv7K46OjoqXl5cSGhqqpKSkFBjDzJkzC/181x8xMTGFvobiYskbf0lec0musaTSTLPSKIqilMP3AKtLS0vD3d2d1NRUtR9Ap4O774bTp+Htt+Gll6wdohCigkhOTqZ27docPXrUvFY8UWXkyz0msLlpVuVGq4WXX1Z/nj8fbt2ybjxCiArj6NGjODo60qJFC2uHIiqxUiXorKws4uLiOHPmDMnJyZaKqfw89RT4+UFSEkRFWTsaIUQFcfToUVq2bImDg4O1QxGVmMkJ+saNGyxdupRevXrh5uZGw4YNadGiBXXq1KFBgwY8/fTTZbZwusU5OMC0aerP774LWVnWjUcIUSFMmjSJw7KmvyhjJiXoDz74gIYNGxIVFUXfvn3ZvHkzR44c4ezZs+zZs4eZM2eSnZ1Nv379ePDBB/nrr7/KKm7LCQ6GunXh4kX4d4ShEEIIYW0mDRIbMWIEr732GnfffXeR12VmZhIVFYWjoyPjxo0rdZCWUGRHfUSEOje6ZUs4dkztnxZCCCFKqTSDxKruKO7cUlPVvui0NNi0Cf6dGC+EEEKUhtVHcV+9epUtW7bwzTffFLhgus1zd4fQUPXnt95SN9MQQgghrKjUCfqrr76icePGzJo1i5kzZxIQEEBURRwRPWkSODvDH3/AL79YOxohhBBVnMkJOveybQCzZs1i//797N+/n8OHD7NhwwZeffVViwVYburWhZAQ9ediFrgXQgghyprJCbpDhw58/fXXhtv29vZcuXLFcDspKQlHR0fLRFfepk4Fe3u1Br1vn7WjEUIIUYWZPEjs4sWLhIaG4ujoyJIlSzh//jxPPvkkOTk5ZGdno9VqWblyJQ899FBZxWyWEnfUjx0Lq1bBoEGweXN5hSeEEKISKtdBYg0bNmTLli088cQT9OrViyNHjnDu3Dm2bdvGzz//TGxsrM0lZ5O89BJoNPD113DihLWjEUIIUUWZPUhsxIgR/PHHHxw9epTevXuj0+lo27Ytzs7Oloyv/LVoAY89pv4s21AKIYSwErMS9NatW3n//fc5cOAAn3zyCfPmzeN///sf06ZN41Zl2HRixgz13/Xr4cgRq4YihBCiajI5QU+ZMoWgoCD++OMPnn32WebMmUOvXr04dOgQzs7OtGvXju+//74sYi0/gYHw5JPqz6+9Zt1YhBBCVEkmDxKrXbs2P/30Ex06dCA5OZmuXbty9uxZw/0nT57k2Wef5ddff7V4sKVhckf92bPq0p85OfDbb9C9e9kHKYQQolIp10Fi1atXJyYmBoC4uLh8fc4tW7a0ueRslqZNIShI/fmVV2R1MSGEEOXK5AQdERHB6NGj8fHxoVevXsypzAOpZswAJyfYtQt++sna0QghhKhCzNos49q1a1y4cIEmTZpQs2bNMgjL8sxuZggLg/nzoX17OHBAnYIlhBBClEC5NXHHxsYCaj90p06dik3Oly9fNikYmxQeDjVqwKFDsHGjtaMRQghRRZiUoDt16sSzzz7LH3/8Ueg1qampfPzxx7Rq1Yqvvvqq1AFaXZ06ai0a1BHdOTnWjUcIIUSVYG/KxSdPnmTu3Lk88MADODs706FDB3x8fHB2dub69eucPHmSEydO0L59e+bNm1exVxTLLSwMFi+G06fh00/V5UCFEEKIMmRWH/StW7fYsmULv/32G3///Te3bt3C09OTdu3a0b9/f1q1alUWsZZKafoBAHj3XZg+HRo0gDNn1MFjQgghRBFKk3vMStAVUakT9K1b0LgxxMfDokUwfrzlgxRCCFGplOs86CqrWjV4/XX15zffhIwM68YjhBCiUpMEbYpx48DfH5KSIDLS2tEIIYSoxCRBm8LREWbPVn+eNw+uX7duPEIIISotSdCmevJJaNUKUlLgvfesHY0QQohKyqwE/ffff/Ptt9+yf/9+S8dj++zs1D5ogAULIDHRquEIIYSonExO0J9//jlNmzZl0KBBdOvWjY4dO/LPP/+URWy2a+BA6NIFbt6Et96ydjRCCCEqIZMT9KxZs3jqqac4ffo0P/27gcTLL79s8cBsmkbzX2Jetgz+/tu68QghhKh0TJ4H7ejoyNmzZ2nYsCEAp0+fpkOHDmTY+LSjUs+DLkjfvhAdrW5LuWKFZcoUQghRaZTrPOjs7GxcXFwMt5s3b45OpyOxKvbFzp2r/rtqFZw4Yd1YhBBCVCpmDRJbtWoVv//+O+np6QDY29tz8+ZNiwZWIXTpAkOGgE4HkydD1ViUTQghRDkwOUH37NmTN998kx49elCzZk2aNGnC7du3Wb58Odu3b+fGjRtlEaftevdddX70tm3w3XfWjkYIIUQlYfZa3H/99RcHDx7k0KFDhiMlJQWtVkuTJk04deqUpWMtlTLpg9Z7+WV45x1o0gSOH1cTthBCiCrPZjbLiImJ4cCBAxw+fJi3bGz6UZkm6Bs3oGlTdU70e+/BlCmWLV8IIUSFZDMJ2paVaYIGiIpS1+p2c4O//oK6dS3/HEIIISqUCr+bVUREBJ06dcLV1ZW6desyePBgzpw5Y3TN7du3CQ0NpXbt2tSoUYOhQ4eSlJRkpYgLMGYMdOgAaWn/7XolhBBCmMkmEvTOnTsJDQ1l7969bNu2jaysLPr162c0t3ry5Ml8++23bNiwgZ07dxIfH8+QIUOsGHUeWq269CfAJ5/A0aNWDUcIIUTFZpNN3P/88w9169Zl586d3HvvvaSmplKnTh3Wrl3L448/DqgLpLRo0YI9e/bQtWvXfGVkZmaSmZlpuJ2Wloavr2/ZNXHrPfkkrF8PvXvDL7+oq44JIYSokvRN3HFxcUa5x8nJCScnpyIfaxM16LxSU1MB8PDwAODgwYNkZWXRt29fwzXNmzfHz8+PPXv2FFhGREQE7u7uhsPX17fsAwd1G0pnZ9ixAzZtKp/nFEIIYdN8fX2NclJERESxj7EvzRNGR0cTHR3NlStX0Ol0RvetMHPpS51Ox6RJk+jevTutWrUCIDExEUdHR2rWrGl0bb169QpdwSw8PJywsDDDbX0Nusz5+cG0aTBnDkydCg89pCZsIYQQVVZBNejimJ2gZ82axezZs+nYsSPe3t5oLNSUGxoayvHjx/ntt99KVU5Jmg/KzEsvqWtzx8So/dJVbTMRIYQQRtzc3EzuXjU7QS9btoyVK1cyatQoc4vIZ/z48Xz33Xfs2rWL+vXrG857eXlx584dUlJSjGrRSUlJeHl5Wez5LaZ6dXj7bRg1St07eswY8Pa2dlRCCCEqELP7oO/cucM999xjkSAURWH8+PFs2rSJX375hUaNGhnd36FDBxwcHIiOjjacO3PmDLGxsXTr1s0iMVjcU0+pa3VnZMArr1g7GiGEEBWM2Qk6JCSEtWvXWiSI0NBQPvvsM9auXYurqyuJiYkkJiZy69YtANzd3QkODiYsLIzt27dz8OBBgoKC6NatW4EjuG2CVgsLF6o/r1wJBw5YNRwhhBAVi9nTrCZOnMjq1atp06YNbdq0wcHBwej+Dz74oORBFNJ/HRUVxdixYwF1oZIpU6bw+eefk5mZSf/+/fnwww9L3MRd5iuJFWbUKPjsM+jeHX79VaZdCSFEFWKVpT7vu+++wgvVaPjll1/MKbbMWC1BX76srtN98yZ8/rk6T1pUaTqdjrS0NNzc3NBqbXKmoxDCQmQt7hKwWoIGdcrVjBlQvz6cOgU1apTv8wubkp2dzaVLl6hfvz729qWa6SiEsHEVfi3uSm/qVGjUCC5dgpkzrR2NEEKICqBUX99TUlJYvny5Ye/nli1bEhwcjLu7u0WCqzSqVYMlS9RFSxYsgJEjoV07a0clhBDChpldgz5w4AABAQHMnz+f5ORkkpOTmT9/PgEBARw6dMiSMVYOAwbA8OGg08Ezz0BOjrUjEkIIYcPM7oPu2bMnjRs35uOPPzb0o2VnZxMSEsKFCxfYtWuXRQMtLav2QeslJECLFpCaCpGRMGGCdeIQViV90EJUHVbpgz5w4AAvvfSS0QeMvb0906dP54DM+S2Yt7e6whjAq6+qfdJCCCFEAcxO0G5ubsTGxuY7HxcXh6ura6mCqtSeeQa6doUbN2DiRGtHI4QQwkaZnaCHDx9OcHAw69evJy4ujri4ONatW0dISAgjRoywZIyVi1YLH30E9vawcSN88421IxJCCGGDzO4Ae++999BoNIwePZrs7GwAHBwceP7553lb34wrCta6NUyZAu+8A+PHQ58+MjdaCCGEkVIvVHLz5k3Onz8PQEBAAC4uLhYJzNJsYpBYbjdvwt13w8WLEBYG779v7YhEOZFBYkJUHVZdqMTFxYXWrVvTunVrm03ONsnFBT78UP15wQI4fNiq4QghhLAtJn19DwsLY86cOVSvXp2wsLAirzVls4wqSz83ev16dfDY3r1gZ2ftqIQQQtgAkxL04cOHycrKMvxcmMJ2pxIFmD8ffvhB3Y7yww9lbrQQQgigFH3QsbGx1K9fP99uPIqiEBcXh5+fn0UCtBSb64PObdkyeP55cHWFkyfVTTVEpSV90EJUHVbpg27UqBFXr17Ndz45OZlGjRqZW2zVlHtu9IsvWjsaUQidTkdKSgo6nc7aoQghqgCzE3RhFe/09HScnZ3NDqhKyj03etMmtU9a2BxJ0EKI8mRy+5p+cJhGo2HGjBlGI7dzcnLYt28fbdu2tViAVUbr1vDKKzB7NrzwAvTqBV5e1o5KCCGElZicoPWDwxRF4dixYzg6Ohruc3R0JDAwkKlTp1ouwqrk1VfVlcWOHIFnn4XNm0EG3AkhRJVkcoLevn07AEFBQSxcuND2BlxVZI6OsHo1dOigJupPP4XRo60dlRBCCCswuw86KipKknNZaN0a3nhD/fnFF2XHKyGEqKLMTtARERGsWLEi3/kVK1bwzjvvlCqoKm/6dOjcWd03OiQESrcaqxBCiArI7AT9f//3fzRv3jzf+bvvvptly5aVKqgqz94eVq4EJyf48UdYvtzaEQkhhChnZifoxMREvL29852vU6cOCQkJpQpKAC1awNy56s9hYfD339aNRwghRLkyO0H7+vqye/fufOd3796Nj49PqYIS/5o0Cbp3VxcwGTcOZP6tEEJUGWavM/j0008zadIksrKy6NOnDwDR0dFMnz6dKVOmWCzAKs3OTm3qbtMGfvkFli6F0FBrRyWEEKIcmJ2gp02bxrVr13jhhRe4c+cOiqJQrVo1XnrpJV5++WVLxli1NW4M77yjjuiePh3691fPCSGEqNTM3ixDLz09nVOnTlGtWjWaNGmCk5OTpWKzKJveLKM4Oh307Qvbt0OPHrBjh2xLaQWW2uRCNssQouooTe4p9adDbGws165d486dO1y4cMFwfuDAgaUtWuhptbBihTpH+rff4P331dq0EEKISsvsBH3hwgUee+wxjh07hkajMWyeod8LOicnxzIRClXDhrBggTov+tVX4b77oFMna0clhBCijJg9invixIk0atSIK1eu4OLiwvHjx9m1axcdO3Zkx44dFgxRGIwbB8OGQXY2jBgBaWnWjkgIIUQZMTtB79mzh9mzZ+Pp6YlWq8XOzo4ePXoQERHBi7KncdnQaNRtKRs0gPPnZUS3EEJUYmYn6JycHFxdXQHw9PQkPj4egAYNGnDmzBnLRCfyq1kT1q5VB4l99pm6oYYQQohKx+wE3apVK44ePQpAly5dmDdvHrt372b27Nn4+/tbLEBRgHvu+W9DjRdegL/+smo4QgghLM/sBP3aa6+h+3dlq9mzZxMTE0PPnj3ZunUrkZGRFgtQFCI8HHr1gvR0tT/6zh1rRySEEMKCzB7F3b9/f8PPjRs35vTp0yQnJ1OrVi3DSG5RhvRN3IGBcPCgOrL73XetHZUQQggLMasGnZWVxf33389feZpWPTw8JDmXp/r1ISpK/fm999Sdr4QQQlQKZiVoBwcH/vzzT0vHIswxcCCMH6/+PHo0JCVZNx4hhBAWYXYf9MiRI1luoX2Kd+3axaOPPoqPjw8ajYbNmzcb3T927Fg0Go3R8eCDD1rkuSuFd99VVxm7cgXGjJFdr4QQohIwuw86OzubFStW8PPPP9OhQweqV69udP8HH3xQ4rIyMjIIDAxk3LhxDBkypMBrHnzwQaL0zblgs2t+W4WzM6xbBx07qs3cb78Nr7xi7aiEEEKUgtkJ+vjx47Rv3x6As2fPGt1naj/0gAEDGDBgQJHXODk54eXlZVqQVUnLlrB4MQQHw+uvq8uAPvCAtaMSFqbT6UhLS8PNzQ2t1uwGMCFEBWB2gt6+fbsl4yjWjh07qFu3LrVq1aJPnz68+eab1K5du9DrMzMzyczMNNxOqwrLYo4bB3v2wCefqFOvDh0CPz9rRyUsSKfTkZKSQo0aNSRBC1GB5M1BTk5OxbYEm/wXfuHCBUq5Q6XJHnzwQVavXk10dDTvvPMOO3fuZMCAAUVuyBEREYG7u7vh8PX1LceIrWjRIujQAa5dg8cfh1xfUoQQQliHr6+vUU6KiIgo9jEm7wdtZ2dHQkICdevWBWD48OFERkZSr14986LOG5BGw6ZNmxg8eHCh11y4cIGAgAB+/vln7r///gKvKagG7evrWzH3gzbVxYtqkk5Ohueeg6VLrR1RpWAL+0HLXtJCVCz6/aDj4uKMck+Z1KDz5vOtW7eSkZFhajGl4u/vj6enJ+fOnSv0GicnJ9zc3IyOKqNhQ1izRt1cY9kyWL3a2hGJMqBv7tbJqH0hbF7efFSSgc4VshPr0qVLXLt2DW9vb2uHYrsefBBmzlR/fvZZ+HfddFF5SIIWonIzOUHr5yHnPVca6enpHDlyhCNHjgAQExPDkSNHiI2NJT09nWnTprF3714uXrxIdHQ0gwYNonHjxkbLjYoCvP46DBgAt2/DkCFw/bq1IxJAfHw8CxcuNOwAJ4QQBTG5E0tRFMaOHWuont++fZvnnnsu3zzojRs3lrjMAwcOcN999xluh4WFATBmzBiWLl3Kn3/+yapVq0hJScHHx4d+/foxZ84cmQtdHK1WXa+7Qwe4cEFdaezrr9XzwmqWL1/O6dOnWb58ObNmzbJ2OEIIG2Vygh4zZozR7ZEjR5Y6iN69exc5MvxHWWPafB4e8NVX6haV330Hc+eqNWthNcHBwaSlpREcHGztUIQQNszkBJ17NS9RQbRvDx9+qC5iMmMGtGoFjz1m7aiqLB8fHyZOnIiPj4+1QxFC2DBp66wqxo2DCRPUn0eOhH/7+4UQQtgmSdBVyQcfqMt/3ryp7oIlO18JIYTNkgRdldjbw/r10LQpxMWpzdyy0pgQQtgkSdBVTa1a8O23ULOmum73M89AOS/dKoQQoniSoKuipk3hiy/Azk5dZey996wdUYUg85eFEOWpVIv5RkdHEx0dzZUrV/KtZrRixYpSBSbK2AMPwIIF6sCxl16CFi3gkUesHZVNk/nLQojyZHYNetasWfTr14/o6GiuXr3K9evXjQ5RAYSGqsuAKoq6PeXx49aOyKYFBwfTvHlzmb8shCgXZtegly1bxsqVKxk1apQl4xHlSaNRt6c8cwZ27ICHH1b7pWV+boFk/rIQojyZXYO+c+cO99xzjyVjEdbg4ABffqn2S8fGqkn6xg1rRyWEEFWe2Qk6JCSEtWvXWjIWYS21a8P330PduuoCJo8/DllZ1o5KCCGqNLObuG/fvs1HH33Ezz//TJs2bXBwcDC6/4MPPih1cKIc+fvDli3Qqxf89JM6/WrFCrUZXAghRLkzO0H/+eeftG3bFoDjeQYXlXb7SWElHTuq068GDoSVK6FBA3jjDWtHJYQQVZLZCXr79u2WjEPYiocfhqVL1dHds2aBr6+6yYYQQohyJQuViPyeeQZefVX9+dln4YcfrBuPEEJUQaVaqCQlJYXly5dz6tQpAFq2bElwcDDu7u4WCU5Y0Zw56nrdq1erg8a2b4dOnawdlRBCVBlm16APHDhAQEAA8+fPJzk5meTkZObPn09AQACHDh2yZIzCGjQa+Phj6NsXMjLgwQfhxAlrRyVykaVHhajczE7QkydPZuDAgVy8eJGNGzeyceNGYmJieOSRR5g0aZIFQxRW4+gIGzdC586QnAz9+kFMjLWjEv/KvfSoTqcjJSUl35K7QoiKq1Q16Jdeegl7+/9aye3t7Zk+fToHDhywSHDCBri6qnOk774b4uPVGnVCgrWjEhgvPSoJWojKx+wE7ebmRmxsbL7zcXFxuLq6liooYWM8PNS50f7+cOGCWpNOTrZ2VJVaYQk3d7O2LD0qROVmdoIePnw4wcHBrF+/nri4OOLi4li3bh0hISGMGDHCkjEKW+DjA9u2gbe3uqnGQw9Berq1o6q0CkvQuZu1hRCVm9mjuN977z00Gg2jR48mOzsbAAcHB55//nnefvttiwUobIi/v1qTvvde2LcPBg9WVx9zcrJ2ZFVGcHAwaWlpsqOWEFWA2TVoR0dHFi5cyPXr1zly5AhHjhwxjOR2kg/syqtVK7VPunp1iI6GYcPgzp0yearCapEVsb/VUjFLs7YQVUepFypxcXGhdevWtG7dGhcXF0vEJGxdly7wzTfg7AzffgvDh5fJ5hqSoMuercYlhDCxiTssLIw5c+ZQvXp1wsLCirxWNsuo5Pr0ga+/Vtft3rwZRoyAzz9Xt68UFYY+QdeoUQOtVhYWFMKWmJSgDx8+TNa/NaXDhw8Xep1sllFF9OunzpN+7DH46isYNQo++wzsS7VAnRBCCExM0Lk3yFi1ahX169fP961bURTi4uIsE52wfQ89pCbnIUNg/Xqws1OXB7Wzs3ZkQghRoZndptWoUSOuXr2a73xycjKNGjUqVVCignnkEdiwQa05r10LQUGQk1OmT1nefafSVyuEKG9mJ2hFUQo8n56ejrOzs9kBiQpq0KD/atCffgpjx8K/0+/KgqIopKamlmuCTk1NLfT/vRBCWJrJnYX6wWEajYYZM2YYjdzOyclh3759tG3b1mIBigpkyBBYt04dMPbZZ3DrllqjdnS0dmRVik6nIy0tDTc3Nxn4JUQFZnKC1g8OUxSFY8eO4Zjrw9fR0ZHAwECmTp1quQhFxfL44+pI7ieeUPumhw5Vm7+lVaXcyMhsISoHkxO0fqBYUFAQCxcuxM3NzeJBiQpu0CB1nvTgwfDdd/Doo+pUrOrVrR2ZEEJUGGbPh4mKigLg5MmTxMbGcifPalIDBw4sXWSiYuvfH7ZuVZPzzz/DgAHqsqCykYoQQpSI2Qk6JiaGwYMHc+zYMTQajWHwjH4OdE4Zj+IVFcB996lrdw8YAL/+Cg88oCZtDw9rRyaEEDbP7A6qF198kUaNGnHlyhVcXFw4ceIEu3btomPHjuzYscOCIYoK7Z571DW7a9VSN9i49164fNnaURlUlulTubehFEJUDmYn6D179jB79mw8PT3RarVotVp69OhBREQEL774oiVjFBVdx46wa5e6ZeWJE2rSPnvW2lEBlSdByzaUQlQ+ZifonJwcXP/tT/T09DR8c2/QoAFnzpyxTHSi8mjVCnbvhiZNIDYWuneHAwesHVWlERwcTPPmzWUbSiEqEbMTdKtWrTh69CgAXbp0Yd68eezevZvZs2fj7+9vsQBFJdKwIfz2G7RvD1evqn3U0dHWjqpAFa1mLdtQClH5mJ2gX3vtNcOH1+zZs4mJiaFnz55s3bqVyMhIk8ratWsXjz76KD4+Pmg0GjZv3mx0v6IozJgxA29vb6pVq0bfvn3566+/zA1dWFPdurB9u7obVnq6upb3hg3WjiofayfokvQpa7VaatasWaK5ztZ+PUII05mdoPv378+QIUMAaNy4MadPn+bq1atcuXKFPn36mFRWRkYGgYGBLFmypMD7582bR2RkJMuWLWPfvn1Ur16d/v37c/v2bXPDF9bk5qaO5h46FO7cURc1ef99kGU0ATWZLlmyhFOnThXZp2yJBC2Dy4SwXRZdZsjDw8OsrSYHDBjAm2++yWOPPZbvPkVRWLBgAa+99hqDBg2iTZs2rF69mvj4+Hw1bVGBODmpa3dPmKDenjoVXnyxzDfZqAh0Oh2DBg2iWbNmZd6nLIPLhLBdZifoPn36MGvWrHznr1+/bnINuigxMTEkJibSt29fwzl3d3e6dOnCnj17Cn1cZmYmaWlpRoewMXZ2sHChWnsGWLxYXc87I8O6cdmAevXq8eKLL5Z5n7IMLhOifOTNR5mZmcU+xuwEvWPHDhYvXszgwYPJyPWBeufOHXbu3GlusfkkJiYC6gdWbvXq1TPcV5CIiAjc3d0Nh6+vr8ViEhak0UBYmNoP7eSkLhF6332QlGTtyMqMTqcjPT3dJvqDZXCZEOXD19fXKCdFREQU+5hSNXH//PPPJCYm0rVrVy5evFiaoiwuPDyc1NRUwxEXF2ftkEQBDH2jQ4aoI7pr14Y//oAuXeDYMWuHVyZsKUHbUixCVGZxcXFGOSk8PLzYx5QqQXt7e7Nz505at25Np06dymQFMS8vLwCS8tSokpKSDPcVxMnJCTc3N6ND2B6jwUvdu8Pvv0NAAPz9N3b33ks1G52GVVlIghaifOTNR05OTsU+xuwErR8M5uTkxNq1a5k4cSIPPvggH374oblFFqhRo0Z4eXkRneuDOi0tjX379tGtWzeLPpewAU2bqkuC9u5NQno6S0NCSHjjjUo1wlur1drUVpA6nY7U1FRJ0kLYGLM/IZQ8H5ivvfYaa9as4X39gB8TpKenc+TIEY4cOQKoA8OOHDlCbGwsGo2GSZMm8eabb/LNN99w7NgxRo8ejY+PD4MHDzY3fGHLateGn37io3btOAasmDsXgoPVKVlWpCiKWbMUwHg6k60laEVRJEELYYPM/oSIiYnB09PT6NzQoUPZu3cvK1asMKmsAwcO0K5dO9q1awdAWFgY7dq1Y8aMGQBMnz6dCRMm8Mwzz9CpUyfS09P54YcfcHZ2Njd8YescHBi3aRONe/UiRKOBqCh1cZMiBgaWtdy7tpmqrKYzyQIkQlReZm832aBBgwLPt2rVilatWplUVu/evYv84NNoNMyePZvZs2ebVK6o2HzuuovnV67E68QJeOopdS3vjh3V+dNFjD+wJTqdjrS0NIKCgkhLS7P4dCZ9gpYvq0JUPiYl6LCwMObMmUP16tUJCwsr8toPPvigVIEJoaf076+O7B48GE6dwq5PH1xnzeLGsGHWDq1Y+gTq5eVV4ulM+oFbLi4u+e4zZfUwIUTFZlKCPnz4MFlZWYafC2NuP50QhWraFPbuhTFjSNi8mY/Dwxm+bx+sWQP2ZjcE2aRLly4RGRnJ2LFjadSokdF9+gQthKj8TPpk2759OwBZWVlotVqWLVtGkyZNyiQwIfJxc4OvvuL/7r+f4zt2wObNzLjvPvjiCyiky6UiioqK4vz582zYsEFmKghRhZnVTubg4MCff/5p6ViEKJ5WS/CqVdR54AFSHR1J2L8f2rWDb7+1dmRmyTvIS6fTMXToUAICAhhWAZrwhRBlx+yOrJEjR8oC+8IqfHx8cGrRglNdu/KRjw9cvw4DB8K0afBvF0xp5d3lKT4+nsjIyHwL5pSWTqcjOTmZ5ORkdDodOp0OFxcXxo4dm29526Lk7puWHaqEqBzMTtDZ2dksXbqUjh078uyzzxIWFmZ0CFGWhg8fTtO2bRm3axdMmqSefO89uPdeiI0tdfl5p0UtX76cM2fOsH79+lKXnTeBWmIecu4EbeqUrqtXr7J48WJJ6ELYGLMT9PHjx2nfvj2urq6cPXuWw4cPGw79giNClBXDbk8NGsD8+bBxI7i7qwPJ2rWD774rVfl5d3kKDg6mWbNmDB8+3Og6c+Yhl/UWj6buUPXNN99w5swZaRETwsaYPfxVP2BMCJvw2GPQti088QQcOACPPqo2ec+dCw4OJheXd5cnHx8fXnzxxXzX6RO0KSuDBQcHl8mcaD1TdqjSarU8+eSTfPvtt7LlpBA2RiZTisqjUSP47TfQJ9J334XevS3S5G1J1tjisbB+aa1Wi7+/v2w5KYQNKvUE0pMnTxIbG8udPOskDxw4sLRFi0pKv7pWmeww5uQECxdCr14wbpy6O1br1rBkCfzvf+r+0yWIz9rLZ2o0Gtzc3Cy2IEnuZvVZs2ZZpEwhRNkyO0FfuHCBxx57jGPHjhmtUaxfpCQnJ8cyEYpKJ3ezcJkZMkRt8n7qKXV3rFGj4JtvYNky8PAwXKavWU6ePBk/Pz9DfKmpqWUXWwloNBqLrhhW1s3qQgjLM/uvf+LEiTRq1IgrV67g4uLCiRMn2LVrFx07diyTfaGFMJm/v9rkPXs22NnBhg3QqhX8+KPhkrIesGWOpKQkVqxYYdFR1dZoVhdClI7ZCXrPnj3Mnj0bT09PtFotWq2WHj16EBERUeBgGiEKUuZzdu3t4fXX1dHdzZpBQgI8+CBMmAA3b5o84rk8bNiwgfPnzxMVFVVmz2ELzfhCiKKZnaBzcnJwdXUFwNPT0/AB26BBA86cOWOZ6ESlV2412I4d4dAhNTEDLF4M7dvjEx9vczXLJ554ghYtWpTplwZJ0ELYPrMTdKtWrTh69CgAXbp0Yd68eezevZvZs2fj7+9vsQBF5VauNVgXF4iMVJu4fXzgzBnsevbEPTIS8gxyLC8FrVDm7e3N1KlTqV+/vlViEkLYBrMT9GuvvWb49j179mxiYmLo2bMnW7duJTIy0mIBisrNKn2j/frBsWMwbBia7GxqzZ+PXZcu6paWFqbfOrKwmmruFcqsVavVaDS4u7vLFpZC2Biz/yL79+/PkCFDAGjcuDGnT5/m6tWrXLlyhT59+lgsQCHKhIcHrF9PzqefkuPhgeb4cejaVV3c5OZNiz3NpUuXWLhwIZcvXy7w/twrlOl0Os6ePcuiRYssvuZ3YeLj41m0aBE3b96UBC2EjbHoX6SHh4fsBS0qDo0G5cknufzTT+iefBJ0OnjvPezat8d5z558lyclJREZGWnSgLaoqCjOnTvHypUrC7xfv0KZfmOM9evXl3rN7+Jq7Xnjs7VR7EIIlUnzoE3ZBOODDz4wORghrEFXuza6Tz9FO3IkPPccmvPnuWvUKNKeeEJd4MTTE/gveZqy2EdQUBCJiYmMHTu2RNcPHz6cdevW5Vvz2xSmJOigoCDS09OLHAOQe2GZgmrZxd0vhDCPSQn68OHDJbpOatGiQnr4YThxAt1LL6Fdtgy3L75A+fVXNUk/8gjDhw9n/fr1Jg1o8/HxITg4GG9vb65evQoY7zyVN4nWq1eP8ePHY2dnl6+svIkwOzubP/74g3Xr1jFlyhTDQisloZ/eNmHChGLHAORdbzxvHOasRy6EKJ5JCVo2yBCVnpsbukWLSLj3Xuq8+iqOMTEwZAjaBx/krmnT1B20ciWzglYiK44+QQMmDQjLmwizs7P55JNPuHTpkslLeOqnt0VFRZk8gl4SshDlo1R/XSkpKbz//vuEhIQQEhLC/Pnzrb5Eoqg8ctc0y9vtzp259N136MLDwdER7Q8/4DtgALUWLoRbtwzXWXslsoEDB9KsWTNDktUvU1pY4tcn16CgIJo3b05QUFB5hiuEMIHZn3wHDhwgICCA+fPnk5ycTHJyMh988AEBAQEcOnTIkjGKKsqaCVqj0eBWty68+SYcO4aub1+0d+7gsWgRdm3bwpYtgGXmcZvSZ5yXp6cn48ePN9Tq9QlavzZ+3ufR/616eXnZ3AItQghjZn/yTZ48mYEDB3Lx4kU2btzIxo0biYmJ4ZFHHmHSpEkWDFGI8mf05aBpU3Rbt5K4aBHZXl5oLlyARx6BwYPxuXOnRImuqC8biYmJREZGFjoVy1IURSEtLa3A5A2yupgQtqZUNeiXXnoJe/v/urHt7e2ZPn06Bw4csEhwQtgMjYaMAQOI/fFHdFOnqmt8f/01dq1bU/P99yE9vciHF5Wg9Wtv66diWWvhkPJK0PJFQIiSMfsTwM3NjdjY2Hzn4+LiDGt0C1HZKNWro4uIgD//hPvuQ3P7NjUXL8auZUtuLlmCLjvb5DKHDRtGQECAYSqWNZv2cyurRCoJWoiSMfsTYPjw4QQHB7N+/Xri4uKIi4tj3bp1hISEMGLECEvGKCqZMt/Bqjy0aAHR0eR88QVZfn5oEhJwGT8eTadOsHOnSUXVqVOHESNGGBYrsRWSSIWwLrMT9HvvvceQIUMYPXo0DRs2pGHDhowdO5bHH3+cd955x5IxikrG2iOfi1PiLxAaDcpjj3H5p5/IeustdK6uaI4cgd69YcgQOH+e+Ph4li9fbiiroKSnKAq3bt2yWCIsi/2kCyIJXIiyZXaCdnR0ZOHChVy/fp0jR45w5MgRkpOTmT9/Pk5OTpaMUVQytrgHc24mf4FwciJn8mQubd+O7tlnQauFTZugRQtWDR/O+ZMnDXs7501qWq0WV1dXsxb3iY+PZ8WKFfnW7S7JftKWaEaXBC1E2TL5r1On0/HOO+/QvXt3OnXqxOzZs2ncuDGtW7fGxcWlLGIUlYxVdrAygblfIHS1a6NbvBiOHlV3zMrK4tkDB+i5fz9Pp6ZCWlq+x5QmUUZFRXH+/Hk2bNhgOBcfH096ejre3t5Gc5y1Wi3u7u5cuXKFyMhIEhMTbaKfWwhROJP/OufOncsrr7xCjRo1uOuuu1i4cCGhoaFlEZsQhSrLkc6l/gLRqpW65/TPP+PVvj2zsrPxW7QIAgLQLFwImZlGl+t0OjIzM7l8+XK+vaFzy9v0HhQUREBAAMOGDTNcs3z5ci5fvoyLi4tR/PovAhs2bDCsJ65/bnPnYJurUoxBEKIcmPzptnr1aj788EN+/PFHNm/ezLfffsuaNWukmUuUG61Wi4eHBx4eHrZdA7z/fu789htn585F16QJXL2K3dSp1O/TB82qVfDviG99gl69ejVnzpzh66+/LvB15W169/HxYdy4cUaDy/TbVz755JMFljF8+PB8K48VlqD1ifTw4cNFJlRTE66tj0EQwlaY/OkWGxvLQw89ZLjdt29fNBqNfBsWZjGnidca05D0TcQm9xVrNCT36cOdQ4fg449R7roL+/h47EJC1JHgq1cbEvWoUaNo0aIFoaGhBb62kjS962v//v7++crQ6XRUr17daOWxougT6csvv1xkQjV1y0pbH4MghK0w+RMuOzsbZ2dno3MODg5kZWVZLChRddjKnN/ilDpOe3sICSHn1CmSw8NR6tSBc+dgzBhudOzIzx98ADk5RTatl7bpXV9bLmwlsbw1YX0iffvtt4tMqPp1vUuacG19DIIQtsKk3axAnRIyduxYo5Hat2/f5rnnnqN69eqGcxs3brRMhEJUJtWqkfbMM7hNn479Rx/BvHms+Ptvrv/9Nxvuu4/xr78Ozz+vJnQL04/6Hjt2LP7+/vnuz930PGvWLEMi9fLywtvb22i979zN4pJwhSgbJlcHxowZQ926dXF3dzccI0eOxMfHx+icEEKteRe4LWONGjB9Oly8yMjJk2np6Miz169TJywMu5YtYelSo12zLKGgUd+5lbTp+dKlSyxcuLDQtcNNnX4l07WEKJjJX9OLmlspREVT1k3shSZovRo1qD1jBoGtW+N+7hw5y5ZhFxMDL7wAs2bBpElqjbqAL73F1YjzCgoKIj4+3mjUd24lrQlHRUVx7tw5Vq5cydNPP53vflP3i5b9pYUomPw1iCotd4I2pyanf3xiYmKppg5lV6vGrcmTufTbb+QsWAB+fpCUBOHh4OeH8vLLpJ45YxRbcTViPf3r8vLyyjfq2xxBQUE0btzYsHa4EKJsVIgE/cYbb6DRaIyO5s2bWzssUckUlqD154F8tW19gi7JSGadTsfp06dZsGCBUSJPSEjgyy+/JDExEaVaNZTQUHUA2erV0LIlpKWheecdarVvT+3XX4fz54GC50EX97ryzh/PuxRpSfj4+BAcHIy3t3eJHyOEMF2FSNAAd999NwkJCYbjt99+s3ZIooooKkHrlaT/VqfT8fHHH+dL5EuWLGHPnj0sXbr0v4sdHGDUKDh2DL7+GqVLF7R37uC+di12LVrA4MH4nj/PixMm4OXlZRRnYS0ACQkJREVFkZGRYXgNK1eu5OLFi4atLk19T8zpN64oI/eFsLYK8xdib2+Pl5eX4fD09LR2SEIYlLT/dtiwYTRq1MhoGU5QV0YrcPqTVgsDB5Lz669cXrOGm/fei0ZR4OuvcejXj+ZPPYXbV1+RkphIdnY2aWlpuLu7U6tWrXwJcOXKlZw7d85oHElISAitWrUiJCSkgKcuOJHqa9379u0zq1lfErQQJVNh/kL++usvfHx88Pf353//+1+Be1HnlpmZSVpamtEhhLXVrVuXkSNHGmq9AKGhoXTt2pWRI0cWXiPVaLjdpQuXP/mE67t3ozz7LEq1ajidOkW98HBqBQainTULh2vXCl1lbezYsTRu3Njoy0FRXywKS6T6vu933nmHM2fOyMBRIUogbz7KzLPkb0EqRILu0qULK1eu5IcffmDp0qXExMTQs2dPbty4UehjIiIijKZ9+fr6lmPEQqhKUlv09vZmyJAhVK9evdBFRPQURSG5Xj1yFi8m5+JFrk2fTo6PD3bJydi/9Rb1u3dHO2YM/P475CnL29ub4ODgUs9X1vd9v/TSSzRr1ixfa4CpZG1uURX4+voa5aSIiIhiH1MhEvSAAQMYNmwYbdq0oX///mzdupWUlBS++OKLQh8THh5Oamqq4YiLiyvHiIVQlWlzrocHN55/nqzTp7myaBG6rl3RZGWhXbsWuneHtm3V+dQWbj3SrwHesmVLXnzxxWITvr6/Ojs7u8B+a1mbW1QFcXFxRjkpPDy82MdUiASdV82aNWnatCnnzp0r9BonJyfc3NyMDiHKg6kDqLRarWFlvqJ2ltKPwM63HriDAzcfeYQ727cTv3kzujFjwNkZ/vwTXngBOz8/ar/yCtqjR0v1ukoqb424uAQta3OLqiBvPsq9GmdhKmSCTk9P5/z58zLNQ9gkUxO0vb29IfEWtbPUokWLjEZgF+ROYCC6Tz6B+HhYuBBatECTkYHr55/j0qMHrYKDsYuKsvje1Lkfo68Rf/LJJyV6H2SpUCEKViES9NSpU9m5cycXL17k999/57HHHsPOzo4RI0ZYOzRhBhnFa0yr1VKtWjXs7OwKXU1Ln/QiIyPz9dcW2Idbqxa8+CKcOEHsF1/warNmXLa3p8bJkzi88AJ4ecGoUWi2b4d/E2hhv5eS/L5yJ2F9jXjs2LGyhKcQpVAhPiEvXbrEiBEjaNasGU888QS1a9dm79691KlTx9qhCTNUtARtyUFMSUlJrFixosCyiloWVJ/0gHz9tfpFUgpK3mg0LD9+nIMNGzI/JISLL7yArlkzdZ3vzz7Drl8/6vfsifaNNwwLoBQUlym/L6kRC2EZFeITct26dcTHx5OZmcmlS5dYt24dAQEB1g5LVDL6RJyYmGiUkCw1iCk+Pp6pU6dy8uRJk6cm6ZPeiy++mK+/Vr/dI8CpU6dYsmRJvhpts2bNGPTUU+imTEF37Bjs3QvPPYfi7o59fDzauXOhcWPo1QuioqCIGRL6LxKm7o0to7WFME2FSNBClAd9Io6KijJK0JYaxLR8+XKcnZ3JzMw0e2pSQbXT3Mm7WbNmDBo0yChB+/j4GEZb+/n5Ye/gAF26oFuyhOsnTpD+8cco/fqBRgO7dsG4cWoT+FNPwTffQJ75msVuAJJL7v74kiyHKoT4jyRoIf5VWCI2tcm2sJpicHAwgYGBzJw506gsfcL7559/iIyMNLuGqU/EJd0MQ6fTkZaVhfPYsWh+/BFiY+Gtt6BpU7h5Ez7/HAYNAi8vlOBg0jdvRpeVBahN9YsXLyYpKanI51AUhb/++ouFCxfy0EMPFfj+VrQuDyHKi/xFCPEvS/WdFtYk7uPjw/jx4/MtU6vVasnIyGDatGkcPXrUejXM+vUhPBzdyZPc2LYNZeJE8PaGlBQ0K1ZQ47HH0Pj6op00ia8WLODM6dOsX7++0C8k8fHxREZGEhUVxZkzZ9i6dWuB768kaCEKJn8RQuRRmm0ntVptkU3iGo0mX/OwVqvl66+/plq1aty+fbvYDTeKmittCTpF4VrjxuS89x7ExcH27ehCQsipWRNNUhLaJUuYsmEDXfbu5bmkJKJmzeL0qVP5vljoEzNgkRXHLKE0m3wIUd7srR2AELZG/yFe0n5W+C9BQ/HrWxeUoENDQ7l9+zbDhw8vsgafO0GXtsapr/lOnjwZPz+/gi+ys4PevdH16MGl6dPxPX0azfr1eG3ezOwbN2D9ep4H/s/Fhac9PWHbNujRA1AHr6WmpvLwww/j7++Pj48PiYmJpYq5tMz53QphLfI/VAgbUNL+48ISvDlNxCaPTndwQBkwAN3KlZz//XdOzp3Lrccfx9vdnVk3b1J/zRro1w87Hx88p0yh/sGDPPHoo6xfv97qiVmIikgStBAVgD4J29vbGxJ07nPFJeikpKR8A9AKa4rX16wvXbpUaHOw4uzM9V690Hz6KVmXLnFx2TJygoOhbl00KSnU2LgRxyee4OcBA/jn66/55vXX0V66hJubG2lpaRWiiVmaw4W1SYIWohzp19M2tbZbUC25oHMFJWKA9evXc+bMGZYvX25IPF5eXgU2xeeuWReUoOLj41m8eDFXr14FQGdvT2JgIFmLF6tLjO7aBZMmofj5MTAnh/TkZEbs3IljkybU6t0b7euvo9u9G3JyTHoP8irrBCoJWlibJGghylFpRywXtdiHfrDZ6dOnjRYr0Wq1PP3004backGJJ/c5fc26sEFdy5cv5+jRo8yaNSt/HHZ20LMnzJ9PzrlzfPrUUzg1bMhGb28UrRbNsWPU/PBD7O+9V51rPWYMfPEFXLtm8nthTgKVxVJERSIJWohi2NI0oKL6jfWDzZo3b260WIlWq6V58+ZMmjSp0AFouZNdcdPNgoODuX37Nk5OToYV0QpaVUxrZ8eYadPwue8++q9fT87ly+SsXEn6o4+iuLvD1auwejUMHw516kCnTvDKK7B9e77FUSyluH53qTULW2L9TxwhbJypCdrchF6S5u/iVjUrarBZcTXngq7XjxjPXfP08fHhvffeo2XLloZatqIo+R6blpZG06ZNCQ4OVuPx9OTyfffxekAAcQcPwo4dMHUqtGoFigIHDkBEBPTpAx4e8NBDMH8+HD+u3m8Bxb1/ud8PqW0La5MELYSFmbsrVEkSe2kWU1EUhdTUVC5dumRYc7yoQVu5E3Temme9evUYN24cPj4+xMfHs3z5cqNEljvR5f7iYShn9Wp13e9334Vjx+DyZVi1CkaOhHr11JXMvv8ewsKgdWu46y61OXz1avj7b5Nfuznvn6XWYBfCXJKghcijrGpOttJUrl8TO++a4/BfzdfNzc3ofO6ap1arxc3NzdCsHRUVxblz51ixYgUpKSlkZ2cbkrP+2lq1ahW6iItOpyPFxQXdyJHw6aeQkABHj8J770H//lCtmnpu9Wo1STdsCA0bog0KosaGDXDhgsVq2LlZag12IcwlCVqIf+kT6MqVK8u95lSeyVu/+1VBiaewBJ275pk31qCgIBo3bsyYMWMMCTo1NRVFUbC3t6d+/frY29vnK0dfy86d0EFdySzFzw/d5Mnotm4l5cIFdD/9BC+9BF27qgPR/v6bxM8+Y8n06Vxp1gz8/NTa98cfw9mzZifsvE35trJtpvSNV02SoIX4lz7xhISElHvNyZIJurCpVnoFJZ6CWg0KarrOHa9+PraPjw/BwcF4e3ub9LoK6xNPTk4mOTkZnU6nXnP7Nrr77oO334Y9eyAlBX78kY969OBwzZp8pNXCpUuwZg088ww0a6Y2kw8eDPPmwW+/qftfFxMP2G6ztiToqkkStBB52FLNKa+SJPLcc56Lo0/MkZGRnD59msjISEOi1jddF7R3dd6acd4y33rrLSIiIoiPjzcpueh0OkPtu1A1akC/foz99FPueuopxh4/Dj//DK+/DvfeC46O8M8/8PXXaq27Z09wd4cuXdBOmULNbdvQJiQUWLQ0awtbIglaiAqkJAl6+PDhNGvWLF+SKahmra8xAjRv3hzAUIN86KGHSEtL46GHHioyDn1NO+HfpBcVFcXOnTvZsWOH0cIo+j7pvPHnrb3njrOo8QBeXl6MGjUKr0aN4P77YfZs2LkT0tLg99/VPuwhQ9T51llZsH8/LFgATzyh7tzl56dO8XrvPXVE+Y0bRTblW5OMKK+aZLMMISoRrVZL06ZN87UA6BcxOXPmDCtXriQsLMwwaCstLY0JEyag1WrR6XQsWrSI4OBgPv74Y9zc3Ni6dSudOnUq9DlXrlzJxYsXWbx4MVqtlieeeIL4+Hg0Gk2JaqL6QWvLly/n9ddfz9cCoL9v1qxZJXsTnJygWzf1mDJF7Y++eFFtHv/9d/U4elTdqSsuTl0oBUCjwa5ZMzxbtEDTsyfarl2p2bYt2ECCzt30XuL3QVR4kqCFqET0tb709PR850NDQ5k/fz4hISH5dt6qW7cuFy5cwN/f35Dcg4KCSExMLHabyJCQEG7cuIGiKJw+fZqtW7fyyiuvGMrPTV+bdnFxMUzvCgoKIj09neDgYOLj48nIyMDHx4dhw4ZRo0YN0tLSStfkrNFAo0bq8dRT6rn0dNi3T61VHzgAf/wBcXFoTp+mxunTsGmTep2dnTpPu1Mn6NAB2rZVp31Vr25+PGbQf5GqiE3vhQ08FMWTBC1EFVFU33pB+0zrB38V1xevL1df+w4KCiq6D/nf5zt79izr169n4sSJhrjmzJnDpUuXaNq0KdWqVSt0vXD4bxBbeHg4/v7+JXwX/lWjhtosfv/9/51LSkK3fz+Zv/+O87FjaA4cgKQktbZ99Ch88ol6nUYDTZpAYKCasPX/+vio95UBWx4XURzZ4tN88m4JISyiJEkkd1+qvik79yC04OBgmjVrxrBhw4p9vqIGsRU2hatI9eqhffRRqkVEoPnuO3XudVwcbNwI4eHw4IPg7a02mZ89Cxs2wKuvwiOPqH3adepA375qs/rq1XDokLrgShUn/efmkxq0EJVMaQY3mbsKWt5r3d3dSUtLA/77gJ48ebKhLzUyMpL09HTq169PcHCwoWx9kndxceFmMcktODiY5OTkQudzp6Sk4OzsXGztrdAmWI1GTbz168Njj/13/soVtUZ95Ih6HD0Kp0+rG35ER6tH7jIaNYK771aPli3Vf5s3BxeXYt/LykD6z80nCVqIPGxp9K459PGbcm1SUlKRzcWmlunh4YG9vb3x8p7Llxv6UhVF4fLlyzRr1sxoupZ+pPfNmzdZtGgREydOLPR3Ub9+faZNm0b9+vVLFFdhTG6CrVsXHnhAPfRu34YTJ/5L2H/+qd6+elVd6ezCBfj22/+uLyhxt2gBTZuCm1u+p6zI/ycrcv+5tUmCFiIPU5JRRad/rR988AEXL15k5cqVzJ4923Cfuf2Gud/D3B/QufurFy5cyLBhw0hJScHDw8PoeXIvR1qaWld8fDyLFi1i8uTJ+Pn5mV0OFDPYydlZHUTWoYPx+StX1ER98qT6r/64dq3gxA3qIivNmqnJ+t9D27QpNf39bWJEuakqcv+5tUmCFqKKKKoWph+JHRISYjhX1GIkJS0XCv6A1u+6lZOTQ2pqquHx+rKCg4MNI7vNeV59s7qiKIYpW4Ul+txN8EUlcbMGO9Wtqx733fffOUVRF1LRJ2t98j5zRh2Upj927cr7gtV1yJs2/S+BN2kC/v7qnG4Hh5LFVM4qcu3f2jRKccMtK4m0tDTc3d1JTU3FrYAmJCGqsuzsbC5dulSihFyYwmqYucsGuHTpEl5eXly+fJmcnBzs7Oxo0KCB0fOWNp6ZM2eyf/9+mjVrhkajKTL56q/t3LlzkbV1S7xHxUpNhb/+Ugeh5T7OnFGnhhVGqwVfXzVZ5z0aNQJPzxKNMJcpUZZXmtwjNWghhEWUpGtAf429vT21atXi0qVLrFixgrCwsFI3QeeWdwGWoppXS9pHWtKadqm4u0PHjuqRm6JAYmL+pH3+vNpMfvu2ug3n33/D9u35y61RwzhpN2yoJnQ/P/WoXRs0GpkSZWOkBi2EKNPaYVFlF1Z7LW08+sd7eXmRmJhYZDklfa6S1rTLnT55x8T816+d+7h8ufgyqlUDX190vr5k1K6NS/Pm2DVo8F8C9/VVrxEmkxq0EMJmFdUHWdYjfC3Z/2mzo5E1GnV+trc33HNP/vv1tevcSTs29r8jMVHd7evsWRLPnuVj4GkgX5uDp6earO+6S12UJe/h7a3OBZeat8VIDVoIYbW+x8Jqr6WNx5THl7QGXS590NaQmalu1xkXx8wFCzh04gQdXF15o149NYH//TdkZJSsLHt7dXOSghK4PonXq6cmezu7sn1dNkJq0EKIUrG1qWWljcecueDFJfJKOxrZyQkCAiAggGB/f9Lmz2fc5MlqbRnUJvSUFHVVtdhYiI8v+LhyBbKz1WR/6VLRz6nRqEm6bl01YRf3bxVtXpcatBDCaiptrbSCKtXvIytLnR6WkFB4Eo+PV+eAm5p2atT4L2HXrasOaqtdW03y+p9z3/bwUGvzNkBq0EKICqnS1kqrIgeH/5ZGLUp2tpqkk5LUWrf+39w/5/43M1OdYpaero5aLyl394KTuP6chwfUqgU1a6pHrVrqYxwdS/MuWJQkaCGE1dha03pVVy5fmOzt1dpwvXrFX6socONG/kR+7Zq6jOq1a/8d+tvXr6uPTU1VjwsXTIvPxcU4aet/Lu62u7u6TKsFF4yRBC2EEAKwwS9MGo2a9NzcoHHjkj0mJ0dN0kUlcf2Rmqpem5IC/27uws2b6mHu7lvVqv2XrN3dS7UpiiRoIYQQlYedndqE7empLolaUjk5apLWJ+zcR3Hnrl//b2vRW7fUIzGx1C9FErQQQghhZ6c2WdeqZd7js7PVBJ+WptbM9f8mJMAzz5hVpCRoIYQQorTs7dWBZx4exufT0sxO0BVq6OSSJUto2LAhzs7OdOnShf3791s7JCGEEKJMVJgEvX79esLCwpg5cyaHDh0iMDCQ/v37c+XKFWuHJoQQQlhchVmopEuXLnTq1InFixcD6lJ+vr6+TJgwgZdffjnf9ZmZmWRmZhpup6am4ufnR1xcnCxUIoQQolykpaXh6+tLbGws7u7uhvNOTk44OTkV/WClAsjMzFTs7OyUTZs2GZ0fPXq0MnDgwAIfM3PmTAWQQw455JBDDps7Zs6cWWzuqxCDxK5evUpOTg718kxsr1evHqdPny7wMeHh4YSFhRlup6Sk0KBBg3zfYkTpdOrUiT/++MPaYZSKLb2G8oylLJ/LUmVbohxzyzD1cfqakrTSWZYt/X2aQ1EU2rdvz8GDB40WgCm29kwlHsVdWPOBu7u7/PFYkJ2dXYV/P23pNZRnLGX5XJYq2xLlmFuGuY9zc3Ozmf9PlYEt/X2ay8nJyawFYCrEIDFPT0/s7OxISkoyOp+UlISXl5eVohIAoaGh1g6h1GzpNZRnLGX5XJYq2xLlmFuGLf2/qMoqw+/B3NdQoQaJde7cmUWLFgHqIDE/Pz/Gjx9f4CCxvGQ3KyFEWZLPGGFpFaaJOywsjDFjxtCxY0c6d+7MggULyMjIICgoqESPd3JyYubMmSVq9xdCCFPJZ4ywtApTgwZYvHgx7777LomJibRt25bIyEi6dOli7bCEEEIIi6tQCVoIIYSoKirEIDEhhBCiqpEELYQQQtggSdBCCCGEDZIELYQQQtigKpmg4+Li6N27Ny1btqRNmzZs2LDBcN9jjz1GrVq1ePzxx60YoRCiIivsM6aozx4h8qqSo7gTEhJISkqibdu2JCYm0qFDB86ePUv16tXZsWMHN27cYNWqVXz55ZfWDlUIUQEV9hmTlpZW6GePEHlVyRq0t7c3bdu2BcDLywtPT0+Sk5MB6N27N66urlaMTghR0RX2GVPUZ48QeVXIBL1r1y4effRRfHx80Gg0bN68Od81S5YsoWHDhjg7O9OlSxf2799fYFkHDx4kJycHX1/fMo5aCFFRlMdnjHz2iOJUyASdkZFBYGAgS5YsKfD+9evXExYWxsyZMzl06BCBgYH079+fK1euGF2XnJzM6NGj+eijj8ojbCFEBVHWnzHy2SNKpNgdo20coGzatMnoXOfOnZXQ0FDD7ZycHMXHx0eJiIgwnLt9+7bSs2dPZfXq1fnK3L59uzJ06NAyi1kIUXFY+jOmqM8eIXKrkDXooty5c4eDBw/St29fwzmtVkvfvn3Zs2cPoG6gPXbsWPr06cOoUaOsFaoQogIqzWeMfPYIU1S6BH316lVycnKoV6+e0fl69eqRmJgIwO7du1m/fj2bN2+mbdu2tG3blmPHjgHQt29fhg0bxtatW6lfv77hD04IIaB0nzFFffYIkVeF2W7Sknr06IFOpyvwvp9//rmcoxFCVDZFfcYUdl6IvCpdDdrT0xM7OzuSkpKMziclJeHl5WWlqIQQlYV8xojyUukStKOjIx06dCA6OtpwTqfTER0dTbdu3awYmRCiMpDPGFFeKmQTd3p6OufOnTPcjomJ4ciRI3h4eODn50dYWBhjxoyhY8eOdO7cmQULFpCRkUFQUJAVoxZCVBTyGSNsgpVHkZtl+/btCpDvGDNmjOGaRYsWKX5+foqjo6PSuXNnZe/evdYLWAhRochnjLAFVXItbiGEEMLWVbo+aCGEEKIykAQthBBC2CBJ0EIIIYQNkgQthBBC2CBJ0EIIIYQNkgQthBBC2CBJ0EIIIYQNkgQthBBC2CBJ0EIIIYQNkgQthBBC2CBJ0EIIIYQNkgQthBDlKC4ujt69e9OyZUvatGnDhg0brB2SsFGyWYYQQpSjhIQEkpKSaNu2LYmJiXTo0IGzZ89SvXp1a4cmbIzUoEWV1rt3byZNmmTtMGxeZXmfrl27Rt26dbl48aLVYvD29qZt27YAeHl54enpSXJysuH+J598kvfff99K0QlbIglaVFiJiYlMmDABf39/nJyc8PX15dFHHyU6OtraoVnE2LFj0Wg0aDQaHBwcqFevHg888AArVqxAp9OVaywbN25kzpw5htuWTNi9evVCo9Hw+eefG51ftGgRPj4+FnkOvblz5zJo0CAaNmxo0XLNdfDgQXJycvD19TWce+2115g7dy6pqalWjEzYAknQokK6ePEiHTp04JdffuHdd9/l2LFj/PDDD9x3332EhoZaOzyLefDBB0lISODixYt8//333HfffUycOJFHHnmE7OzscovDw8MDV1dXi5erKAqHDx/G29ubr776yui+gwcP0r59e4s9182bN1m+fDnBwcEWK7Mwbdu2pVWrVvmO+Ph4wzXJycmMHj2ajz76yOixrVq1IiAggM8++6zM4xQ2ThGiAhowYIBy1113Kenp6fnuu379uqIoinL79m1lwoQJSp06dRQnJyele/fuyv79+42u7dWrlzJx4kTD7QYNGijz5883uiYwMFCZOXOm4frx48crEydOVGrWrKnUrVtX+eijj5T09HRl7NixSo0aNZSAgABl69at+Z5nwoQJyrRp05RatWop9erVM5RZmDFjxiiDBg3Kdz46OloBlI8//thwLicnR3nrrbeUhg0bKs7OzkqbNm2UDRs2mPT8GzZsUFq1aqU4OzsrHh4eyv333294f3O/T2PGjFEAoyMmJkZZtWqV4uHhody+fduo3EGDBikjR44s8DWeOXNGAZQlS5YoLi4uSkZGhuG+1q1bKzNmzCjyPTLFhg0blDp16uQ7n5OTo7zzzjtKQECA4ujoqPj6+ipvvvmmoijm/76Lc/v2baVnz57K6tWrC7x/1qxZSo8ePUx/kaJSkRq0qHCSk5P54YcfCA0NLXBgTc2aNQGYPn06X331FatWreLQoUM0btyY/v37G/X3mWPVqlV4enqyf/9+JkyYwPPPP8+wYcO45557OHToEP369WPUqFHcvHkz3+OqV6/Ovn37mDdvHrNnz2bbtm0mP3+fPn0IDAxk48aNhnMRERGsXr2aZcuWceLECSZPnszIkSPZuXNniZ4/ISGBESNGMG7cOE6dOsWOHTsYMmQISgFjSBcuXEi3bt14+umnSUhIICEhAV9fX4YNG0ZOTg7ffPON4dorV66wZcsWxo0bV+BrOXjwIM7OzoSEhODm5sb3338PwO3btzl16pRFa9C//vorHTp0yHc+PDyct99+m9dff52TJ0+ydu1a6tWrZ7jf3N93YRRFYezYsfTp04dRo0YVeE3nzp3Zv38/mZmZ5r1YUTlY+xuCEKbat2+fAigbN24s9Jr09HTFwcFBWbNmjeHcnTt3FB8fH2XevHmGc+bUoHPXbLKzs5Xq1asro0aNMpxLSEhQAGXPnj1Gz5O3RtSpUyflpZdeKvQ1FFaDVhRFGT58uNKiRQtFUdTamIuLi/L7778bXRMcHKyMGDGiRM9/8OBBBVAuXrxY4PPlfZ/y3tZ7/vnnlQEDBhhuv//++4q/v7+i0+kKLHfq1KlK586dDY998sknFUVRlL179yqAEhsbW+DjzDFo0CBl3LhxRufS0tIUJycno9aI3Mz9fRfl119/VTQajRIYGGg4/vzzT6Nrjh49WuTvQ1QN9tb8ciCEOZQSzAw8f/48WVlZdO/e3XDOwcGBzp07c+rUqVI9f5s2bQw/29nZUbt2bVq3bm04p699XblypdDHgTqaN+81JaUoChqNBoBz585x8+ZNHnjgAaNr7ty5Q7t27Ur0/IGBgdx///20bt2a/v37069fPx5//HFq1aplUlxPP/00nTp14vLly9x1112sXLnSMNitIIcOHTLUkocMGcKQIUPIzMzk0KFD1KlTx2jw1AcffEBycjJvvvkmANnZ2dSuXZtr165hb1/8R9mtW7dwdnY2Onfq1CkyMzO5//77C32cub/vwvTo0aPYQX7VqlUDKHGtXFROkqBFhdOkSRM0Gg2nT5+2eNlarTbfF4CsrCyj2w4ODka39aOsc98G8n0IF/Q4c0djnzp1ikaNGgGQnp4OwJYtW7jrrruMrnNycirR89vZ2bFt2zZ+//13fvrpJxYtWsSrr77Kvn37DM9TEu3atSMwMJDVq1fTr18/Tpw4wZYtWwq9/tChQ4wYMQJQR4Y7ODjw448/FjhA7Pjx4zz00ENG70HDhg1LlJwBPD09uX79utE5fSIsirm/79LQd8PUqVPHYmWKikf6oEWF4+HhQf/+/VmyZAkZGRn57k9JSSEgIABHR0d2795tOJ+VlcUff/xBy5YtCy27Tp06JCQkGG6npaURExNj2RdQSr/88gvHjh1j6NChALRs2RInJydiY2Np3Lix0ZG7BlocjUZD9+7dmTVrFocPH8bR0ZFNmzYVeK2joyM5OTkF3hcSEsLKlSuJioqib9++hcZw4cIFUlJSDInY3t6egQMH8tVXX3Ho0KF8/cXHjx8nMDDQcPvIkSNGNdnitGvXjpMnTxqda9KkCdWqVbO5qXnHjx+nfv36eHp6WjsUYUVSgxYV0pIlS+jevTudO3dm9uzZtGnThuzsbLZt28bSpUs5deoUzz//PNOmTcPDwwM/Pz/mzZvHzZs3i5xm06dPH1auXMmjjz5KzZo1mTFjBnZ2duX4yoxlZmaSmJhITk4OSUlJ/PDDD0RERPDII48wevRoAFxdXZk6dSqTJ09Gp9PRo0cPUlNT2b17N25ubowZM6bY59m3bx/R0dH069ePunXrsm/fPv755x9atGhR4PUNGzZk3759XLx4kRo1auDh4YFWq37ff+qpp5g6dSoff/wxq1evLvQ5Dx48iKOjI61atTKcGzp0qGHA1auvvmo4rygKp0+fNmrGv379OuHh4cW+Nr3+/fsTHh7O9evXDU33zs7OvPTSS0yfPh1HR0e6d+/OP//8w4kTJ8plOlZhfv31V/r162e15xe2QRK0qJD8/f05dOgQc+fOZcqUKSQkJFCnTh06dOjA0qVLAXj77bfR6XSMGjWKGzdu0LFjR3788cci+1XDw8OJiYnhkUcewd3dnTlz5li1Bv3DDz/g7e2Nvb09tWrVIjAwkMjISMaMGWNIiABz5syhTp06REREcOHCBWrWrEn79u155ZVXSvQ8bm5u7Nq1iwULFpCWlkaDBg14//33GTBgQIHXT506lTFjxtCyZUtu3bpFTEyMYfEPd3d3hg4dypYtWxg8eHChz3no0CFatWqFo6Oj4dwDDzxATk4Od+7cMWrijomJoWnTphw4cMBwbsCAAbRu3ZoNGzawbds2rl+/zsyZM40Sfm6tW7emffv2fPHFFzz77LOG86+//jr29vbMmDGD+Ph4vL29ee6550rytpWJ27dvs3nzZn744QerxSBsg6zFLYSwuPvvv5+7776byMhIi5T3zTffsGHDBj799FPDuQYNGrBr1y4aNGgAqE3eP/30E9OnTy+0nC1btjBt2jSOHz9u9AXHlixdupRNmzbx008/WTsUYWVSgxZCWMz169fZsWMHO3bs4MMPP7RYucePHzcaO5CRkUFqaqohOet0Oj788ENmzJhRZDkPP/wwf/31F5cvXzapf748OTg4sGjRImuHIWyAJGghhMW0a9eO69ev884779CsWTOLlXv8+HGeeOIJw+3Tp08bErZOp2Py5Mk8//zz1K9fv9iybH3Tj5CQEGuHIGyENHELISq0efPmsXnzZtq0acMDDzxgGN0uREUnCVoIIYSwQbY5SkIIIYSo4iRBCyGEEDZIErQQQghhgyRBCyGEEDZIErQQQghhgyRBCyGEEDZIErQQQghhgyRBCyGEEDZIErQQQghhgyRBCyGEEDZIErQQQghhgyRBCyGEEDbo/wF2CMIGjPKfqQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Polarization fraction as a function of the column density\n",
    "PvN_plot_Step2 = HtoH.PvN(CatP, CatdP, CatN, Nscale=[1e+21,3e+22], Pscale=[0.0,30.0], showfit='true', weighted='false', errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "c49996f1",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAegAAAF2CAYAAAClPAs/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABJBklEQVR4nO3deVwTZ/4H8E8ilyiJAiJQENC6HkURqUVrFavWo6fW3bWuVqja7bZ4UrXQulq1bVBbrbqubg+Ptlpbu2pva0sVeyhWBBUPFIVCK4qKgiByZX5/9Jcs4Uwmk2SSfN6vV16ayWTynSdhvvPM88zzKARBEEBERESyorR1AERERNQQEzQREZEMMUETERHJEBM0ERGRDDFBExERyRATNBERkQwxQRMREckQEzQREZEMMUETERHJEBO0jFRWVuLll19GZWWlrUOxeyxL6bAspcOylJajl6dCDkN9rl+/HuvXr0deXh4A4K677sLChQsxevRoAMDt27fx/PPPY/v27aisrMTIkSPx73//Gx07drRh1NIrLS2FWq1GSUkJVCqVrcOxayxL6bAspcOylJajl6csatBBQUFITk5Geno6jhw5gqFDh+Kxxx7DyZMnAQBz5szB559/jh07diA1NRUXL17E448/buOoiYiILMfF1gEAwCOPPGLw/NVXX8X69etx6NAhBAUF4d1338W2bdswdOhQAMCmTZvQo0cPHDp0CP3797dFyERERBYliwRdV21tLXbs2IHy8nIMGDAA6enpqK6uxvDhw/XrdO/eHZ06dcLBgwebTNCVlZUG7RJarRbFxcXw8fGBQqGw+H6IUVpaavAviceylA7LUjosS2nZQ3kKgoBr167B29sbSuX/Llq7u7vD3d29xTfLwvHjx4U2bdoIrVq1EtRqtfDll18KgiAIW7duFdzc3Bqs369fP2H+/PlNbm/RokUCAD744IMPPviQ3WPRokUt5kXZ1KC7deuGzMxMlJSU4JNPPkFsbCxSU1NFby8pKQkJCQn65yUlJejUqRMKCgocsjMBERHJT2lpKYKDg5Gfnw+1Wq1f3mLtGTK6xO3m5oY777wTABAVFYVffvkFq1evxvjx41FVVYUbN26gXbt2+vUvX74Mf3//JrfX1OUDlUrFBE1ERFalVqtNzj2y6MXdGK1Wi8rKSkRFRcHV1RUpKSn617Kzs5Gfn48BAwbYMEIiIiLLkUUNOikpCaNHj0anTp1w8+ZNbNu2Dfv378c333wDtVqNqVOnIiEhAd7e3lCpVJgxYwYGDBjAHtxEROSwZJGgi4qKMHnyZBQWFkKtVqN379745ptv8MADDwAAVq1aBaVSiXHjxhkMVEJEROSoZDGSmDU4+ogzREQkP+bkHtm2QRMRETkzJmgiIiIZYoImIiKSISZoIiIiGWKCJiIikiEmaCIiIhligiYiIpIhJmgiIiIZYoImIiKSISZoIiIiGWKCJiIikiEmaCIiIhligiYiIpIhJmgiIiIZYoImIiKSISZoIiIiGWKCJiIikiEmaCIiJ7F582YoFAooFAr9siFDhkChUCAuLs52gVGjmKCJiByALtE29ti9ezcAoEOHDoiOjkZ0dLRtg5VQXFwcFAoFhgwZYtT6M2bMQEREBFxcXKBQKODv72/ZAM3gYusAiIhIOm5uboiMjDRY5u3tDQB46KGH8NBDD9kiLNl4//334ebmBm9vb1y5csXW4TSLNWgiIgcSEBCAQ4cOGTwGDx4MoPFL3I2JjY1F165d4eXlBTc3N4SEhGDmzJkoLS3Vr6OruYaGhmLLli0IDg6GSqXC7NmzUVFRgdmzZ0OtVqNTp07YsGGDwfYvXryIKVOmIDAwEG5ubujcuTOWLl2Kmpoa/Tq6KwKTJ0/GokWLEBAQgPbt22PSpEm4efMmAOg/GwBSU1P1+7Z///4m9+3EiRMoKirCgw8+aFK52gJr0EREZODTTz+Fi4sLunTpgps3b+LChQtYu3YtCgsLsWPHDoN1CwsL8dxzzyEgIAA3b97E6tWrsXfvXvz+++9Qq9UoKChAfHw8YmJi0KNHD1y7dg39+/dHQUEBvLy80KNHD5w6dQoLFy5Ebm4uNm7caLD97du3w8PDA76+vrh06RK2bt2KkJAQvPrqq4iMjER5eTmuXr0KLy8v9OzZEwCgUqma3Lfg4GDpC8xCWIMmInIgv/76a4M2aFOlpqbi6tWryMzMxPnz5/HSSy8BAHbv3o3bt28brFtVVYW9e/fi7Nmz+uSXk5ODY8eO4dSpU/Dw8IBWq0VqaioA4F//+hcKCgrQsWNHnD9/HseOHcMnn3wC4I8afk5OjsH2PTw8cPr0aeTk5CAqKgoAkJKSAgDYtWuX/pJ937599VcM+vbta/I+yxFr0EREDqSxNmhTfffdd5g4cSLOnz9vkJBrampw5coVg1po+/btMXDgQABAp06dUFBQgPDwcISGhgIA/Pz8kJ+fj8uXLwMADh8+DAC4fPky/Pz8DD5XEASkpaXhzjvv1C8bOnQo7rjjDgBA9+7dkZ6ert+Wo2OCJiJyILo2aLG2bt2KuXPn6rcVHByMq1ev4sKFCwCA2tpag/XrXk52cXFpsExXgxcEweB9dS9J1+Xp6WnwvF27dg22X39bjoqXuImISE+X3L28vJCbm4u0tDSMGDFCsu3369cPwB/Jdvv27frL0t9++y2ee+45jB071qTt6RJ6eXm5ZDHKBWvQRESk17t3bwDAzZs30blzZ7i7u6OkpESy7cfHx+Odd97B77//jm7duqFHjx64efMmCgoKUF1djcmTJ5u0ve7duwMAjhw5gl69eqFNmzbYt28fWrdu3ej6Q4YMwW+//YaioiIAwNWrV/WX1Ldu3Sqre8RZgyYiIr2pU6ciISEBvr6+uHnzJoYMGYIlS5ZItv0OHTrg0KFDeOqpp+Dj44OTJ0+ioqICgwYNwqpVq0ze3pQpUzBu3Dio1WpkZWUhLS2twWX4uvLy8nD+/Hn9rVq1tbU4f/48zp8/j4qKCtH7ZQkKwUku5peWlkKtVqOkpKTZLvhERERSMSf3sAZNREQkQ0zQREREMsQETUREJENM0ERERDLEBE1ERCRDTNBEREQyxARNREQkQ0zQREREMsQETUREJENM0ERERDIkiwSt0WjQr18/eHl5wc/PD2PGjEF2drbBOkOGDGkwCfk//vEPG0VMRERkWbJI0KmpqYiPj9dPOVZdXY0RI0Y0mD7s6aefRmFhof6xfPlyG0VMRERkWbKYbnLPnj0Gzzdv3gw/Pz+kp6dj8ODB+uWenp7w9/e3dnhERERWJ4sadH26uUe9vb0Nlm/duhW+vr4IDw9HUlISbt261eQ2KisrUVpaavAgIiKyhfr5qLKyssX3yC5Ba7VazJ49GwMHDkR4eLh++d/+9jd88MEH2LdvH5KSkvD+++9j0qRJTW5Ho9FArVbrH8HBwdYIn4iIqIHg4GCDnKTRaFp8j+zmg3722Wfx9ddf48cff0RQUFCT633//fcYNmwYcnJy0KVLlwavV1ZWGpyhlJaWIjg4mPNBExGR1ejmgy4oKDDIPe7u7nB3d2/2vbJog9aZPn06vvjiCxw4cKDZ5AwA0dHRANBkgjZm54mIiKxBpVKZXDmURYIWBAEzZszArl27sH//foSFhbX4nszMTABAQECAhaMjIiKyPlkk6Pj4eGzbtg2ffvopvLy8cOnSJQCAWq1G69atcf78eWzbtg0PPvggfHx8cPz4ccyZMweDBw9G7969bRw9ERGR9GTRBq1QKBpdvmnTJsTFxaGgoACTJk1CVlYWysvLERwcjLFjx2LBggVGXzLQtQOwDZqIiKzFnNwjixp0S+cIwcHBSE1NtVI0REREtie726yIiIiICZqIiEiWmKCJiIhkiAmaiIhIhszqJFZdXY1Lly7h1q1b6NChQ4Oxs4mIiEgck2vQN2/exPr16xETEwOVSoXQ0FD06NEDHTp0QEhICJ5++mn88ssvloiViIjIaZiUoFeuXInQ0FBs2rQJw4cPx+7du5GZmYmzZ8/i4MGDWLRoEWpqajBixAiMGjUK586ds1TcREREDs2kgUomTJiABQsW4K677mp2vcrKSmzatAlubm6YMmWK2UFKgQOVEBGRtZmTe2Qxkpg1MEETEZG12XwksatXryItLQ21tbXo168fJ7AgIiIyk9kJ+r///S+mTp2KP/3pT6iurkZ2djbWrVuHp556Sor4iIiInJLJvbjLysoMni9evBiHDx/G4cOHkZGRgR07duCll16SLEAiIiJnZHKCjoqKwqeffqp/7uLigqKiIv3zy5cvw83NTZroiIiInJTJncTy8vIQHx8PNzc3rFu3DufPn8cTTzyB2tpa1NTUQKlUYvPmzXjwwQctFbMo7CRGRETWZtVOYqGhofjyyy/x4YcfIiYmBjNnzkROTg5ycnJQW1uL7t27w8PDw9TNEhERUR2ix+KeMGECfvnlFxw7dgxDhgyBVqtFnz59mJyJiIgkIKoX91dffYXTp08jIiIC77zzDlJTUzFx4kSMHj0aS5YsQevWraWOk4iIyKmYXIN+/vnn8dRTT+GXX37BM888g6VLlyImJgZHjx6Fh4cHIiMj8fXXX1siViIiIqdhcicxHx8f7N27F1FRUSguLkb//v1x9uxZ/eunTp3CM888gx9++EHyYM3BTmJERGRt5uQek2vQbdq0QW5uLgCgoKCgQZtzz549ZZeciYiI7I3JCVqj0WDy5MkIDAxETEwMli5daom4iIiInJqoyTKuXbuGCxcuoGvXrmjXrp0FwpIeL3ETEZG1We0Sd35+PoA/2qH79evXYnL+/fffTQqGiIiI/mBSgu7Xrx+eeeYZ/PLLL02uU1JSgrfffhvh4eH473//a3aAREREzsik+6BPnTqFV199FQ888AA8PDwQFRWFwMBAeHh44Pr16zh16hROnjyJvn37Yvny5bIb7pOIiMheiGqDrqiowJdffokff/wRv/76KyoqKuDr64vIyEiMHDkS4eHhlojVLGyDJiIiazMn94hK0PaICZqIiKzNqvdBExERkeUxQRMREckQEzQREZEMMUETERHJkKjpJn/99VccP34cHTt2xD333CN1TERERE7P5AT94YcfIi4uDtXV1VAoFPrpJTt06GCJ+IiIiJySyZe4Fy9ejL/97W84c+YM9u7dCwBITEyUPDAiIiJnZvJ90G5ubjh79ixCQ0MBAGfOnEFUVBTKy8stEZ9keB80ERFZm1Xvg66pqYGnp6f+effu3aHVanHp0iVTN0VERERNENWLe8uWLfj5559RVlYGAHBxccGtW7ckDYyIiMiZmdxJbNCgQXjllVdw8+ZNKJVKhIWF4fbt23j33XcxfPhw3H333fDy8rJErERERE5D9Fjc586dQ3p6Oo4ePap/3LhxA0qlEl27dsXp06eljtUsbIMmIiJrMyf3iLoPGgC6du2Krl274oknntAvy83NxZEjR5CRkSF2s0RERASJRxILCwvDX/7yF7z22msmvU+j0aBfv37w8vKCn58fxowZg+zsbIN1bt++jfj4ePj4+KBt27YYN24cLl++LGX4REREsiGLoT5TU1MRHx+PQ4cO4dtvv0V1dTVGjBhhcOvWnDlz8Pnnn2PHjh1ITU3FxYsX8fjjj9swaiIiIsuR5XzQV65cgZ+fH1JTUzF48GCUlJSgQ4cO2LZtG/785z8D+OP+6x49euDgwYPo379/g21UVlaisrJS/7y0tBTBwcFsgyYiIqvRtUEXFBQY5B53d3e4u7s3+15Z1KDrKykpAQB4e3sDANLT01FdXY3hw4fr1+nevTs6deqEgwcPNroNjUYDtVqtfwQHB1s+cCIiokYEBwcb5CSNRtPie0R3ErMUrVaL2bNnY+DAgQgPDwcAXLp0CW5ubmjXrp3Buh07dmxygJSkpCQkJCTon+tq0ERERNbWWA26JWYl6JSUFKSkpKCoqAhardbgtY0bN4raZnx8PLKysvDjjz+aE5pRlw+IiIisQaVSWe82q8WLF2PJkiW4++67ERAQAIVCIXZTetOnT8cXX3yBAwcOICgoSL/c398fVVVVuHHjhkEt+vLly/D39zf7c4mIiORGdILesGEDNm/ejCeffNLsIARBwIwZM7Br1y7s378fYWFhBq9HRUXB1dUVKSkpGDduHAAgOzsb+fn5GDBggNmfT0REJDeiE3RVVRXuvfdeSYKIj4/Htm3b8Omnn8LLy0vfrqxWq9G6dWuo1WpMnToVCQkJ8Pb2hkqlwowZMzBgwIBGe3ATERHZO9G9uKdNm4Zt27ZJEsT69etRUlKCIUOGICAgQP/46KOP9OusWrUKDz/8MMaNG4fBgwfD398fO3fulOTziYiI5Eb0fdCzZs3Ce++9h969e6N3795wdXU1eH3lypWSBCgVjsVNlqLValFaWgqVSgWlUpZ3LhKRjdhkLO7jx4+jT58+AICsrCyD16ToMEZkL7RaLW7cuIG2bdsyQRORZEQn6H379kkZBxEREdXB030iIiIZMmugkhs3buDdd9/Vz/3cs2dPTJ06FWq1WpLgiIiInJXoGvSRI0fQpUsXrFq1CsXFxSguLsaqVavQpUsXHD16VMoYiYiInI7oGvScOXPw6KOP4u2334aLyx+bqampwbRp0zB79mwcOHBAsiCJiIicjegEfeTIEYPkDAAuLi6YP38+7r77bkmCIyIiclaiL3GrVCrk5+c3WF5QUAAvLy+zgiIiInJ2ohP0+PHjMXXqVHz00UcoKChAQUEBtm/fjmnTpmHChAlSxkhEROR0RF/ifv3116FQKDB58mTU1NQAAFxdXfHss88iOTlZsgCJiIickeihPnVu3bqF8+fPAwC6dOkCT09PSQKTGof6JEupqanBb7/9hqCgIIM+GURENhnqU8fT0xO9evUydzNERERUh0kJOiEhAUuXLkWbNm2QkJDQ7LpymyyDiIjInpiUoDMyMlBdXa3/f1M4WQYREZF5TErQdSfI2LJlC4KCghrM3iMIAgoKCqSJjoioHk7vSc5C9K87LCwMV69ebbC8uLgYYWFhZgVFRNQU3fSeWq3W1qEQWZToBN1U5++ysjJ4eHiIDoiIiIhE9OLWdQ5TKBRYuHChwW1VtbW1SEtLQ58+fSQLkIiIyBmZnKB1ncMEQcCJEyfg5uamf83NzQ0RERGYO3eudBESERE5IZMTtK6j2FNPPYXVq1dz0A8iIiILED1QyaZNm6SMg4iIiOoQ3UlMo9Fg48aNDZZv3LgRy5YtMysoIiIiZyc6Qf/nP/9B9+7dGyy/6667sGHDBrOCIiIicnaiE/SlS5cQEBDQYHmHDh1QWFhoVlBERETOTnSCDg4Oxk8//dRg+U8//YTAwECzgiIiInJ2ojuJPf3005g9ezaqq6sxdOhQAEBKSgrmz5+P559/XrIAiYiInJHoBD1v3jxcu3YNzz33HKqqqiAIAlq3bo0XXngBiYmJUsZIRETkdEQnaIVCgWXLluGf//wnTp8+jdatW6Nr165wd3eXMj4iIiKnJDpB6+Tn5+PatWuoqqrChQsX9MsfffRRczdNRETktEQn6AsXLmDs2LE4ceIEFAqFfvIM3VzQtbW10kRIRETkhET34p41axbCwsJQVFQET09PZGVl4cCBA7j77ruxf/9+CUMkIiJyPqJr0AcPHsT3338PX19fKJVKtGrVCvfddx80Gg1mzpypn1SDiIiITCe6Bl1bWwsvLy8AgK+vLy5evAgACAkJQXZ2tjTREREROSnRNejw8HAcO3YMYWFhiI6OxvLly+Hm5oa33noLnTt3ljJGIiIipyM6QS9YsADl5eUAgCVLluDhhx/GoEGD4OPjg48++kiyAImIiJyR6AQ9cuRI/f/vvPNOnDlzBsXFxWjfvr2+JzcRERGJI6oNurq6GsOGDcO5c+cMlnt7ezM5ExERSUBUgnZ1dcXx48eljoWIiIj+n+he3JMmTcK7774rZSxERET0/0S3QdfU1GDjxo347rvvEBUVhTZt2hi8vnLlSqO3deDAAaxYsQLp6ekoLCzErl27MGbMGP3rcXFx2LJli8F7Ro4ciT179ogNn4iISNZEJ+isrCz07dsXAHD27FmD10xthy4vL0dERASmTJmCxx9/vNF1Ro0ahU2bNumfc1IOIiJyZKIT9L59+yQLYvTo0Rg9enSz67i7u8Pf31+yzySiP2i1WpSWlkKlUkGpFN3qRUQSM/mv8cKFC/qJMaxp//798PPzQ7du3fDss8/i2rVrza5fWVmJ0tJSgwcRNaTVanHjxg1otVpbh0LksOrno8rKyhbfY3KC7tq1K65cuaJ/Pn78eFy+fNnUzZhk1KhReO+995CSkoJly5YhNTUVo0ePbnbGLI1GA7VarX8EBwdbNEYiIqKmBAcHG+QkjUbT4ntMvsRdv/b81VdfGfVB5njiiSf0/+/Vqxd69+6NLl26YP/+/Rg2bFij70lKSkJCQoL+eWlpKZM0ERHZREFBAVQqlf65Mf2oRLdB21Lnzp3h6+uLnJycJhO0u7s7O5IRSYzt1UTiqFQqgwRtDJP/whQKRYNe2tYePey3337DtWvXEBAQYNXPJXJ2bK8msh5Rl7jj4uL0tdPbt2/jH//4R4P7oHfu3Gn0NsvKypCTk6N/npubi8zMTHh7e8Pb2xuLFy/GuHHj4O/vj/Pnz2P+/Pm48847DcYDJ7KVjIwMzJw5E2vWrEG/fv1sHQ4ROQiTE3RsbKzB80mTJpkdxJEjR3D//ffrn+vajmNjY7F+/XocP34cW7ZswY0bNxAYGIgRI0Zg6dKlvIRNsjBv3jzk5uZi3rx52L9/v63DISIHYXKCrjtYiFSGDBnS7K1b33zzjeSfSSSVFStWYObMmVixYoWtQyEiB2KXncSI5CQyMhIffvghgoKCbB0KETkQdsMkIiKSISZoIiIiGWKCJiIikiEmaCIiIhkyq5NYSkoKUlJSUFRU1GDggo0bN5oVGBFRY3jfOTkL0TXoxYsXY8SIEUhJScHVq1dx/fp1gwcRkSXUve+cyJGJrkFv2LABmzdvxpNPPillPEREzeJ95+QsRCfoqqoq3HvvvVLGQkTUIt53Ts5C9CXuadOmYdu2bVLGQkRERP9PdA369u3beOutt/Ddd9+hd+/ecHV1NXh95cqVZgdHRETkrEQn6OPHj6NPnz4AgKysLIPXrD39JBERkaMRnaD37dsnZRxERERUBwcqISIikiGzBiq5ceMG3n33XZw+fRoA0LNnT0ydOhVqtVqS4IiIiJyV6Br0kSNH0KVLF6xatQrFxcUoLi7GqlWr0KVLFxw9elTKGImIiJyO6AQ9Z84cPProo8jLy8POnTuxc+dO5Obm4uGHH8bs2bMlDJGI5CIjIwMTJkxARkaGrUMhcnhm1aBfeOEFuLj87yq5i4sL5s+fjyNHjkgSHJEj02q1uHHjRoNx7OWMw2wSWY/oBK1SqZCfn99geUFBAby8vMwKisgZ2GOCXrFiBcLCwjjMJpEViE7Q48ePx9SpU/HRRx+hoKAABQUF2L59O6ZNm4YJEyZIGSORQ5BrQjblsrVumM3IyEgrREbk3ET34n799dehUCgwefJk1NTUAABcXV3x7LPPIjk5WbIAiRyFLkG3bdsWSqV87nCse9l6//79tg6HiP6f6KOEm5sbVq9ejevXryMzMxOZmZn6ntzu7u5SxkhEFsTL1kTyZNZ90ADg6emJXr16SRELEdkAZ4cyjlarRWlpKVQqlegrIFJsg5yHSQk6ISEBS5cuRZs2bZCQkNDsupwsg6hpcm2PlhO5JTMpmijk2sxB8mRSgs7IyEB1dbX+/03hZBlEzdNqtSgpKbF1GLLGZEbOzqQEXXeCjC1btiAoKKjBH44gCCgoKJAmOiIiIicl+rQ0LCwMV69ebbC8uLgYYWFhZgVFRETk7EQnaEEQGl1eVlYGDw8P0QEREUnNUm3+7EtAlmRyL25d5zCFQoGFCxfC09NT/1ptbS3S0tLQp08fyQIkIjKXpdqz2U5OlmRygtZ1DhMEASdOnICbm5v+NTc3N0RERGDu3LnSRUhEFqXValFWVsZaIJHMmJygdR3FnnrqKaxevRoqlUryoIjIepigieRJ9EAlmzZtAgCcOnUK+fn5qKqqMnj90UcfNS8yIiIiJyY6Qefm5mLMmDE4ceIEFAqFvtOY7h7o2tpaaSIkIiJyQqJ7NcycORNhYWEoKiqCp6cnTp48iQMHDuDuu+/mgPtERERmEl2DPnjwIL7//nv4+vpCqVRCqVTivvvug0ajwcyZM42auo6IiIgaJ7oGXVtbCy8vLwCAr68vLl68CAAICQlBdna2NNERObCsrCxMnDiRJ7MEgPdUU0OiE3R4eDiOHTsGAIiOjsby5cvx008/YcmSJejcubNkARI5Ko1Go5+HmYgJmuoTnaAXLFig/yEtWbIEubm5GDRoEL766iusWbNGsgCJHFVSUhLnYSaiJolugx45cqT+/3feeSfOnDmD4uJitG/fnrNZERkhPDwcW7duRUhIiK1DIRuQ23SaJD+S/iq8vb1FJecDBw7gkUceQWBgIBQKBXbv3m3wuiAIWLhwIQICAtC6dWsMHz4c586dkyhqIuvIyMjAhAkT7LrNWalUol27dkwoEuAlbWqJ6L+yoUOHYvHixQ2WX79+HUOHDjVpW+Xl5YiIiMC6desafX358uVYs2YNNmzYgLS0NLRp0wYjR47E7du3RcVOZAktHXDnzZunb3POyMjAxIkTkZWVZeUozcMETWQ9ov/K9u/fj3/9618YM2YMysvL9curqqqQmppq0rZGjx6NV155BWPHjm3wmiAIePPNN7FgwQI89thj6N27N9577z1cvHixQU2byJZaStArVqzQtznrkrVGo7FylM6NtVayJ2adBn/33Xe4dOkS+vfvj7y8PIlCMpSbm4tLly5h+PDh+mVqtRrR0dE4ePBgk++rrKxEaWmpwYPIliIjI/Hhhx8iMjJSn6yTkpJsHZZTYYImW6mfjyorK1t8j1kJOiAgAKmpqejVqxf69etnkRHELl26BADo2LGjwfKOHTvqX2uMRqOBWq3WP4KDgyWPjUisyMhIbN26FeHh4bYOxW4x2ZI9CQ4ONshJxlw9E92LW9cZzN3dHdu2bcMrr7yCUaNG4YUXXhC7SUklJSXp564G/jh7YZImchxymYtZF4e3tzfb5qlJBQUFBrM/uru7t/ge0QlaNzmGzoIFC9CjRw/ExsaK3WSj/P39AQCXL19GQECAfvnly5fRp0+fJt/n7u5uVAEQOSrexmMdgiCgpKSEneeoWSqVyuTpmUX/mnJzc+Hr62uwbNy4cTh06BA2btwodrMNhIWFwd/fHykpKfplpaWlSEtLw4ABAyT7HCKx5NqzmZeApeUIt8mRfRF9RAkJCWn0gBQeHm5yLbqsrAyZmZnIzMwE8Efyz8zMRH5+PhQKBWbPno1XXnkFn332GU6cOIHJkycjMDAQY8aMERs+kWTkmqBNodVqUVJSIqtk3lRCtFV5171NzpJ4YkU6Jl3iTkhIwNKlS9GmTRuD9t3GrFy50ujtHjlyBPfff7/B5wBAbGwsNm/ejPnz56O8vBx///vfcePGDdx3333Ys2cPPDw8TAmfyC7Y4tK07jKtj4+PbE406ibEuh1QdQm6pqbGqvGsWLECM2fOtPjQrHJpWyfbMylBZ2RkoLq6Wv//ppg6mtiQIUMatGnX396SJUuwZMkSk7ZLZI94gP6DtRKisXS3yQUFBdk6FHISJiXoffv2AQCqq6uhVCqxYcMGdO3a1SKBEZF1KJVKtG3b1upj6Ld0pcAZEmLdHuBE9Yk6PXd1dcXx48eljoWIbECXoK1dW7dFW6utO3rV/3xd0wLbm6kxov8iJ02ahHfffVfKWIiILCoxMdGmc3Bbq6MZOQbR90HX1NRg48aN+O677xAVFYU2bdoYvG5KJzEick7WHuQjOTkZCQkJkrVrZ2RkYMaMGUhKSjJq2lBz2tV5X7vzEf0tZ2VloW/fvvDy8sLZs2eRkZGhf+hulyKipikUCqjVaqc+2FrrEq/u0jIA/XjoUjB10pO647GbOqMZb79yPqJr0LoOY0TORMpajFzun87IyMCUKVMAABs3bkS/fv1sGo8l6BJpYmIitmzZItl2V6xYoa9Bi41Jo9HgoYcekiwmchzOe+pOJIIj1mLmzZuH7OxsZGdn21XbqCkdvnSzhyUnJ0sagzmTnnBGM2qJ6Bq0zqlTp5Cfn4+qqiqD5Y8++qi5myYiK1ixYoW+Bi2Xe46N0dRAJo3RXVr29/dvdhY8a9Il97oc8QSQxBOdoC9cuICxY8fixIkTUCgU+oFGdPdS1tbWShMhkQOSy+Vt4I9EsXv3bgAwqqOTXMhtIBMpCIKA0tJStGrVytahkAyIPjrMmjULYWFhKCoqgqenJ06ePIkDBw7g7rvvtsi80ERyZ8olVzklaHtVt8OVJbFWS7Yi+uhw8OBBLFmyBL6+vlAqlVAqlbjvvvug0Wgwc+ZMKWMksgst3ePKpGyfrJWgdb26T548afR7ePLg2EQfKWpra+Hl5QUA8PX1xcWLFwH8cYksOztbmuiI7Iiu009Tl1yZoKk5uhM8UzqyMUE7NtFHivDwcBw7dgwAEB0djeXLl+Onn37CkiVL0LlzZ8kCJLIX1rrkSo4jKysLEydOREZGhv4ELzEx0dZhkUyITtALFizQn7UtWbIEubm5GDRoEL766iusWbNGsgCJyHE4Yo2vbpI1lUaj0TeL6Hp133XXXRaIkuyR6AQ9cuRIPP744wCAO++8E2fOnMHVq1dRVFSEoUOHShYgETkOWydoSzQz1E2ypkpKSmq2WYScm6SNYd7e3lafso6IyFiWSNDmJNnw8HBs3brVqGYRW8/ERdZn0n3QCQkJRq/LyTKI7IdCobDrSRhsOZGELsnq7iG3VCymDMxCjsGkBG3smRtr0US2l5GRgZkzZ2LNmjUtjq9tix7mUo4Brrt0bot5ra0ViyMOzELNMylBc4IMIvsh9xqXbgxwAJg7dy4+/fRTu67FS+XkyZNYtmwZ1q5da3DSortLICgoyIbRkTWZ9Zdw48YNvPHGG5g2bRqmTZuGVatWoaSkRKrYiJyCpTpOtXRftq2tWLEC3bp1Q7du3ZCcnGzRzmNarRZlZWV20Xs8OTlZdKczciyiE/SRI0fQpUsXrFq1CsXFxSguLsbKlSvRpUsXHD16VMoYiazCVj2M63+uVJ2B5Hhfdt19040Bvnv3blExmnJZvqamBleuXEFVVZVJ37EtfhOJiYmyPrEi6xGdoOfMmYNHH30UeXl52LlzJ3bu3Inc3Fw8/PDDmD17toQhElmHrW8B0mlpyFB7JuW+mZKgtVotKioqUFNTY3KCLikp0U8GZA133XWX0T27ybGZVYN+4YUX4OLyv2ZsFxcXzJ8/H0eOHJEkOCJnJPdL0+YwZ99MPYGSywmXKbKysvDkk08iKyvL1qGQDIhO0CqVCvn5+Q2WFxQU6MfoJpKaPR50TSXHS9NSMWffbJ2gpe7prlQqoVarDe560Wg0yMvLg0ajcYrfOjVP9C9t/PjxmDp1Kj766CMUFBSgoKAA27dvx7Rp0zBhwgQpYyTSs/VByx4Hi7B1mYklt8lFLJGg628vKSkJoaGhSEpKstvvjaRj0m1Wdb3++utQKBSYPHkyampqAACurq549tlnTZqNhcieGHvrki0HzmgsFrncI2wKXQJzJuHh4Xj//ffRqlUrW4dCMiD6r9XNzQ2rV6/G9evXkZmZiczMTBQXF2PVqlVwd3eXMkYi2TC2DdURaz+OuE+2JrerBLbC31bjTP5VaLVaLFu2DAMHDkS/fv2wZMkS3HnnnejVqxc8PT0tESORbDhy+3BLnPUg2tJsVebcY80E/Qdn/W21xORfxauvvooXX3wRbdu2xR133IHVq1cjPj7eErERERnNUgf5lmarEgRBskFQ7GlAFbI8kxP0e++9h3//+9/45ptvsHv3bnz++efYunUrf1BE1KKmaoxSdL4TBAElJSWSHosyMjJQXl4OPz8/o28NM2dfdMm+sfuuGys7e+y0SMYzOUHn5+fjwQcf1D8fPnw4FAoFLl68KGlgRPaGlytb1lQZJSYmtjiAiS2S0bx581BUVIQ2bdoY3axhqYFmGis7Rx7UhkQk6JqaGnh4eBgsc3V1RXV1tWRBEdkjOSRoe23LS05ObrHzXXPJKCMjAxMnTjQY4EOKhK7rFJiUlGTye6wx0IwjD2pDIhK0IAiIi4vD448/rn/cvn0b//jHPwyWETkzW116tNcEbUznu+aSkS55azSaBsvmzZuHzMxMJCYmYseOHSZ9L5GRkdi6dSvCw8MbfV2pVEKlUhkMNtLcvphzEtfYd+vMnRadgcm/ktjYWPj5+UGtVusfkyZNQmBgoMEyImdmyqVHW7cjGps0LBGnQqGAWq1u8rPrJqXmklFjNd26Cf3FF1/ExYsXJb8krFQq4e3tbfT97lInaHJsJg9UsmnTJkvEQeRQVqxYgZkzZxp16dHW8zYbOyCIJeJs6bPrDrLSHF1Nt/4y3fzJy5Ytw/Tp0xEfH4+3335b0kvCSqXSIoPAZGRkICEhAWvWrDGYF5qcB3uzkF2xdW3TWKZcerR0O6JUZWbP7Z2RkZFYt24d/vrXvzb6vcixdmpMxzlybEzQZFccsdeqpdsRpTrQ23N7Z0ZGBuLj45GZmdno63JM0MZ0nCPHxgRNdsWStThdTTM9PV12B2uxtFotXnrpJYSGhsruQG/NXu+JiYnIz8/Hiy++aPHPkoo9nxCRNJigya5Y8qClq53Pnz/fYRK0IAgIDQ3F1q1bZXegt2aCTk5ORqdOnfDaa69Z/LNMpau9A2gw/SQ1TY5XPaTGBE30/3S188ZmY3OGg4Ej07VB9+nTR/Jtt9TG39Jvp26Cbt++fYPbtpyBmH4SWq0WxcXFKC4udti/S7tI0C+//DIUCoXBo3v37rYOixyIVqtFWFhYkzVNUxO0rQYtsZdOdMaquz9yGAimMS31izDlt3Ps2DE888wzKCgokN1+WpLYviWWGN5VTuzmF3DXXXehsLBQ//jxxx9tHRI5EK1Wi5KSkkbHQBbDVsmk7oFOqVTazSXTpsqr/v40V6YtzTolNq6WyrBuvwhzr7TMmzcPeXl5WLJkiUXGK5cre75DwJLsJkG7uLjA399f//D19W12/crKSpSWlho8iBxd3QOd1CcJdbfXXCISk6SaitWUA3f9WadM2f+mkp8x26jbL8LcEz1jRktzpDsYdJyhQ1z9fFRZWdnie+wmQZ87dw6BgYHo3LkzJk6ciPz8/GbX12g0BiObBQcHWylSItux5IHOkgm6KabsT1JSkkFyMyVBWzv5NVVGxoyW1tzJCvtKyFdwcLBBTqo7LG1T7CJBR0dHY/PmzdizZw/Wr1+P3NxcDBo0CDdv3mzyPUlJSSgpKdE/CgoKrBgx2TNbX0qUa1ur3IWHh4vurS7lJdaWhi8FxCVSY05WmKDlq6CgwCAnGTMBi8lDfdrC6NGj9f/v3bs3oqOjERISgo8//hhTp05t9D3u7u5wd3e3VojkQOoO7CHnoTedkVarRWlpqdFjXxur7rCgzX32jRs3WvxsnmC1zFLfo5ypVCqoVCqT3mOXJdOuXTv86U9/Qk5Ojq1DIQfU3AhO9nbwbWwaRntmbg2Rk1X8jy33x9HK0lLs4yhTT1lZGc6fP4+AgABbh0J2wpQDgqWmC7Qm3f7OnTsXubm5+k5jzko3oYWLi4vsvj9zOqiZg0lS/uTzK23G3LlzkZqairy8PPz8888YO3YsWrVqhQkTJtg6NLIT5h6MLNkubYkDsW5/k5OT0blzZ6xatUrS7du6nb6+lm6HMmfGKUsnyqY6qIn5XLl9L2Qeu0jQv/32GyZMmIBu3brhr3/9K3x8fHDo0CF06NDB1qGRlcnh/mKpWXKfLNWru7nysEWSaKkMtVotysrKRJ2gSXG7VnOk7KDmyLdiNcYS977LiV0k6O3bt+PixYuorKzEb7/9hu3bt6NLly62DotswFYJmgMpGLLW/br1v2+xyb+lBG3sLWQtaWzfW4rZnJOo+uXjTL/TjIwMPPnkkzh79qzDnpDYRYImsjV7GUjBWicwxt6vK3WnLkvVEKVK0I0lyOZilrp87OV3KoV58+bB1dUVNTU1DntCwgRNBMPezvbSEawxcoi9/shaUnZEWrZsGTp16oRly5ZJsj2pNZYgm6vVmjLymNTfrS3bq6X47BUrVqBr1654//33HfaExP6OQEQitHRA0NVyNBqNLJKcuRy1s1BERATWrFmDiIgIq3+2sWVa/6SksaQt5vuR+ndpy/ZqKT47MjISW7duRXh4uISRyYv9HoHIKYmtkbV0QNDVcowZ3cce6PY3MTHR1qFIRlfbFNvRS2wvbh1je1sbM8NSU9+PUqmESqVCaWmpJFcdmvt7sWV7df3PdoSTYktgaZBdEZugWzoYWfts3NL3oDY3t7W5bHX/rDkTUUiRoJv6DYlJLrptxcXFGQwkY80Ebcv26vqfzQTdOJYGOQW5dZ6xdJKz5P6aMkyjI11ql7JMddvavHmzvmlFCkx0joXfIknKnkcnMmaSg8bIeZ8tdcA2drtSt3MqFIoma8JS1z7FMuXeXN2wsq+99pok3xETtGPht0iSknOyaonYg5sU+2ypcrPFAbvuZ0rdztncpermErQpNXlza/3156VuTmRkJLZv347Bgwc3uNc7PT3dbv+WSBpM0GRXbHXJlOMi/09L30HdspJyikRzRgMzpSZvTq1fqVTitddeM+ukRPf58+fPN/s34UhNDM6ICZrsitiDp7kJ1t7uQW0qXilOBKS+bG1qghbTScyUmrw5tX6lUonBgwdj+/btotuqpezg52xDfzoaJmiSlDlDMRpzkBZ78LRkDVjMPlviwFm3DC2ZoM1JYKZ8fmPrnj17Fk8++aTJvy9TOniZ2xnMnJq+FJ9fV3PflS3bq6X6bLH9RuyFY+4V2YzYxGPsgVtuvbEBcftc/8ApRY3aWpfJzfkOzE3Qa9euRV5enqxrhPUTdHPJqLnXmvtNSPH34ggJ2tE7xTnmXpHNONNg/Tpi9rn+gdNRL0VKeQDNyMjArVu30KFDByxbtgzFxcUoLi5uNklJedJi7L7U78gmNkEnJiZabAxvsg9M0CQpsbUre+7MIkWtXsoJJqTUWCymJF0pE3RiYiKuXbuGNm3aICIiAtevX8f169eNTtBST0zR1OeYOyiK7nOWL1/udCe7ZIgJ2snI6eBfl6PWII1lyQkm6jL1RMjcBG3MttPT042KKTk5GcHBwc0OX9rc/hkzBKdYpn5nxtSqo6KiZNecIyeOfnkbYIJ2OnJN0MZeJrb1wBtyZkzyldOJkFarRXFxMRISEnDhwgVJYpLD/hnzPTjC783WnKEMHXfPyK4Ye5lY7gna2BMgS+yHMcnJ1n0E6icvQRAwd+7cRmOqv25iYiIKCgqQnJxsMD1oXc11vjNlhC9z9ksOJwk6zpDEHJrgJEpKSgQAQklJia1DsanDhw8L/fv3Fw4fPmyR7dfW1grXr18XamtrTXpfdXW1kJubK1RXV1skLkuqu89i96PuNsR+R8a8z9T4pP5eYmJihI4dOwoxMTFCdXW1kJOTI2RnZws5OTkNPqPuuoIgCD///LMQGRkp7Ny5Uxg8eLDg5+cnREdHG7yvfry6bQwePFiIjo4W/Pz89NtripjfcExMjODn5ycMHDhQOHTokCR/Y83FYem/Y5KOObmHCdoBNfeHXf+gJxf2nKDrkmI/xH5Hxnx2S+vU/+1I/b3UTSyHDx8WoqOjhZ07dzaaoOsnoYqKCuGHH34Q0tPThZ9++kkIDw8XwsPDhZSUlCbj1W3j0KFDQkpKihAdHW2RpKbbl927dwsVFRUW/y3L9e+YGmKCNoIzJejmDqpyPfMWW/OWGykSmtjvSIoE3VSCk+r3Unf7MTExQvv27QUvLy9h9+7dDWKqH0tFRYVw4MABITU1Vbh586a+Rly3Fl3/PXWfW/IkUHc1ICcnxyoJWq5/x9SQObmHDRMOqLlOKnIc6AOw77YyqTveyek7smR76ooVK1BTUwMXFxdR0y0mJSUhNDQUSUlJkscmd3L6jZDl2N/RkFokp04qzkCuPePNodunZcuWtdipTOz+R0ZG4v3330fXrl2NTrJ1p5sMDw/H+++/j/Dw8CbXr3viZ+mTQEt2QiPnxATtgOQ6/q4zkGLAFUt+R8ZuW5d0IyIiJJuNqjHh4eH44IMPMHDgQKP299y5c3j22Wdx4sQJqNVqKBSKZte3VoLWTa2Zm5vb7H3aRKbgUdoByXX8XUcl9S02Yr8jY97X0jq2GNHN2P11cXHBf/7zH+Tl5eGFF15Au3btcOrUKVnUWpVKJVatWoXOnTtj+fLl/BsjSfAXRGSmuknZlvcZS3HypdsXU2qBpszeZMq69ffHxcUFa9euRefOnfXlq9FoDE6IbHkCqjsxjoqKYoImSfAX5IBYS7auuknZ3jvv6PYlLi7O6Jp0eno6YmNjkZ6ebtJnKZVKtG/fHu3btzd6Ion65ZuUlGRwQuQsv31n2U9npxAEEbOf26HS0lKo1WqUlJRApVLZOhyL0Gq1KC0thUql4h+uFdXU1OC3335DUFAQAOj/7+LiYuPITKfbl9jYWGRnZ6Nbt27YsmVLs/sTExODU6dOoWfPnkhNTTVq+2LLp+77lUoliouLUVpaik6dOtm8vPn3R40xJ/fwV+RAHLE3MdlGcnIywsLCjGpPjYuLw+3btxEXF9fids1t47Zmr2xTyS0esn/8JRGZyREPzKa0p27evBkeHh7YvHlzi9s1txOdI5Y1UVP4K3cg9jynsj2Tc63OGpKTk9GpUyckJye3uK6tJ+sgsidsg3YgQ4YMwZkzZ9C9e3fs37/f1uGQHdK1o7Zt2xZlZWVGtafevn0bmZmZ6NOnDzw8PJpd19w26KbiZbsvyRXboAkAaydkPt0VABcXF6OvBLi4uCAwMNAmnbSc8YoFOQ/762ZKTdK1G+p6ExNZg4uLCzp16mTrMIgcDk87ichqWOMlMh5r0A6EBz+SO91vlIhaxgTtQHjwIyJyHKxqERERyRATNBERkQwxQRMREcmQXSXodevWITQ0FB4eHoiOjsbhw4dtHRIREZFF2E2C/uijj5CQkIBFixbh6NGjiIiIwMiRI1FUVGTr0IiIiCRnNwl65cqVePrpp/HUU0+hZ8+e2LBhAzw9PbFx40Zbh0ZERCQ5u7jNqqqqCunp6UhKStIvUyqVGD58OA4ePNjoeyorK1FZWal/XlJSAuCPcVGJiIisQZdzdDlIx93dHe7u7s2+1y4S9NWrV1FbW4uOHTsaLO/YsSPOnDnT6Hs0Gg0WL17cYHlwcLBFYiQiImpK/eFwFy1ahJdffrnZ99hFghYjKSkJCQkJ+udarRbFxcXw8fGBQqGwYWRNKy0tRXBwMAoKChx2xi1rYVlKh2UpHZaltOyhPAVBwLVr1+Dt7W0wymNLtWfAThK0r68vWrVqhcuXLxssv3z5Mvz9/Rt9T2OXD+xllC2VSiXbH5u9YVlKh2UpHZaltORenmq1WtT77KKTmJubG6KiopCSkqJfptVqkZKSggEDBtgwMiIiIsuwixo0ACQkJCA2NhZ333037rnnHrz55psoLy/HU089ZevQiIiIJGc3CXr8+PG4cuUKFi5ciEuXLqFPnz7Ys2dPg45j9szd3R2LFi0yqm2CmseylA7LUjosS2k5enkqBEEQbB0EERERGbKLNmgiIiJnwwRNREQkQ0zQREREMsQETUREJENM0DZw4MABPPLIIwgMDIRCocDu3bsNXi8rK8P06dMRFBSE1q1b6ycHoYY0Gg369esHLy8v+Pn5YcyYMcjOzjZY5/bt24iPj4ePjw/atm2LcePGNRj0hlouy+LiYsyYMQPdunVD69at0alTJ8ycObPBGMNk3O9SRxAEjB49utFjARlflgcPHsTQoUPRpk0bqFQqDB48GBUVFTaIWDpM0DZQXl6OiIgIrFu3rtHXExISsGfPHnzwwQc4ffo0Zs+ejenTp+Ozzz6zcqTyl5qaivj4eBw6dAjffvstqqurMWLECJSXl+vXmTNnDj7//HPs2LEDqampuHjxIh5//HEbRi1PLZXlxYsXcfHiRbz++uvIysrC5s2bsWfPHkydOtXGkcuPMb9LnTfffFO2ww/LgTFlefDgQYwaNQojRozA4cOH8csvv2D69OkGQ2vaJYFsCoCwa9cug2V33XWXsGTJEoNlffv2FV566SUrRmafioqKBABCamqqIAiCcOPGDcHV1VXYsWOHfp3Tp08LAISDBw/aKky7UL8sG/Pxxx8Lbm5uQnV1tRUjsz9NlWVGRoZwxx13CIWFhY0eC6ihxsoyOjpaWLBggQ2jsgw7P71wTPfeey8+++wz/P777xAEAfv27cPZs2cxYsQIW4cme7rLrd7e3gCA9PR0VFdXY/jw4fp1unfvjk6dOjU5VSn9oX5ZNrWOSqWCi4vdjHlkE42V5a1bt/C3v/0N69ata3JOAWqoflkWFRUhLS0Nfn5+uPfee9GxY0fExMTgxx9/tGWYkmCClqG1a9eiZ8+eCAoKgpubG0aNGoV169Zh8ODBtg5N1rRaLWbPno2BAwciPDwcAHDp0iW4ubk1mCilY8eOuHTpkg2itA+NlWV9V69exdKlS/H3v//dytHZl6bKcs6cObj33nvx2GOP2TA6+9JYWV64cAEA8PLLL+Ppp5/Gnj170LdvXwwbNgznzp2zZbhm42mvDK1duxaHDh3CZ599hpCQEBw4cADx8fEIDAw0qAmSofj4eGRlZTnEmbOttVSWpaWleOihh9CzZ88W57R1do2V5WeffYbvv/8eGRkZNozM/jRWllqtFgDwzDPP6OdmiIyMREpKCjZu3AiNRmOTWKXABC0zFRUVePHFF7Fr1y489NBDAIDevXsjMzMTr7/+OhN0E6ZPn44vvvgCBw4cQFBQkH65v78/qqqqcOPGDYNadHNTlTq7pspS5+bNmxg1ahS8vLywa9cuuLq62iBK+9BUWX7//fc4f/58gys748aNw6BBg7B//37rBmoHmirLgIAAAEDPnj0N1u/Rowfy8/OtGqPUeIlbZqqrq1FdXd2g92GrVq30Z4r0P4IgYPr06di1axe+//57hIWFGbweFRUFV1dXg6lKs7OzkZ+fz6lK62mpLIE/as4jRoyAm5sbPvvsM3h4eNggUvlrqSwTExNx/PhxZGZm6h8AsGrVKmzatMkGEctXS2UZGhqKwMDABrdenT17FiEhIdYMVXKsQdtAWVkZcnJy9M9zc3ORmZkJb29vdOrUCTExMZg3bx5at26NkJAQpKam4r333sPKlSttGLU8xcfHY9u2bfj000/h5eWlb1dWq9Vo3bo11Go1pk6dioSEBHh7e0OlUmHGjBkYMGAA+vfvb+Po5aWlstQl51u3buGDDz5AaWkpSktLAQAdOnRAq1atbBm+rLRUlv7+/o1ewenUqVOjJ0bOrKWyVCgUmDdvHhYtWoSIiAj06dMHW7ZswZkzZ/DJJ5/YOHoz2bQPuZPat2+fAKDBIzY2VhAEQSgsLBTi4uKEwMBAwcPDQ+jWrZvwxhtvCFqt1raBy1Bj5QhA2LRpk36diooK4bnnnhPat28veHp6CmPHjhUKCwttF7RMtVSWTf1uAQi5ubk2jV1ujPldNvYe3mbVkLFlqdFohKCgIMHT01MYMGCA8MMPP9gmYAlxukkiIiIZYhs0ERGRDDFBExERyRATNBERkQwxQRMREckQEzQREZEMMUETERHJEBM0ERGRDDFBExERyRATNBERkQwxQRMREckQEzQRkYmuXbsGPz8/5OXlmfS+J554Am+88YZlgiKHwwRNDisuLg4KhQIKhQKurq7o2LEjHnjgAWzcuFHSqTuHDBmC2bNnN7uOLo6mHi+//LJk8ciNMeVjSd98802L5b93716Ttvnqq6/iscceQ2hoqH5ZXFwcxowZY7DeJ598Ag8PD31SXrBgAV599VWUlJSYu1vkBJigyaGNGjUKhYWFyMvLw9dff437778fs2bNwsMPP4yamhqrxVFYWKh/vPnmm1CpVAbL5s6da7VYpFJVVWUXnzd48GCDsvbx8cE///lPg2XDhg0zenu3bt3Cu+++i6lTpza73jvvvIOJEydi/fr1eP755wEA4eHh6NKlCz744ANR+0JOxtbTaRFZSmxsrPDYY481WJ6SkiIAEN5++21BEAQhJCREWLVqlcE6ERERwqJFi/TPd+zYIYSHhwseHh6Ct7e3MGzYMKGsrEyIjY01eerFTZs2CWq1utHXamtrhddee00IDQ0VPDw8hN69ews7duzQvx4TEyNMnz5dmDVrltCuXTvBz89PeOutt4SysjIhLi5OaNu2rdClSxfhq6++MthuTEyMEB8fL8THxwsqlUrw8fERFixYoJ/CtKXPrbuNWbNmCT4+PsKQIUMEQRCEr7/+Whg4cKCgVqsFb29v4aGHHhJycnL030FT5dNSuTf1ecbE2pTffvtNACDs3bvXqPUbs2PHDqFDhw4Nltf9vS1btkzw8PAQdu7c2WC9xYsXC/fdd5/ozyfnwRo0OZ2hQ4ciIiICO3fuNGr9wsJCTJgwAVOmTMHp06exf/9+PP744xAEAatXr8aAAQPw9NNP62tjwcHBomPTaDR47733sGHDBpw8eRJz5szBpEmTkJqaql9ny5Yt8PX1xeHDhzFjxgw8++yz+Mtf/oJ7770XR48exYgRI/Dkk0/i1q1bBtvesmULXFxccPjwYaxevRorV67EO++8Y/Tn6rbh5uaGn376CRs2bAAAlJeXIyEhAUeOHEFKSgqUSiXGjh0LrVZrdvk09nnGxtqYjIwMAEDfvn2NjqG+H374AVFRUU2+/sILL2Dp0qX44osvMHbs2Aav33PPPTh8+DAqKytFx0BOwtZnCESW0lQNWhAEYfz48UKPHj0EQWi5Jpeeni4AEPLy8hrdVkxMjDBr1iyj42qqBn379m3B09NT+Pnnnw2WT506VZgwYYL+s+rWvmpqaoQ2bdoITz75pH5ZYWGhAEA4ePCgQYw9evTQ15gFQRBeeOEFoUePHkZ9rm4bkZGRLe7flStXBADCiRMn9O9rrHyMqUHX/zxjY23K4sWLheDg4BbXGzZsmHD27NlGX3vssceEKVOmNFgeGxsruLm5CQCElJSUJrd97NixZn9PRDouNj07ILIRQRCgUCiMWjciIgLDhg1Dr169MHLkSIwYMQJ//vOf0b59e0ljysnJwa1bt/DAAw8YLK+qqkJkZKT+ee/evfX/b9WqFXx8fNCrVy/9so4dOwIAioqKDLbTv39/g30eMGAA3njjDaM/F0CjNcdz585h4cKFSEtLw9WrV/Ud8PLz8xEeHm7Uvjel/ueZEmtjjh49alTt+dy5c+jSpUujr1VUVMDDw6PR13r37o2rV69i0aJFuOeee9C2bdsG67Ru3RoAGlzhIKqPCZqc0unTpxEWFgYAUCqVEATB4PXq6mr9/1u1aoVvv/0WP//8M/bu3Yu1a9fipZdeQlpamn4bUigrKwMAfPnll7jjjjsMXnN3d9f/39XV1eA1XS/1us8BGN1T3djPBYA2bdo0eP8jjzyCkJAQvP322wgMDIRWq0V4eHiLnbpaKvfGPs+UWBtz9OhRTJs2rcHykydPYsqUKaioqMD48ePh7+8PpbLxFkBfX19cv3690dfuuOMOfPLJJ7j//vsxatQofP311/Dy8jJYp7i4GADQoUOHFuMl58YETU7n+++/x4kTJzBnzhwAfxwoCwsL9a+XlpYiNzfX4D0KhQIDBw7EwIEDsXDhQoSEhGDXrl1ISEiAm5sbamtrzY6rZ8+ecHd3R35+PmJiYszeXn1paWkGzw8dOoSuXbua9bnXrl1DdnY23n77bQwaNAgA8OOPPxqs01T5GFPu9ZkT69WrV1FQUNCgBl1ZWYnx48fjww8/RK9evTBmzBiDqxT1RUZGNtsLOyQkBKmpqfokvWfPHoMknZWVhaCgIPj6+poUPzkfJmhyaJWVlbh06RJqa2tx+fJl7NmzBxqNBg8//DAmT54M4I9OY5s3b8YjjzyCdu3aYeHChWjVqpV+G2lpaUhJScGIESPg5+eHtLQ0XLlyBT169AAAhIaGIi0tDXl5eWjbti28vb2brH01x8vLC3PnzsWcOXOg1Wpx3333oaSkBD/99BNUKhViY2PNKov8/HwkJCTgmWeewdGjR7F27Vq88cYbZn1u+/bt4ePjg7feegsBAQHIz89HYmKiwTpNlU9L5S51GR09ehRAww5iu3fvRkxMjL6ZoEePHggMDGxyOyNHjkRSUhKuX7/eZDNHcHAw9u/fj/vvvx8jR47Enj17oFKpAPzRyWzEiBHN7icRwARNDm7Pnj0ICAiAi4sL2rdvj4iICKxZswaxsbH6JJqUlITc3Fw8/PDDUKvVWLp0qUFNTqVS4cCBA3jzzTdRWlqKkJAQvPHGGxg9ejQAYO7cuYiNjUXPnj1RUVGB3NxcgwEsTLF06VJ06NABGo0GFy5cQLt27dC3b1+8+OKLZpfF5MmTUVFRgXvuuQetWrXCrFmz8Pe//92sz1Uqldi+fTtmzpyJ8PBwdOvWDWvWrMGQIUP06zRVPi2Ve1PExpqRkYGOHTs2SL4nTpxAnz599M/T09MxatSoJrfTq1cv9O3bFx9//DGeeeaZJtcLCgoySNLffPMN3NzcsHv3buzZs6fF/SRSCPUbgYjI4QwZMgR9+vTBm2++aetQZGflypX49ddfsXr1auzduxejRo3CtWvXmu0E+OWXX2LevHnIysoy6WrJ+vXrsWvXLpNHLiPnxBo0ETm1SZMmYfTo0YiMjER4eDhCQkJa7KH/0EMP4dy5c/j9999Nuq/b1dUVa9euNTdkchKsQRM5AdagiewPEzQREZEMcahPIiIiGWKCJiIikiEmaCIiIhligiYiIpIhJmgiIiIZYoImIiKSISZoIiIiGWKCJiIikiEmaCIiIhligiYiIpKh/wMTnoBh791HWQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Polarization fraction as a function of the dust column density\n",
    "PvT_plot_Step2 = HtoH.PvT(CatP, CatT, Tscale=[17,27], Pscale=[0.0,30.0], CatdP=CatdP, errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "491d0620",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Saving figures\n",
    "# PNG\n",
    "PvI_plot_Step2.savefig('Plots/'+SourceName+'_PvI_plot_Step2.png',dpi=300)\n",
    "PvN_plot_Step2.savefig('Plots/'+SourceName+'_PvN_plot_Step2.png',dpi=300)\n",
    "PvT_plot_Step2.savefig('Plots/'+SourceName+'_PvT_plot_Step2.png',dpi=300)\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "PoLiteWIP",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
